AD-A092  718 


MICHIGAN  UN IV  ANN  ARBOR  DEPT  OF  NAVAL  ARCHITECTURE  —ETC  F/8  13/10 
ADAPTIVE  PATH  CONTROL  OF  SURFACE  SHIPS  IN  RESTRICTED  WATERS. (U) 

AUG  80  M  &  PARSONS*  H  T  CUONG  D0T-C6-7A843-A 


AD  A  0  9  2  7  1  8 


ftoport  No,  CfHH5-a) 


ADAPTIVE  fWK  OOKTRDL 
OF 

SURFACE  SHIPS 
IN 

RESTRICTED  WATERS 

mom  G.  RARSOfS 
HUA7UCU0NG 


FINAL  REPORT 
AUGUST  1980 

Document  it  available  to  the  public  through  the 
National  Technical  Information  Service, 
Springfield,  Virginia  22161 


Proporod  for 

U.S.  DEPARTMENT  OF  TRANSPORTATION 
United  States  Coast  Guard 
Offlos  of  Roaooreh  and  Dovotopmont 
Washington,  D.C.  20990 


Mcmci 


7y*  k  dtoaamfcw fad  mtdm  th§  s&cmordnip  of  fho  Dcpanmoot 

of  lmwportK*m  in  *m  mmel  of  Sfttormitkw  «*cton«»»  Th*  IMftd 
i&ta«  Gowtrnrnoot  mwiw  rs©  WsUh*  lot  Hi  wrnmm  or  s*»  tfionof. 

Th*  eontmti  of  rmert  do  not  iritiaiaMrfv  Rifteet  th*  crffteW  ntnwr 
or  poiiry  of  *#**  Comc  Ctoard;  md  th*y  H&  mi  aaoMtat*  *  tfmdurd, 
s^M^fka^ksn,  or  ittfutotkm. 

This  fteort,  or  portfom  tmwof  mm  not  te  mors  tor  adrtrdiint  or 
*t*t  prixn&WMi  pwrpam,  Chotk**  of  wadi  mm t  md  iimwtoc'oirw't 
dk*s*  oof  wwmttwi*  or  apptvm  of  *«#•»  prodnet*. 


rcdtrmmf  UffVti itliM  Nomt  and 

Department  of  Naval  Architecture  and 
Marine  Engineering 
The  University  of  Michigan 
Ann  Arbor,  Michigan  48109 


12  SfMttnng  Agency  N«m«  and  Add'ltl 

U.S.  Department  of  Transportation 
U.S.  Coast  Guard 

Office  of  Research  and  Development 
Washington,  D.C.  20593 


15  **•••» 


10.  Unit  Nn.  (TSAIS) 


IwJ  **  / 


DClY-CG-74843-A  h 


•port  and  Period  Coworod 


,  f 

/  Final  i 


14  AbitiMct 

<  * 

The  feasibility  and  effectiveness  of  using  an  integrated  system  for  the 
adaptive  control  of  a  surface  ship  along  a  prescribed  path  in  restricted 
waters  is  investigated.  The  prescribed  path  is  restricted  to  straight  line 
segments  which  would  extend  between  isolated  turning  maneuvers.  The  controller 
consists  of  four  major  components  arranged  in  two  loops:  an  inner  or  control 
loop  and  an  outer  or  gain  update  loop.  The  inner  loop  consists  of  a  Kalman 
state  estimator  and  an  optimal  stochastic  control  law  which  would  provide 
effective  control  when  subjected  to  disturbances  and  measurement  noise.  The 
outer  loop  estimates  the  parameters  of  the  system  equations  of  motion  to  enable 
the  controller  to  adapt  itself  for  the  changes  in  ship  characteristics  which 
take  place  due  to  changes  in  the  ship  operating  condition  and  environment. 

The  Brownian  motion  process  is  shown  to  be  an  effective  model  for  bias  distur¬ 
bances  in  the  design  of  the  Kalman  filter.  A  design  criterion  is  proposed  for 
the  selection  of  the  appropriate  value  of  the  diffusion  coefficient  to  achieve 
the  desired  character  of  the  disturbance  model.  A  control  loop  designed  using 
this  approach  is  shown  to  perform  very  effectively  with  typical  disturbances. 
The  transient  response  of  the  control  loop  is,  however,  sensitive  to  variations 
in  two  of  the  coefficients  in  the  equations  of  motion. V  Weighted  Least-Squares 
(WLS)  and  Minimum  Variance  (MVE)  parameter  estimation  algorithms,  which  might 
be  used  to  estimate  these  coefficients  on-line,  are  investigated.  For  effec¬ 
tive  parameter  estimation,  the  ship  must  be  excited  to  a  sufficient  level  of. 
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motion  using  an  open-loop  rudder  coronand.  Adaptive  path  control  would,  there¬ 
fore,  have  to  consist  of  alternate  periods  of  open-loop  control  and  closed- 
loop  path  control.  The  WLS  algorithm  is  shown  to  be  more  appropriate  for 
off-line  estimation  of  time-invariant  parameters.  The  MVE  algorithm  shows 
promise  as  an  effective  on-line  ship  parameter  estimator  for  use  in  an 
adaptive  system. 
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Abstract 


Itie  feasibility  and  effectiveness  of  using  an  integrated  system  for  the 
adaptive  control  of  a  surface  ship  along  a  prescribed  path  in  restricted 
waters  is  investigated.  The  prescribed  path  is  restricted  to  straight  line 
segments  which  would  extend  between  isolated  turning  maneuvers.  The  controller 
consists  of  four  major  components  arranged  in  two  loops :  an  inner  or  control 
loop  and  an  outer  or  gain  update  loop.  The  inner  loop  consists  of  a  Kalman 
state  estimator  and  an  optimal  stochastic  control  law  which  would  provide 
effective  control  when  subjected  to  disturbances  and  measurement  noise.  The 
outer  loop  estimates  the  parameters  of  the  system  equations  of  motion  to  enable 
the  controller  to  adapt  itself  for  the  changes  in  ship  characteristics  which 
take  place  due  to  changes  in  the  ship  operating  condition  and  environment. 

The  Brownian  motion  process  is  shown  to  be  an  effective  model  for  bias  distur¬ 
bances  in  the  design  of  the  Kalman  filter.  A  design  criterion  is  proposed  for 
the  selection  of  the  appropriate  value  of  the  diffusion  coefficient  to  achieve 
the  desired  character  of  the  disturbance  model.  A  control  loop  designed  using 
this  approach  is  shown  to  perform  very  effectively  with  typical  disturbances. 
The  transient  response  of  the  control  loop  is,  however,  sensitive  to  variations 
in  two  of  the  coefficients  in  the  equations  of  motion.  Weighted  Least-Squares 
(WLS)  and  Minimum  Variance  (MVE)  parameter  estimation  algorithms,  which  might 
be  used  to  estimate  these  coefficients  on-line,  are  investigated.  For  effec¬ 
tive  parameter  estimation,  the  ship  must  be  excited  to  a  sufficient  level  of 
motion  using  an  open-loop  rudder  command.  Adaptive  path  control  would,  there¬ 
fore,  have  to  consist  of  alternate  periods  of  open-loop  control  and  closed- 
loop  path  control.  The  WLS  algorithm  is  shown  to  be  more  appropriate  for 
off-line  estimation  of  time-invariant  parameters.  The  MVE  algorithm  shows 
promise  as  an  effective  on-line  ship  parameter  estimator  for  use  in  am  adaptive 
system. 
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1.  Introduction 


The  problem  of  controlling  surface  ships  alonq  prescribed  paths  in 
maneuvering  situations  is  becoming  increasingly  imtxirtant  from  operational, 
safety,  and  environmental  viewpoints.  In  these  situations,  ships  are 
subjected  to  short-term,  essentially  zero-mean  disturbances  due  to  passinq 
ships,  current  and  wind  variations,  waves,  and  bank  and  bottom  changes.  Ships 
are  also  subject  to  more  long-term,  non-zero-mean  disturbances  due  to  current, 
wind,  second-order  wave  forces,  and  banks.  The  dynamic  characteristics  of  the 
ships  also  change  significantly  depending  on  depth-under-keel,  draft,  trim, 
and  speed.  These  maneuvering  situations  can  place  severe  demands  on  pilots, 
conning  officers,  and  helmsmen  thus  making  some  form  of  automated  control 
desirable  or  perhaps  even  necessary  in  the  future. 

In  this  report,  we  investigate  the  use  of  two  types  of  adaptive  control¬ 
lers  for  the  path  control  of  surface  ships  in  restricted  waters.  In  this 
context,  restricted  waters  include  straits,  channels,  harbors,  canals,  rivers, 
and  harbors  where  the  ship  must  steer  a  prescribed  path  and  not  just  maintain 
a  prescribed  heading.  In  our  earlier  work,1  we  investigated  the  use  of  non- 
adaptive,  optimal  stochastic  controllers  for  this  purpose.  These  control 
systems  consisted  of  a  steady-state  Kalman  filter  and  a  steady-state  optimal 
state  feedback  controller.  The  Kalman  filter  uses  noisy  measurements  to 
qenerate  an  unbiased  estimate  of  the  state  which  is  then  used  by  the  controller 
to  generate  the  rudder  command.  These  controllers  were  shown  to  provide 
effective  control  when  a  ship  is  subjected  to  short-term,  essentially  zero- 
mean  disturbances.  Our  earlier  work  did  not  produce  controllers  which  could 
accommodate  more  long-term  disturbances  without  a  mean  offset  from  the  desired 
path.  We  also  showed  the  desirability  of  an  adaptive  system  which  could 
automatically  account  for  changes  in  the  dynamic  characteristics  of  the  ship 
due  to  changes  in  depth-under-keel  or  other  operating  conditions. 

In  this  work,  we  develop  surface  ship  path  controllers  which  can  accom¬ 
modate  long-term  disturbances  with  a  zero-mean  offset  from  the  path  and  which 
can  adapt  for  changes  in  the  dynamic  characteristics  of  the  ship.  Our  overall 
approach  is  shown  schematically  in  Fig.  1.  The  control  system  consists  of 
four  major  components  which  are  arranged  into  two  loops.  The  inner  or  control 
loop  consists  of  a  steady-state  Kalman  filter  or  state  estimator  and  a  steady- 
state  optimal  state  feedback  controller.  This  inner  loop  is  similar  to  the 


control  systems  used  in  our  earlier  work1  except  that  the  design  used  here 
can  accommodate  long-term  disturbances.  The  outer  or  gain  update  loop  con¬ 
sists  of  an  on-line  parameter  estimator  and  a  second  function  which  recalcu¬ 
lates  the  steady-state  filter  and  controller  gains  using  the  latest  estimates 
of  the  ship  parameters.  The  gain  update  loop  would  be  implemented  batch-wise 
or  at  a  slower  rate  than  the  control  loop.  Thus,  the  control  loop  uses  piece 
wise  constant  gains.  The  control  and  gain  update  loops  would  both  be  imple¬ 
mented  in  an  onboard  digital  computer. 


Figure  1.  Overall  Schematic  of  Adaptive  Fath  Controller 

A  number  of  approaches  can  be  taken  to  develop  an  adaptive  ship  path 
controller.  The  preferred  approach  would  be  a  design  which  would  be  robust 
enough  to  provide  effective  control  under  all  ship  operating  conditions  with¬ 
out  adjustment.  In  general,  this  is  not  feasible.  For  years,  ship  heading 
autopilots  have  included  manual  deadband  and  gain  adjustments  to  accommodate 
changing  sea  states.  This  approach  can  provide  only  limited  adaptation  and 
requires  proper  operator  action.  An  alternative  approach  would  be  to  have 


controller  gains  preprogrammed  as  functions  of  a  few  operating  conditions  such 
as  depth-to-draft  ratio  H/T  and  draft.  This  approach  would  be  feasible  but 
could  not  cover  all  variables  and  would  require  more  extensive  knowledge  of 
ship  characteristics  than  is  usually  available.  The  final  approach,  which 
we  use  here,  is  to  have  the  controller  automatically  establish  the  predomi¬ 
nant  part  of  the  ship  dynamics  on-line  thus  permitting  adaptation  to  changing 
conditions . 

Work  has  been  underway  in  recent  years  on  the  development  of  adaptive 
ship  heading  autopilots.  Development  in  The  Netherlands  began  with  simplified 
versions  of  Bech's  second-order  turninq  rate  (iJj)  equation;  i.e., 

ip  +  —  +  —  ip  +  - -  (ai<3  +  bij;)  =  - -  (T^  +  6)  .  (1) 

L  1  2j  12  12 

Honderd  and  Winkelman2  assumed  that  '3=0  and  that  only  b  varies  with  water 
depth.  They  designed  a  simple  model  reference  adaptive  controller  for  heading 
control  using  a  sensitivity  model  approach,  van  Amerongen  and  Udink  ten 
Cate3*4  developed  a  second  model  reference  adaptive  controller  using  a 
stability  (Liapunov)  approach.  van  Amerongen,  Nieuwenhuis,  and  Udink  ten 
Cate5  also  reported  the  development  of  a  model  reference  adaptive  controller 
using  a  gradient  based  method.  In  this  latter  work,  the  <5  term  in  eq.  (1) 
was  retained  and  both  a  and  b  were  assumed  to  vary  with  depth.  Some  of 

these  heading  autopilots  have  been  successfully  tested  at  sea.  van  Amerongen 
and  van  Nauta  Lemke6  reported  the  most  recent  Dutch  work  on  model  reference 
type  adaptive  heading  autopilots  in  ld78. 

In  contrast  to  the  Dutch  work,  a  number  of  recent  heading  autopilot 
designs  include  adaptive  features  which  are  preprogrammed;  i.e.,  open-loop 
schedules  which  are  prescribed  in  advance.  Oldenbuig7  describes  an  autopilot 
which  adapts  for  speed,  depth-under-keel  and  sea  state  in  this  manner.  Ware, 
Fields,  and  Bozzi8  describe  an  autopilot  which  includes  a  pair  of  parallel 
notch  filters  in  which  the  notch  frequencies  are  scheduled  as  a  function  of 
vessel  speed  in  order  to  provide  adaptation  to  changes  in  wave  encounter 
frequency.  Sugimoto  and  Kojima8  report  another  example  of  a  recent  autopilot 
design  with  preprogrammed  gain  adjustment.  Their  autopilot  is  designed  to 
adapt  for  loading  condition,  speed,  and  sea-state. 


A  third  direction  in  adaptive  headinq  autopilot  development  has  evolved 
from  the  work  of  Astrttm  in  Sweden.  AstriJm1  u  recently  demonstrated  the  need 
for  adaptive  ship  heading  controllers  and  noted  that  Kallstrom  had  used  the 
approach  presented  by  Xstrttm,  et  al1 1 ' to  desiqn  adaptive  headinq  auto¬ 
pilots  for  a  number  of  larqe  tankers.  He  noted  that  these  have  been  tested 
successfully  and  that  one  had  been  operatinq  for  more  than  one  year  at  sea 
(1976).  Kallstrdm  and  his  colleagues13  reported  on  this  work  in  1977.  In 
this  approach,  the  system  is  modeled  as  a  discrete,  sinqle-input ,  single- 
output  system  (auto  regressive  movinq  averaqe  formulation)  disturbed  by  white 
noise.  The  unknown  system  parameters  are  then  estimated  by  recursive  least- 
squares,  extended  least-squares  or  maximum  likelihood  schemes.  A  minimum 
output  variance  control  law  is  used.  Volta  and  Tiano11'  have  also  developed 
adaptive  heading  controllers  based  on  the  work  of  Xstrom  and  Wittenmark . 1 1 
Simulation  of  these  controllers  was  reported  by  Brink,  Baas,  Tiano,  and 
Volta15  in  1978. 


In  this  work,  we  aie  concerned  with  path  control  instead  of  heading  control 
of  the  ship.  Our  approach  also  user;  mu  1 1 i - var iabl e  instead  of  single-input, 
single-output  methods.  In  restricted  waters  where  path  control  is  of 
interest,  the  conditions  which  iltei  the  rhuracter  ist  ios  of  the  ship,  such  as 
water  depth,  can  change  fairly  i  idly  rttp.ired  to  the  dynamics  of  the  ship. 
This  creates  a  difficult  on-  i  in-  ;  a  turn*  t .  r  *•  t  im.it  i  »i ,  ptonlom.  Off-line 
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which  might  be  50  to  100  s.  and  a  ship  length  or  two.  This  makes  practical 
on-line  parameter  estimation  from  noisy  measurements  very  difficult  in 
restricted  waters. 

This  report  is  presented  in  four  principal  parts.  In  Section  2,  we 
formulate  the  ship  path  control  problem  including  a  development  of  the  linear 
equations  of  motion,  a  discussion  of  measurement  selection,  and  a  development 
of  the  design  process  disturbances  which  are  used  in  the  evaluation  of  system 
performance  using  digital  simulation.  In  Section  3,  we  present  the  develop¬ 
ment  and  evaluation  of  the  inner  or  control  loop  part  of  the  path  controller. 
This  presentation  also  includes  the  optimal  state  feedback  controller  and 
Kalman  filter  gain  calculation  methods  used  in  the  gain  update  calculation 
function  of  the  outer  loop  in  Fig.  1.  This  section  includes  a  review  of  the 
shaping  filter  approach  for  process  disturbance  modeling  which  we  used  in  our 
earlier  work1  and  which  cannot  accommodate  more  long-term  process  disturbances 
without  a  non-zero  mean  offset  from  the  path.  We  then  present  an  alternative 
approach  using  a  random-walk  process  disturbance  model  which  can  accommodate 
the  more  long-term  disturbances.  This  section  closes  with  a  demonstration  of 
the  need  for  adaptation  in  a  surface  ship  controller  used  in  restricted 
waters.  Sections  4  and  5  present  the  development  and  evaluation  of  two 
separate  approaches  for  the  parameter  estimation  function  of  the  outer  or  gain 
update  loop.  Section  4  presents  a  batch  type,  weighted-least-squares  para¬ 
meter  estimation  scheme  using  a  moving  window  of  ship  motion  data.  This  is  a 
statistical  approach  which  makes  no  use  of  the  probable  measurement  noise 
levels.  Section  5  presents  a  recursive  type,  minimum  variance  scheme.  This 
is  a  probabilistic  approach  which  utilizes  probabilistic  estimates  of  the 
measurement  noise  levels.  The  final  section,  Section  6,  presents  our  con¬ 
clusions  based  upon  this  work. 


2.  Problem  Formulation 


In  this  section,  we  formulate  the  surface  ship  path  control  problem  as 
a  linear,  state-variable  control  problem.  The  selection  of  measurements  and 
typical  process  disturbances  are  also  discussed. 


2 . 1  Equations  of  Motion. 


The  development  of  the  linearized,  state-variable  equations  of  motion 
for  a  surface  ship  moving  in  the  horizontal  plane  presented  here  is  based  on 
the  formulation  by  Fujino24  and  is  presented  in  more  detail  in  our  earlier 
work.1  The  coordinate  system  for  the  problem  is  shown  in  Fig.  2.  The  0-£r\ 
system  is  fixed  in  space  with  the  desired  ship  path  along  the  £-axis.  We 
assume  here  that  the  desired  path  is  a  straight  line.  This  simplification 
does  not  alter  the  essential  character  of  the  control  problem  and  is  typical 
in  maneuvering  situations  where  the  ship  is  to  follow  a  series  of  straight 
path  or  leading  line  segments.  The  CG-xy  system  is  fixed  at  the  center  of 
gravity  of  the  ship.  The  positive  sense  of  the  drift  angle  6,  heading  angle 
y,  yaw  rate  r,  and  rudder  angle  <5  are  shown.  Neglecting  the  effects  of  pitch 
and  roll,  the  ship  motion  can  be  described  by  coordinates  x,  y,  and 


Figure  2.  Coordinate  System  for  Path  Control 

The  exact  equations  of  motion  of  the  ship  are  integro-differential  equa¬ 
tions  in  which  convolution  integrals  represent  the  memory  effect  of  the  fluid 
to  previous  motion.25  An  alternative  formulation  yields  differential  equations 
with  frequency  dependent  coefficients.  Fujino25  has  shown  that  for  the  maneuver 


of  interest  here  the  frequency  dependence  is  negligible  and  constant-coefficient 
differential  equations  can  be  utilized.  This  assumption  becomes  less  and  less 
valid  as  H/T+l.  When  the  equations  of  motion  are  linearized  about  the  nominal 
path,  the  equation  in  the  x-coordinate  decouples  so  that  the  ship  motion  can 
be  given  by. 


dv 


(m+mv)  —  =  Y  v  +  (-mU+Yr ) r  +  Y*r  +  Y„6  +  Y  , 
r  dt  v  r  6 


dr 


(Izz+Jzz>  IE  =  Nvv  +  Nrr  +  Nvv  +  V  +  N  , 


(2) 

(3) 


(4) 


which  are  valid  for  small  deviations  from  the  equilibrium  straight-line, 
constant  speed  U  condition.  An  external  sway  force  Y  and  an  external  yawing 
moment  N  are  included  to  account  for  disturbances  which  act  on  the  ship,  it 
is  common  and  convenient  to  utilize  drift  angle  6  instead  of  the  lateral 
velocity  v  so  we  can  use, 

v  =  -UsinB  =  -U8  ,  (5) 


to  express  eq.  (2)  and  (3)  in  terms  of  the  drift  angle.  These  equations  can 
then  be  nondimensionalized  as  shown  in  the  Nomenclature  to  yield, 


dip ' 
dt' 

r*  , 

(6) 

-(m'+my')iV  = 

V8' 

+  (-m'+Y  , ) r ’  +  Yf.r'  +  Yu)'  +  y' 
r  r  6 

(7) 

•  dr ' 

V8' 

. 

(izz+j;z>  at*- = 

+  N_ i r '  +  N*  8'  +  N,,6'  +  N'  , 

(3  o 

(8) 

8’ 


(9) 


d<5 ' 
dt' 


Tr 


(5'  -  6’) 
c 


(10) 


where  we  have  now  included  a  first-order  model  for  the  steering  gear  dynamics. 
The  control  is  the  commanded  rudder  angle  6C’.  The  unit  of  nondimensional 
time  t'  is  the  time  it  takes  the  ship  to  travel  one  ship  length. 


Equations  (6)  through  (10)  can  be  transformed  into  state-variable  form; 

i.e.  , 
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or, 

n*l  m*l  qx  1 

x'  =  Fg  x'  +  Gg  tj  +  T s  w  .  (12) 

The  coefficients  of  the  open  loop  dynamics  matrix  f„  and  the  disturbance 
distribution  matrix  Y  —  are  algebraic  combinations  of  the  stability  deriva¬ 
tives  and  mass  and  inertia  terms  in  eq.  (7)  and  (8) . 1  For  this  study,  we 
utilize  data  obtained  by  Fujino2 4 ' 2 2 ' 2 8  for  a  model  of  the  290  m.  tanker 
Tokyo  Marti.  Fujino  conducted  planar  motion  mechanism  (PMM)  and  oblique  tow 
tests  of  the  model  at  various  water  depth-to-draf t  ratios  H/T.  Selected 
characteristics  for  this  vessel  are  shown  in  Table  1.  The  coefficients  f— 
and  Y^j  obtained  for  the  Tokyo  Maru  at  12  knots  full-scale  at  H/T  values  of 
1.30,  1.50,  1.89,  2.50,  and  00  are  given  in  Table  2.  As  shown  by  Fujino24 
this  vessel  is  course  unstable  for  the  intermediate  depth-to-draft  ratios  from 
about  3.0  down  to  1.75  as  is  typical  of  many  large  tankers. 


characteristic 

Fuj ino 1 s 
model 

prototype 

linear  scale  ratio,  X 

145.0 

- 

length  between  perpendiculars ,  m 

2.000 

290 

breadth,  m 

.  3276 

47.5 

draft,  m 

.1103 

16.0 

displacement 

58.4  kg 

179,100  LT 

block  coefficient 

0.8054 

0.8054 

rudder  area 

3,390.0  mm2 

71.29  m2 

propeller  diameter 

53.8  mm 

7.80  m 

P/D 

0.740 

0.740 

expanded  area  ratio 

0.619 

0.619 

number  of  blades 

5 

5 

Table  1.  Characteristics  of  Tokyo  Maru  Model  and  Prototype. 
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H/T 

- - - 

1.30 

1.50 

1.89 

2.50 

GO 

*22 

-1.6508 

-1.7136 

-1.7657 

-1.8177 

-1.9515 

f  23 

9.3157 

6.6235 

5.7359 

4.6112 

3.1591 

*25 

-0.55543 

-0.79235 

-0.88074 

-1.0416 

-1.0410 

Y21=*26 

346.69 

385.98 

477.68 

536.00 

567.13 

Y22=*27 

4.8040 

-2.2145 

-5.0043 

-5.8625 

2.3365 

f  32 

0.02974 

0.13890 

0.17199 

0.23621 

0.31507 

*  33 

-1.0388 

-0.71895 

-0.52766 

-0.54560 

-0.63651 

*35 

-0.09995 

-0.12092 

-0.15607 

-0.16639 

-0.16163 

Y  31_f 36 

11.825 

14.230 

21.141 

21.942 

16.844 

Y32=*37 

-19.216 

-23.123 

-28.233 

-31.490 

-37.384 

Table  2.  Coefficients  of  Tokyo  Maru  versus  H/T  at  Fn=0.116 

(12  knots  full-scale) 


2 . 2  Measurement  Selection. 

All  of  the  states  in  the  ship  path  control  problem  as  formulated  in  eq.  (11) 
are  available  for  measurement.  The  heading  ^ '  can  be  obtained  from  a 
compass;  the  yaw  rate  r'  can  be  obtained  from  a  rate  gyro;  the  drift  angle 
S'=-v'  can  be  obtained  from  a  doppler  sonar;  the  rudder  angle  5'  can  be  obtained 
from  the  rudder  stock  or  less  accurately  from  the  steering  gear  rams.  The 
lateral  offset  from  the  desired  path  g*  must  be  obtained  using  navigation  aids 
such  as  DECCA  Hi-Fix  or  radar.  Each  of  these  measurements  may  be  subject  to 
bias  and  zero-mean  measurement  errors  and  system  transmission  noise.  In  the 
presence  of  this  measurement  "noise"  and  with  the  measurement  of  only  selected 
states,  the  complete  state  vector  can  be  estimated  using  a  Kalman  filter1'29,30 
provided  all  of  the  states  are  observable  with  the  chosen  measurements. 

The  authors  have  previously  studied  the  observability  of  the  ship  path 
control  problem. 1  It  was  shown  that  it  is  necessary  to  measure  the  lateral 
offset  n’ .  Additional  measurements  improve  the  ability  of  a  Kalman  filter  to 
estimate  all  the  states  and  thus  improve  the  effectiveness  of  an  optimal  state 
feedback  controller.  The  yaw  rate  r’  is  the  next  most  effective  measurement. 

The  heading  \ji'  is  readily  available  and  is  the  next  most  effective  measure¬ 
ment.  The  drift  angle  8'  measurement  was  shown  to  add  little  to  the  effec¬ 
tiveness  of  a  ship  path  controller  which  already  measures  n',r’,  and  i|i  1  .  With 
the  steering  gear  model  used  in  Section  2.1,  there  is  little  need  to  measure 
the  rudder  angle  since  the  state  is  known  exactly  given  any  initial  condition 
<5'(t0)  and  the  subsequent  rudder  command  history  <V(t),  t^t0.  In  any  practi¬ 
cal  application,  the  steering  gear  would  have  its  own,  separate  feedback 
system;  the  first-order  model  included  in  eq.  (11)  is  just  a  means  of  intro- 
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ducing  a  realistic  rudder  time  response  into  our  study.  For  the  purposes  of 
the  inner  or  control  loop  design,  it  is  reasonable  to  assume  a  measurement 
vector  consisting  of  measurements  of  t|i',  r',  and  n  ’  each  contaminated  by 
Gaussian,  white  noise;  i.e., 


p*l 

ro 

o 

o 

o 

"vi 

z  = 

0  10  0  0 
_0  0  0  1  o_ 

X  + 

v2 

Lv3J 

The  concept  of  observability  relates  to  the  feasibility  of  estimating 
the  states  of  a  system  from  a  particular  set  of  measurements.  In  an  adaptive 
system  (or  in  off-line  parameter  estimation) ,  a  second  objective  is  the 
estimation  of  the  parameters  of  the  system  from  a  particular  set  of  measure¬ 
ments  and  the  system  input  history.  Identif iability  is  the  dual  of  observ¬ 
ability  in  this  context.  Astrom  and  Kallstrom18  have  studied  the  parameter 
identif iability  of  the  open-loop  dynamics  matrix  parameters  f ^  in  the  r' 
and  the  6'  equations  in  eq.  (11).  They  showed  that  these  parameters  are 
identif iable  provided  both  4* •  and  v‘  are  measured.  This  conclusion  requires 
that  the  ship  is  controllable  with  the  rudder  which  is  known  to  be  generally 
true  and  which  has  been  demonstrated  theoretically  by  the  authors.1  In  the 
path  control  problem,  ip'  and  n'  can  provide  the  same  identif  iability  as  ip' 
and  v'  as  can  be  seen  from  eq.  (9).  Regardless,  it  is  advisable  to  use  all  fea¬ 
sible  measurements  in  the  parameter  estimation  so  we  assume  that  all  states 
are  measured  for  use  in  the  outer  or  gain  update  loop  in  Fig.  1.  This  gives, 


n*l 

z '  =  x  +  v' 


(14) 


This  full  measurement  vector  is  also  used  in  some  of  the  inner  or  control 
loop  designs  which  are  presented  in  the  following  sections. 


The  final  part  of  measurement  definition  is  to  establish  reasonable 

levels  for  the  measurement  noise  v  and  v'.  The  white  noise  power  spectral 

density  needed  in  our  continuous  system  design  approach  can  be  estimated  by 

assuming  the  noise  to  be  exponentially  correlated  with  an  RMS  noise  level  a 

and  a  correlation  time  t . .  The  t .  should  be  much  faster  than  the  time  con- 

1  3 

stants  of  the  ship  and  less  than  the  system  sampling  time  for  the  white 
noise  model  to  be  valid.  The  power  spectral  density  can  then  be  estimated 
by, 
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(15) 


To  evaluate  the  control  system  effectiveness  in  this  study,  we  use  digital 
simulation  with  a  f ixed-stepsize  Euler  integration  scheme.  This  has  the 
effect  of  approximating  the  continuous  Gauss-Markov  process,  eq.  (11)  and 
eq.  (13),  by  a  discrete  Gauss-Markov  process.  In  these  simulations,  the 
covariance  of  the  computer  generated  random  measurement  noise  must  be  selected 
to  be  consistent  with  the  design  noise  power  spectral  density.  To  provide 
equivalent  state  estimate  error  covariances,  it  is  necessary  that  the  simula¬ 
tion  measurement  noise  variance  given  by. 


where  At  is  the  integration  stepsize . 1 ' 3 1  This  can  also  be  considered  from 
a  more  direct  viewpoint.  If  the  controller  is  implemented  digitally  in  an 

onboard  computer  with  the  system  sampled  each  At,  the  measurement  noise  will 

,  ? 

be  a  white  sequence  with  variance  Oj 


The  reference  measurement  noise  levels  used  in  this  work  are  shown  in 

Table  3.  In  view  of  our  earlier  comments  about  the  rudder  model  included  in 

eq.  (11),  we  assume  exact  knowledge  of  the  rudder  angle.  Astrom  and  Kallstrom18 

note  that  all  sensors  have  dynamics  with  time  constants  less  than  1  sec.  and 

that  the  measurement  errors  are  about  0.1°  in  ip ,  0.02°/s  in  r  and  0.01  m/s 

in  v.  Millers32  uses  RMS  errors  of  0.2°  in  \p ,  0.01°/s  in  r  and  10  m.  in  g. 

Canner33  states  that  DECCA  Hi-Fix  crosstrack  errors  are  as  low  as  lm.  when 

the  baseline  is  along  the  desired  path  as  is  done  at  the  entrance  to  Europoort. 

Astrom  and  Kallstrom18  and  Bystrom  and  Kallstrom1 ^  have  found  errors  in  r  of 

less  than  0.002  °/s  in  systems  identification  of  full-scale  experiments.  In 

view  of  this  data,  the  reference  levels  assumed  in  Table  3  seem  reasonable. 

The  values  for  r^^  and  a.'  are  nondimensional .  The  o/  are  calculated  by 
33  j  3 

eq.  (16)  assuming  At '=0.005  (At=.24  s)  which  we  use  in  our  simulations. 

A  noise  correlation  time  of  0.1  s.  and  a  sampling  time  of  .24  s.  imply  some 


correlation  but  the  resulting  values  for  o  j '  yield  reasonable  covariances 


for  a  white  measurement  noise  sequence, 
percent  of  the  assumed  values  for  Oj. 


Dimensional  1 y , 


the  it.  ' 
3 


are  about  93 
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measurement 

source 

0  . 

3 

T  . 

-  —  3 _ 

r 

33 

0  .  ' 

3 

compass 

0.1° 

0 .  Is 

1 ,298xl0~8 

1.611*10-3 

r’ 

rate  gyro 

0 . 01 °/s 

0.1s 

2.H60*10~7 

7 . 563*10-3 

8’ (v‘) 

doppler  sonar 

0.01  m/s 

0.1s 

1. 118*1 0-8 

1.459*10" 3 

n  1 

DECCA  Hi -Fix 

3m 

0.1s 

4  .  559*10~7 

9. 549*10_3 

Table  3.  Reference  Measurement  Noise  Characteristics 


2 . 3  Process  Disturbances. 

C/ 

While  operating  in  restricted  waters,  a  ship  can  be  subjected  to  a 
wide  range  of  disturbances.  Many  of  these  can  be  characterized  as  being 
short-term  relative  to  the  time  constants  of  the  ship  and  as  having  essentially 
a  zero  mean  value.  First-order  wave  forces,  wind  gusts,  bottom  ripples  as 
observed  outside  Europoort , 34  and  passing  ships  can  be  included  in  this 
category.  Other  disturbances  remain  long  enough  relative  to  the  time  constants 
of  the  ship  that  they  must  be  considered  to  have  non-zero  mean  value.  Second- 
order  wave  forces  and  the  effect  of  a  lateral  current,  bank,  or  steady  wind 
are  included  in  this  category.  For  the  purposes  of  this  study,  we  define  a 
number  of  typical  or  design  process  disturbances  which  are  used  in  simulations 
to  evaluate  the  performance  of  the  various  path  controllers.  These  design 
disturbances  are  defined  briefly  here. 

Passing  Ship.  The  lateral  force  Y'  and  yawing  moment  N'  due  to  a 
passing  ship  was  selected  as  a  typical  short-term, essential ly  zero-mean  dis¬ 
turbance.  The  assumed  design  disturbance  is  shown  in  Fig.  3.  This  disturbance 
is  based  on  results  originally  presented  by  Newton'  for  two  Mariner  vessels 
passing  in  deep  water.  These  results  are  considered  to  be  representative 
forces  and  moment  histories  and  thus  are  reasonable  for  use  in  comparisons 
here.  Yung36  and  Abkowitz,  Ashe,  and  Fortson37  show  that  the  magnitude  of 
the  disturbances  increase  in  shallow  water  as  H/T-*]  so  the  magnitudes  in  Fig.  3 
are  known  to  be  low  at  the  shallower  depths.  In  Fig.  3,  the  nondimensional 
time  scale  is  in  ship  lengths  and  the  ships  are  beam-to-beam  at  t,=0.  These 
lateral  force  and  yawing  moment  histories  are  assumed  to  be  independent  of 
depth  under  keel  in  our  simulations. 


Figure  3.  Design  Passing  Ship  Disturbance. 


Lateral  Current.  The  effect  of  a  lateral  current  was  selected  as  a 
typical  long-term,  nonzero-mean  disturbance  for  use  in  our  ship  path  controller 
simulations.  In  a  steady  current,  the  ship  assumes  an  equilibrium  condition 
with  6'=  0  and  i//'=8'  so  that  the  effective  drift  angle  relative  to  the  water 
Se'  is  zero.  In  this  condition,  there  is  no  external  hydrodynamic  lateral 
force  or  yawing  moment  on  the  ship.  In  our  ship  maneuvering  equations,  we 
have  assumed  to  this  point  that  the  drift  angle  8'  is  with  respect  to  the 
earth.  In  shallow  water,  a  doppler  sonar  would  actually  measure  lateral 
velocity  relative  to  the  bottom.  In  a  lateral  current  vc'  without  an  additional 
disturbance,  eq.  (7),  (8),  and  (9)  should  properly  be  written. 


Now  if  the  drift  angle  relative  to  the  earth  H‘  is  introduced, 


=  8 '  t  v( 


(JO) 


eq.  (17),  (18),  and  (19)  become  as  follows  in  a  steady  current : 


<m'+my,,dt'=  Y*lP'  +  (-m'+Y^Jr'  t  Y;,r*  *  Y^.S'  *  Y,,vt 


(21) 


Vs'  *  Vr'  *  *  V6'  *  "c'c'  ■ 


(72) 


dn’ 
dt 1 


=  *' 


-  IV 


(2  1) 


Thus  when  using  the  drift  angle  with  respect  to  the  earth  in  eq.  (11),  a 
steady  current  has  the  effect  of  applying  an  external  lateral  force  and 
yawing  moment  given  by, 

Y’  =  Y,’v  ’  (24) 

8  c 

and, 

N 1  -  Np'vc'  .  (25) 


For  design  evaluation  pur|»oses,  we  have  used  eq.  (24)  and  C 15)  to 
establish  the  lateral  force  and  yawing  moment  produced  by  a  1  knot  lateral 
current  on  the  Tokyo  Maru  moving  at  17  knots,  in  an  intermediate  water  depth 
H/T=1.89.  This  disturbance  was  assumed  to  be  constant  for  15  ship  lengths 
and  then  to  reduce  linea'ly  to  one  half  this  value  at  20  ship  lengths.  This 
design  disturbance  is  shown  in  Fig.  4.  These  lateral  force  and  yawing  moment 
histories  are  assumed  to  be  independent  of  depth  under  keel  in  our  simulation;,. 


i 

) 

i 
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Figure  4.  Design  Lateral  Current  Disturbance 
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ABC  Harbor  Entrance.  As  a  final  design  evaluation  disturbance  we 
selected  an  approximation  to  the  lateral  current  disturbance  experienced  at 
the  entrance  to  Europoort  (Rotterdam).  This  lateral  current  was  included 
in  the  ABC  harbor  definition  recommended  by  the  SNAMF  H-10  Controllability 
panel  for  use  in  ship  controllability  studies. 3ft  This  proves  to  be  a  very 
demanding  test  for  a  path  controller.  There  is  about  a  1  knot  lateral  current 
outside  the  entrance  which  increases  up  to  about  2  knots  at  about  0.4  km 
from  the  harbor  entrance.  The  lateral  current  then  reduces  to  zero  at  the 
harbor  entrance.  In  this  situation,  the  lateral  current  produces  two  types 
of  forces  and  moments  on  the  ship  in  addition  to  those  resulting  from  the 
use  of  the  drift  angle  with  respect  to  the  earth  in  eq.  (11).  First,  the 
changing  current  will  cause  the  ship  to  accelerate  with  respect  to  the  earth 
so  additional  inertial  forces  and  moments  will  appear  in  eq.  (7)  and  (8). 
Second, the  spatial  current  change  over  the  length  of  the  ship  will  produce  a 
changing  lateral  force  along  the  hull  which  can  be  integrated  stripwise 
similar  to  the  approach  used  by  Newman3 ^  to  yield  a  total  lateral  force  and 
yawing  moment  on  the  ship. 

All  of  the  effects  of  the  steady  and  changing  lateral  current  at  the 
entrance  to  the  ABC  harbor  can  be  moved  to  the  right  side  of  eq.  (7)  and  (8) 
and  combined  as  net  lateral  force  and  yawing  moment  terms.  With  &'  assumed 
to  be  with  respect  to  the  earth,  this  net  result  is  a  disturbance  history 
similar  to  that  shown  in  Fig.  5.  This  disturbance  history  is  based  on 
similar  results  presented  by  van  Oortmerssen4 0  and  has  been  verified  by  the 
authors  using  a  simplified  strip  theory  approach.  In  Fig.  5,  the  ship's  bow 
reaches  the  harbor  entrance  (vc'=0)  at  t'=11.0;  the  stern  thus  reaches  the 
entrance  at  t'=12.0.  The  steady  values  for  t'<6  in  Fig.  5  have  approximately 
the  same  magnitude  (1  knot  lateral  current)  as  those  in  Fig.  4  for  t<15. 
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Figure  5.  Lateral  Current  Disturbance  at  ABC  Harbor  Entrance 


3.  Control  Loop  Design 


In  this  section,  we  discuss  the  design  of  the  inner  or  control  loop  m 
Fig.  1.  This  naturally  includes  the  methods  used  for  the  calculation  of 
controller  gains  C  and  the  Kalman  filter  gains  K  performed  in  the  outer  or 
gain  update  loop.  This  section  involves  an  extension  of  our  earlier  work1 
and  some  of  the  introductory  discussion  included  there  and  elsewhere* q ' ** 1 
will  not  be  repeated  here.  We  begin  with  a  discussion  of  the  optimal,  steady- 
state  controller  design  and  then  proceed  with  development  of  the  Kalman  filter 
which  generates  the  state  estimate  used  by  the  controller.  The  section 
closes  with  a  look  at  the  need  and  requirements  for  an  adaptive  feature  in  a 
path  controller  used  in  restricted  waters. 

3. 1  Optimal,  Steady-State  Controller  Design. 

The  control  problem  is  to  keep  the  ship  along  the  desired  track  usinq 
acceptable  levels  of  rudder  in  the  presence  of  disturbances  while  usinq  noisy 
measurements.  The  ship  model  is  given  by  eq.  (12)  where  in  this  subsection 
we  will  assume  the  disturbance  vector  w  consists  of  a  Gaussian  white  noise 
lateral  force  and  yawing  moment.  Equation  (12)  describes  a  Gauss-Markov 
process  in  which  the  state  of  the  ship  is  completely  described  by  its  mean 
value  state  vector  x  which  follows  eq.  (12)  with  w=0  and  by  its  covariance 
matrix  X  where, 

x  (t)  =  E  (x  (t)  ’  ,  (26) 

X (t )  =  E ( (x (t)  -  x  (t) )  (x (t)  -  x(t))T)  ,  (27) 

and  E [ . . . ]  is  the  expected  value  or  ensemble  average  over  the  many  possible 
observations  at  time  t.  We  will  assume  that  the  white  noise  disturbance  has 
zero  mean  and  covariance  given  by, 

E[w(t)w(T)TJ  =  Q (t ) 6 (t— t )  ,  (28) 

where  Q  is  the  power  spectral  density  matrix  and  6(t-i)  is  the  Dirac  delta 
function.  We  also  assume  there  is  no  correlation  between  the  initial  condition 
of  the  system  x(t0)  and  the  disturbances  w. 

An  optimal  control  design  problem  can  be  stated  as  establishing  the 
linear  state  feedback  control  law, 

u  =  Cx  ,  (29) 
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which  minimizes  the  linear  quadratic  eo-.t  t  un.  1 1<  iti , 


(xTAx  ♦  uTBu)dt 


where  A  and  B  are  weighting  matrices  at  the  designer's  disjxisal.  In  tins 
problem,  we  wish  to  minimize  n‘  while  usmq  acceptable  levels  of  rudder  a’ 
and  6C'  .  We  can  therefore  assume  all  terms  in  the  A  matrix  to  be  zero 
except  the  a  weighting  on  n'?  and  the  a  weightinq  on  * ' '  and  assume-  a 


nonzero  B=b^  weighting  of  5C'^.  To  yield  a  reasonable  weighting  among  these 
terms  and  to  accomplish  the  needed  scaling  of  terms  in  eq.  (30)  we  assume. 


a44  -  <^)"2  =  772 . 4b  3 


a55  =  <*o> 


=  131.332 


U  =  <‘V)'2  =  131.3  32  , 


based  on  a  dimensional  nQ  10.43  m.  (slightly  less  than  one-quarter  beam)  and 

a  dimensional  'S  =6r  =5°.  This  indicates  that  we  are  willing  to  commit 
o  <-o 

approximately  5°  rudder  to  path  cortrol  when  the  ship  deviates  n  from  the 
desired  path. 

The  Separation  Theorem2"3  states  that  the  optimal  control  desired  here 
is  obtained  by  assuming  exact  knowledge  of  the  state  for  use  in  cq.  (29) 
and  by  neglectinq  the  disturbance  w  in  eq.  (12).  When  the  matrices  Fs ,  Gs, 
A,  and  B  are  constant,  the  control  which  minimizes  eq.  (30)  is  given  by, 

C  =  -B-1GsTS  ,  (34) 

CO 

where  is  the  steady-state  (S=0)  solution  of  the  matrix  Riccati  equation, 

S  =  -SFs  -  FsTS  +  SGsB-1C,sTS  -  A  .  (  3 r> ) 


An  efficient  way  to  obtain  is  by  Potter's  algorithm  which  utilizes  tin- 
technique  of  eigenvector  decomposition  to  obtain  a  closed-form  solution  for 
S^.  MacFarlane42  and  Potter43  first  proposed  this  solution  technique.  It 
was  developed  into  a  practical  design  tool  by  Bryson  and  Hall4'*  who  utilized 
the  QR  algorithm  to  obtain  the  eigensystem  in  their  OPTSYS  computer  program. 
Potter43  proved  that  can  be  expressed  in  terms  of  the  eigenvectors  of  the 
2n*2n  partitioned  matrix, 
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which  can  be  shown  to  be  the  coefficient  matrix  of  the  2n  Eulei -Lagrange 
equations  for  minimizing  eq.  (30)  subject  to  eq.  (12) .  The  eigenvalues  of 
M  are  in  pairs  which  are  symmetric  about  the  imaginary  axis  in  the  complex 
plane.  If  the  eigenvectors  associated  with  the  eigenvalues  with  positive 
real  parts  are  designated  with  the  subscript  +  and  those  associated  with  the 
eigenvalues  with  negative  real  parts  are  designated  with  the  subscript  -  , 
the  eigenvectors  can  be  arranged  as  follows: 


(37) 


The  steady-state  solution  of  eq.  (35)  is  then  given  by, 

S  =  A_(X_)_1  .  (38) 

oo 

Equations  (34)  and  (38)  are  used  to  calculate  the  controller  gain  matrices 
in  this  work.  The  QR  algorithm  is  used  to  obtain  the  eigensystem  of  M  as 
defined  by  eq.  (36) . 


3 . 2  First-Order  Shaping  Filter  Approach. 

In  Section  3.1,  we  assumed  that  the  disturbance  vector  w  in  eq.  (12) 
was  composed  of  Gaussian  white  noise.  In  our  earlier  work1,  we  developed 
path  controllers  capable  of  handling  short-term,  essentially  zero-mean 
disturbances  such  as  the  passing  ship  disturbance  shown  in  Fig.  3.  This 
disturbance  acts  in  a  correlated  manner  over  about  three  ship  lengths  so  the 
assumption  of  Gaussian  white  noise  disturbances  is  poor  for  this  problem.  The 
usual  approach  when  a  disturbance  is  correlated  with  characteristic  times 
comparable  to  those  of  the  system  under  consideration  is  to  model  the  distur¬ 
bance  as  the  output  of  a  shaping  filter  which  is  driven  by  Gaussian  white 
noise.  The  passing  ship  disturbance  in  Fig.  3  can  be  reasonably  modeled  as 
exponentially  correlated  processes  with  zero  means.  These  can  be  produced  by 
first-order  shaping  filters.  We  therefore  introduced  the  external  yawing 
moment  N'  and  the  external  lateral  force  Y'  as  additional  states  defined  by: 


_  dN’ 

tn  ^  =  -N  +  wi 


(39) 
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and 
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where  correlation  times  T^=T  -1  ship  lenqth  are  reasonable  assumptions  based 

on  Fig.  3.  These  shaping  filter  equations  are  driven  by  Gaussian  white 
T 

noise  w=(w^,w2)  which  has  power  spectral  density  Q.  This  model  is  equally 
valid  for  the  disturbance  produced  as  the  ship  passes  localized  fixed  objects 
and  other  short-term  disturbances. 

To  incorporate  the  disturbance  model  eq.  (39)  and  (40)  into  the  controller 
design,  the  state  vector  x  is  augmented  by  the  addition  of  states  N'  and  Y' 
to  yield  an  augmented  state  vector  x' .  Equation  (11)  then  becomes. 
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x'  =  Fe  x'  +  Ge  u  +  re  w  .  (42) 

The  upper  right  (5*2)  partition  of  I’e  in  eq.  (41)  is  I's  in  eq.  (11);  i.e., 

f  => etc.  The  measurement  eq.  (13)  now  becomes, 

2  d  21 

z  =  He  x'  +  v  ,  (43) 

where  He  is  just  Hg  in  eq.  (13)  with  a  (3*2)  block  of  zeros  added  at  the 

right.  The  optimal  control  law  is  now  qiven  by, 

u  =  Cx'  ,  (44) 

where  the  control  gain  matrix  is  obtained  as  described  above  except  the  A 
matrix  must  be  extended  to  (7*7)  with  additional  rows  and  columns  of  zeros 

for  use  with  jc'  in  lieu  of  x  in  eq.  30  and  Ke  and  G0  must  be  used  in  eq.  (34), 

(35),  and  (36)  in  lieu  of  Fs  and  Gg.  The  optimal  control  qains  are  then 
given  by, 


'*&***>■>■■ 


(45) 
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Having  only  the  noisy  measurements  of  three  states  included  in  z,  the 
state  vector  x/  needed  in  eq.  (44)  is  not  directly  available.  The  Separation 
Theorem  states  that  it  is  optimal  to  use  the  control  gains  produced  by  eq.  (45) 
with  a  maximum  likelihood  estimate  of  the  augmented  state  vector  x'  produced 
by  an  optimal  stochastic  state  estimator  or  Kalman-Bucy  filter.  The  control 
law  then  becomes, 

u  =  Cx*  ,  (46) 

where  the  state  estimate  is  produced  by  the  Kalman-Bucy  filter  given  by, 

x'  =  Fe  x'  +  Ge  u  +  K(z-He  x')  .  (47) 

When  the  system  is  stationary;  i.e.,  a  statistical  steady-state  with  Fe ,  Ge, 

Te,  Q,  and  R  constant,  the  steady-state  filter  gain  matrix  K  is  given  by, 

K  =  P  HeTR-1  ,  (48) 

OO 

where  P  is  the  steady-state  (P=0)  solution  of  the  matrix  Riccati  equation, 

OO 

P  =  FeP  +  PFeT  +  reQreT  -  PHeTR-,HeP  ,  (49) 


and  Q  is  the  power  spectral  density  of  the  white  noise  process  disturbance  w 
in  eq.  (42)  and  R  is  the  power  spectral  density  of  the  white  measurement 
noise  v  in  eq.  (43) . 


An  efficient  way  to  obtain  P  is  again  by  Potter's  algorithm  which 

00 

utilizes  eigenvector  decomposition.  In  this  case,  P^  can  be  expressed  in 
terms  of  the  eigenvectors  of  the  2n'*2n'  partitioned  matrix, 


T  -i  1 
CHe  R  lHe  ' 


(50) 


If  the  eigenvectors  associated  with  the  eigenvalues  of  M'  with  positive 
parts  are  designated  with  the  subscript  +  and  those  associated  with  the 
eigenvalues  with  negative  real  parts  are  designated  with  the  subscript  - 
the  eigenvectors  can  be  arranged  as  follows: 

fx;  j  xi] 


real 


(51) 
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The  steady-state  solution  of  eq.  (49)  is  then  given  by, 


P»  =  -x;(A;)-1  .  (52) 

Equations  (48)  and  (52)  are  used  to  calculate  the  Kalman  filter  gain  matrices 
in  this  work.  The  QR  algorithm  is  used  to  obtain  the  eigensystem  of  M1  as 
defined  by  eq.  (50). 

At  this  point,  we  will  present  typical  results  for  a  path  controller 
designed  using  the  first-order  shaping  filter  disturbance  models  eq.  (39) 
and  (40) .  These  results  illustrate  the  effectiveness  and  the  limitations 
of  this  type  of  design.  The  authors  have  previously  shown1  that  if  a  non- 
adaptive  path  controller  is  designed  assuming  constant  ship  characteristics, 
the  best  overall  performance  is  achieved  with  a  controller  designed  for  the 
least  course  stable  ship  characteristics.  For  the  Tokyo  Mar u  this  condition 
is  at  about  H/T=1.89  so  the  characteristics  for  this  depth-to-draft  ratio  in 
Table  2  are  used  in  this  example.  The  measurement  noise  power  spectral 
density  R  is  taken  from  Table  3  assuming  measurements  of  ij> 1  ,  r',  and  n '  in 
eq.  (43).  The  nonzero  terms  in  the  controller  cost  function  weighting 
matrices  are  defined  by  eq.  (31),  (32),  and  (33).  The  remaining  definition 
needed  is  the  process  disturbance  power  spectral  density  Q.  Using  the 
passing  ship  disturbance  in  Fig.  3  as  the  design  disturbance,  the  RMS 
values  of  N'  and  Y'  between  t'=-2  and  t'=1.4  are  oN=8. 798*10"^  and 
ay=21 . 178* lO-8 ,  respectively.  Assuming  correlation  times  TN=TY=1  as  noted 
above,  the  nonzero  terms  in  the  power  spectral  density  matrix  can  be  taken 
as , 

qn  =  2TNaN2  =  1.548*10~8  ,  (53) 

and , 

q22  =  2Tyay2  =  8.970xl0-8  .  (54) 

The  control  loop  was  designed  for  the  Tokyo  Maru  with  the  above 
assumptions  using  the  Michigan  Terminal  System  (MTS)  version  of  Bryson  and 
Hall's  OPTSYS  program . 1  ' 44  The  resulting  controller  gains  and  Kalman- 
Bucy  filter  gains  are  shown  in  Table  4.  These  results  are  different  from 
those  reported  earlier  by  the  authors1  due  to  the  use  of  a  corrected  Froude 
number  of  0.116,  the  use  of  a  zero  weighting  a^  on  \p '  *  in  eq .  (30),  and 
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the  use  of  lower  assumed  measurement  noise  levels.  The  zero  row  in  the  K 
matrix  results  because  with  our  steering  gear  model  the  control  history 
produces  exact  knowledge  of  the  rudder  angle.  Therefore,  the  measurements 
provide  no  additional  information.  All  the  nonzero  gains  in  Table  4  must 
vary  with  water  depth  if  they  are  to  remain  optimal  except  c  = 2.4252  which 
the  authors  have  shown  to  be  vessel,  speed,  and  water  depth  independent. 


T 

controller  gams  C 

Kalman 

-Bucy  filter  gains  K 

5.5421 

4.6335 

0.8400 

-0.0021 

2.6601 

18.5075 

12.2263 

-0.3278 

6.3895 

3.3586 

2.6208 

-4.2929 

2.4252 

-0.0753 

-0.5225 

2.9000 

-0.8499 

0.0000 

0.0000 

O.uOOO 

679.68 

0.1164 

0.1850 

0.0440 

-52.776 

-0.0982 

-0.1604 

0.1654 

Table  4.  Optimal  Gains  for  Tokyo  Maru  at  H/T=1.89  and  Fn=0. 116- 
First-Order  Shaping  Filter  Disturbance  Model 

To  illustrate  the  effectiveness  of  the  design  given  in  Table  4,  we 
simulated  the  Tokyo  Maru  under  the  control  of  this  controller  using  the 
SHIPSIM/OPTSIM  computer  program.31  Figures  6  and  7  show  the  results  of  a 
simulation  which  begins  with  the  ship  one  half  beam  (B/2=23.75  m.)  off  the 
desired  track  with  all  other  states  zero.  The  filter  is  assumed  to  have  been 
operating  for  some  time  and  is  given  exact  knowledge  of  the  state  x=x  at 
t'=0.  In  this  simulation,  the  ship  also  passes  another  ship  beam-to-beam  at 
t'=20.  The  magnitude  of  this  passing  ship  disturbance  is  arbitrarily  taken 
as  four  times  that  shown  in  Fig.  3.  Figure  6  shows  the  resulting  lateral 
offset  response.  The  maximum  overshoot  at  t'=5.3  is  1,6  m. ;  the  maximum 
deviation  due  to  the  passing  ship  is  6.2  m.  Figure  7  shows  the  corresponding 
rudder  angle  response.  The  controller  provides  effective  control  under 
these  conditions.  It  provides  effective  control  when  the  ship  is  subjected 
to  the  short-term,  essentially  zero  mean  disturbances  for  which  the  controller 
was  designed. 

The  weakness  of  the  ship  path  controller  designed  using  first-order 
shaping  filters,  as  defined  above,  to  model  the  process  disturbances  is  that 
the  controller  cannot  accommodate  more  long-term  disturbances  without  a  non¬ 
zero  offset  from  the  desired  path.  Normal  state-variable  feedback  controllers,, 
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Figure  6.  Lateral  Offset  Response  to  B/2  Initial  Offset  and 
Passing  Ship  -  First-Order  Shaping  Filter  Disturbance  Model  Design 


Figure  7.  Rudder  Angle  Response  to  B/2  Initial  Offset  and  Passing 
Ship  -  First-Order  Shaping  Filter  Disturbance  Model  Design 


as  this  approach  produces,  yield  a  proport ional  plus  derivative  (P-D)  type  control 
law.  To  produce  zero  steady-state  errors  with  a  constant  disturbance  some  type 
of  integral  action  is  needed;  i.e.  .1  (P-I-D)  type  control  law  is  needed.  To 
show  the  performance  of  the  controller  given  in  Table  4  with  constant  distur¬ 
bances,  the  Tokyo  Maru  was  simulated  under  control  of  this  controller  while 
being  subjected  to  the  design  lateral  current  disturbance  shown  in  Fig.  4.  The 
initial  condition  is  x=x=0  in  this  simulation.  The  lateral  offset  response  is 
shown  in  Fig.  8.  Steady  offsets  of  about  54  m.  and  27  m.  occur  during  the  two 
periods  of  constant  disturbance.  A  1  knot  lateral  current  produces  a  steady 
offset  of  over  one  ship  beam.  The  corresponding  rudder  angle  response  is  snown 
in  Fig.  9.  The  yawing  moment  estimate  N’  produced  by  the  Kalman-Bucy  filter  is 
shown  in  Fig.  10.  The  filter  produces  approximately  constant  but  biased  estimates 
of  N'  and  Y*  during  the  periods  of  constant  disturbance.  The  state  estimates 
have  much  smaller  or  no  apparent  bias  errors.  Comparing  Fig.  10  with  Fig.  4,  the 
filter  estimates  the  0.00103  yawing  moment  during  the  first  15  ship  lengths  to 
be  about  0.00072.  The  equilibrium  rudder  angles  in  Fig.  9  are  zero  and  the 
equilibrium  heading  and  drift  angles  are  equal.  Thus,  from  sq.  (46)  it  can  be 
seen  that  a  steady  offset  n '  develops  to  "counteract"  the  steady  errors  in  the 
disturbance  estimates  N'  and  Y'. 
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Figure  8.  Lateral  Offset  Response  to  Desiqn  Lateral  Current  - 
First-Order  Shaping  Filter  Disturbance  Model  Design 
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Figure  9.  Rudder  Angle  Response1  to  Design  Lateral  Current 
First-Order  Shaping  Filter  Disturbance  Model  Design 


Figure  10.  Yawing  Moment  Estimate  with  Design  Lateral 
Current  -  First-Order  Shaping  Filter  Disturbance  Model  Design 
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3.3  Random  Walk  Approach. 


t 
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As  noted  in  the  previous  section,  some  type  of  integral  control  action  is 
needed  to  produce  a  zero  offset  of  the  ship  from  the  desired  path  with  a  con¬ 
stant  disturbance.  Attention  has  been  directed  recently  toward  the  develop¬ 
ment  of  stochastic  state  variable  controllers  with  integral  control  properties 
when  using  only  incomplete,  noisy  state  measurements.  Among  this  work  is  that 
of  Holley  and  Bryson45  in  which  they  applied  a  number  of  approaches  to  the 
control  of  an  aircraft  landing  in  a  crosswind.  Alternative  approaches  which 
are  similar  to  that  which  we  present  here  have  been  presented  by  Kwatny46 
and  by  Balchen,  et  al.47 

A  reasonable  first  attempt  at  the  treatment  of  a  constant  disturbance 
might  be  to  model  the  disturbance  as  a  random  bias;  i.e., 

=  0  ,  N'  (0)  =  NQ 1  ,  (55) 

where  N0'  is  an  unknown,  random  quantity.  In  this  problem,  models  of  this 

type  could  be  assumed  for  N'  and  Y*  and  the  state  vector  could  be  augmented  to 
T  T 

become  x'=(x  ,N',Y']  as  done  above  with  the  first-order  shaping  filters. 

Holley  and  Bryson45  have  shown,  however,  that  the  steady-state  Kalman  filter, 
as  used  here,  does  not  exist  for  this  problem.  This  results  because  a  time- 
varying  Kalman  filter  would  yield  statistically  exact  estimates  of  the  random 
biases  so  that  in  the  steady-state  the  measurements  will  contain  no  additional 
information  on  these  disturbances. 

A  number  of  successful  approaches  do,  however,  exist  for  this  control 
problem.  First,  it  is  feasible  to  model  the  disturbances  as  first-order  shaping 
filters  as  in  the  previous  section  but  with  the  assumption  that  the  correlation 
times  are  very  long  compared  with  the  dynamics  of  the  ship.  Disturbance  models 
eq.  (39)  and  (40)  would  be  used  with  the  assumption  that  TN  and  Ty  are  perhaps 
100  ship  lengths.  It  would  be  necessary  to  assume  RMS  values  of  the  distur¬ 
bances  oN  and  Oy  for  use  in  eq.  (53)  and  (54)  and  the  design  could  then  proceed 
as  in  the  previous  section.  This  approach  would  continue  to  produce  biased 
estimates  of  constant  disturbances  but  these  estimates  approach  the  correct 
values  as  the  correlation  times  are  increased.45  Note  that  in  this  approach 
it  is  necessary  to  make  assumptions  of  both  the  correlation  times  and  the 
RMS  disturbance  values. 
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A  second  approach  was  presented  recently  by  Holley  and  Bryson.145  They 
develop  a  state  variable  generalization  of  integral  control  in  which  a  Kalman 
filter  is  used  to  estimate  the  state  vector  x  ignoring  the  disturbances  and 
then  an  integral  of  the  output  deviations  is  added  separately  to  the  control 
law.  They  note  that  this  approach  produces  a  degraded  transient  response 
compared  with  the  use  of  first-order  shaping  filters  with  long  correlation 
times  but  that  the  controller  is  less  sensitive  to  modeling  errors.  The 
authors  will  be  undertaking  an  investigation  of  this  type  of  controller  for 
ship  path  control  during  the  coming  year  and  these  results  will  be  reported 
separately  later. 

A  third  approach  to  the  development  of  a  state  variable  controller  with 
integral  control  properties  is  to  model  the  constant  or  nearly  constant  dis¬ 
turbances  as  a  random  walk.  More  precisely  this  model  is  an  independent- 
increment,  continuous  process  and  should  therefore  be  called  a  Wiener-bevy 
process  or  Brownian  motion.  In  this  approach  the  disturbances  are  assumed  to 
be  the  integral  of  white  noise;  i.e.,  for  this  problem. 


dN' 

dt' 


(56) 


and 


dY' 

dt' 


(57) 


Note  that  now  it  is  only  necessary  to  estimate  the  power  spectral  density  Q 
of  the  disturbance  white  noise  w  =  [wN,wy)T.  Kwatny'*1,  has  used  this  modeling 
approach  to  develop  state  variable  controllers  with  integt.ii  control  proport ie 
Balchen,  et  a]1*7  have  used  a  more  general  "environment 44  disturbance  model. 


w 


Ew  +  n  , 


(58) 


whore  n  is  white  noise.  They  then  consider  the  random  walk  model  as  a  special 
case  of  this  approach;  i.e.  E=0.  We  have  chosen  to  investigate  the  effective¬ 
ness  of  this  disturbance  modeling  approach  in  this  work. 


Using  eq.  (56)  and  (57),  the  state  vector  x  can  be  augmented  with  N' 
and  Y'  so  that  eq.  (11)  now  becomes, 
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x'  =  Fex  +  Geu  +  rew  .  (60) 

In  this  formulation,  c  is  a  small  quantity  introduced  so  that  Potter's  algorithm 
can  be  used  to  design  the  steady-state  controller  and  the  steady-state  filter 
which  estimates  x_' .  With  e=0,  matrix  M  defined  by  eq.  (36)  with  Fe  and  Ge  used 
in  lieu  of  Fg  and  Gs,  respectively,  has  four  zero  eigenvalues;  matrix  M' 
defined  by  eq.  (50)  has  two  zero  eigenvalues.  HacFarlane42  has  shown  that 
under  these  conditions  general  solutions  for  and  P^  in  terms  of  the  eigen¬ 
vectors  of  M  and  M’,  respectively,  do  not  exist.  As  c-*-0 ,  however,  the  steady- 
state  controller  gain  matrix  C  and  the  steady-state  filter  gain  matrix  K 
approach  limits.  As  an  example,  in  the  design  presented  in  detail  below  the 
gain  C3g=995.32  when  e=10-2,  998.50  when  £=10~4,  and  998.54  when  e=10~8.  We 
have,  therefore,  used  e=10~4  here.  It  is  not  necessary  to  include  c/0  in 
the  implementation  of  the  filter  once  the  gains  have  been  obtained. 

The  remaining  definition  needed  to  specify  the  control  loop  design  is  the 
white  noise  w  power  spectral  density  matrix  Q.  General  results  for  Brownian 
motion  can  be  used  to  obtain  a  direct,  physical  way  to  establish  Q.  If  we 
have  a  continuous,  Brownian  motion  process  y(t),  it  has  a  zero  mean  and  the 
following  properties,48 

E[y(t>2]  =  o2t  ,  (61) 

where  a2  is  the  diffusion  coefficient  or  variance  parameter  and, 

E[  (y(t2)-y(t1))2]  =  0 2  ( t ^ — t x )  ,  (62) 

provided  t2>tj.  The  formal  derivative  of  the  Brownian  motion  is  white  noise 
where  the  noise  spectral  density  is  given  by, 
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q  =  a- 


T 


Now  if  a  white  noise  process  is  approximated  by  an  equivalent  white  noise 
sequence  with  time  step  At,  Jazwinski8 8  shows  that  it  must  have  a  covariance 
q  which  is  related  to  the  process  spectral  density  q  by, 


in  order  to  provide  a  constant  noise  "power"  level;  i.e.,  constant  area  under 
the  autocorrelation  curve.  If  At  is  taken  as  t^-t  ,  eqs (62) , (63) ,and  (64)  yield. 


"f  y(t2)-y(t1)' 

2~m 

LI  t2-ti  J 

(65) 


Further,  if  At  is  viewed  as  the  integration  steps ize  used  in  the  Euler 
integration  discretization  of  the  continuous  process  we  can  use  eq.  (64)  and 
(65)  to  approximate  the  white  noise  power  spectral  density  as, 

q  =  At  S2  ,  (66) 

where  S  is  the  maximum  expected  rate  of  change  or  slope  of  the  continuous 
process  being  modeled  by  the  Brownian  motion. 

For  the  design  presented  here,  eq.  (66)  was  used  to  establish  the  distur¬ 
bance  model  power  spectral  densities.  Assuming  that  the  yawing  moment  and 
lateral  force  disturbances  could  change  by  approximately  twice  the  values 
produced  by  a  steady  one  knot  lateral  current  in  one  ship  length  wo  have 
S^=0.002  and  $^=0,005.  With  our  At=0.005,  cq.  ((>6)  then  yields, 

q  =  AtS  2  =  2.000' 10"H  ,  (67) 

11  N 

and , 

q22  =  itSy2  =  1 • 250* 10~7  ■  (68) 


The  measurement  noise  power  spectral  density  and  remaining  assumptions  are 
the  same  as  in  the  previous  section.  With  those  assumptions,  the  control 
loop  was  designed  for  the  Tokyo  Maru  using  the  OPTSYS  program.  The  resulting 
controller  gains  and  Kalman-Bucy  filter  gains  are  shown  in  Table  5.  These 
results  can  be  compared  with  those  shown  in  Table  4  for  t he  design  using  the 
first-order  shaping  filter  disturbance  models.  The  state  feedback  gains 
Cj  j  through 


are  the  same  in  both  approaches . 


T 

controller  gams  C* 

Kalman 

-Bucy  filter  gains  K 

5.5421 

1 

4.6600 

0.8798 

-0.0019 

2.6601 

19.3843 

14.1945 

-0.3561 

6.3895 

3.7193 

3.4862 

-8.1031 

2.4252 

-0.0651 

-0.5677 

4.0000 

-0.8499 

0.0000 

0.0000 

0.0000 

998.50 

0.1358 

0.2366 

0.0907 

-12.571 

-0.1551 

-0.2844 

0.4719 

Table  5.  Optimal  Gains  for  Tokyo  Maru  at  H/T=1.89  and  Fn=0.116  - 
Random  Walk  Disturbance  Model 

To  illustrate  the  effectiveness  of  the  design  given  in  Table  5,  we 
simulated  the  Tokyo  Maru  under  the  control  of  this  controller  using  a  one- 
half  beam  initial  offset  with  all  other  states  zero  and  x(0)=x(0) .  This 
simulation  also  includes  a  ship  passing  at  t'=20  producing  a  disturbance 
four  times  that  shown  in  Fig.  3.  Figure  11  shows  the  resulting  lateral  offset 
response  and  is  directly  comparable  with  Fig.  6.  The  overshoot  at  t'=5.6 
is  1.7  m. ;  the  maximum  deviation  due  to  the  passing  ship  is  3.5  m.  Figure 
12  shows  the  corresponding  rudder  angle  response  and  is  directly  comparable 
with  Fig.  7.  The  controller  provides  control  comparable  to  that  provided  by 
the  first-order  shaping  filter  design  shown  in  Fig.  6  and  7.  The  maximum 
rudder  angle  in  Fig.  12  of  35°  shows  that  the  assumed  disturbance  of  four 
times  the  values  shown  in  Fig.  3  actually  violates  the  validity  of  the  J inear 
model.  To  illustrate  the  integral  control  properties  of  this  design,  we 
simulated  the  Tokyo  Maru  under  the  control  of  this  controller  while  being 
subjected  to  the  design  lateral  current  disturbance  shown  in  Fig.  4.  The 
initial  condition  is  x=x=0  in  this  simulation.  The  lateral  offset  response 
is  shown  in  Fig.  13  which  is  directly  comparable  with  Fig.  8.  The  maximum 
lateral  offset  during  the  startup  of  the  system  is  34.5  m. ;  the  controller 
tl  ’n  returns  the  ship  to  the  desired  track  in  the  presence  of  the  constant  one 
knot  lateral  current.  The  maximum  offset  during  the  ramp  current  change  is 
about  9  m.  The  corresponding  rudder  angle  response  is  shown  in  Fig.  14;  the 
Kalman  filter  estimate  of  the  yawing  moment  disturbance  N'  is  shown  in  Fiq.  15. 
Comparing  Fig.  10  with  Fig.  4,  the  filter  produces  essentially  exact  estimates 
of  the  constant  disturbance  levels  of  0. 00103  and  0.000515.  As  desjjed,  this 
design  provides  effective  control  in  the  presence  of  constant  external  distur¬ 
bances  without  a  steady  offset  from  the  desired  path. 
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11.  Lateral  Offset  Response  to  B/2  Initial  Offset  and  Passing  Shi; 
Random  Walk  Disturbance  Model  Design 
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Figure  14. 


As  a  final  illustration  of  the  rfi'cri  in.:  1  inn  t  it  ion.*;  ot  tin 

controller  given  in  Table  S,  we  simulated  the  Tokiio  Vara  etit  er  mg  tlu'  ABF 
harbor  defined  in  Fig.  S  while  under  control  of  tin:  controller .  The  resulting 
lateral  offset  response  in  shown  in  Fin.  le.  The  miximum  ievi.it  ion  of  .f-  m. 
is  not  excessive  but  the  corresponding  i  uddei  angle  response  shown  in  Fi  ].  1  7  i 
physically  impossible.  The  controllers  developed  here  lust  cannot  n-comm. >dat e 
this  extreme  a  disturbance.  With  the  lateral  current  1"  teasing  from  knot  s 
to  zero  in  just  over  one  ship  length  an  extremely  largo  yawing  moment  t  ■  ;  1  ice  • 
on  the  ship.  This  is  very  difficult  for  a  nonunr ie i putory  (causal)  control 
system,  as  used  here,  to  handle.  a  pilot  won  1  1  most,  likely  create  a  y.iwin  i  rat 
just  prior  to  t  lie  harbor  entrance  in  anticipation  of  tile  1  at  oe  yawing  moment 
which  he  knows  will  be  experienced  at  the  channel  mouth.  A  control  system 
could  be  created  in  a  similar  manner  using  a  prept  ogr  immed  rudder  control 
based  on  time-to-go  to  the  channel  mouth.  This  "local"  control  feature  could 
be  superposed  upon  a  general  control  law  of  the  type  •, tailed  here. 
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Figure  16.  Lateral  Offset  Response  at  ABC  Harbor  Entrance 
Random  Walk  Disturbance  Model  Design 
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3.4  Necessity  for  an  Adaptive  Design. 


In  the  previous  section  we  showed  that  the  path  controller  designed  using 
random  walk  disturbance  models  provided  effective  control  with  essentially 
zero  mean  disturbances  as  shown  in  Fiqure  11  and  with  the  design  lateral 
current  disturbance  as  shown  in  Fig.  13.  In  these  examples,  the  controller 
was  designed  for  the  ship  characteristics  at  its  least  course  stable  depth-to- 
draft  ratio,  H/T=1.89,  and  then  the  ship  was  simulated  to  operate  at  this  same 
depth-to-draf t  ratio.  The  characteristics  of  the  ship  are,  however,  known  to 
change  considerably  with  water  depth  as  shown  in  Table  2.  They  will  also 
change  with  ship  operating  condition  such  as  draft  and  trim.  We  have 
previously  shown  that  if  a  constant  gain  design  is  to  be  used,  the  best  over¬ 
all  performance  is  achieved  by  using  the  design  developed  for  the  ship's  least 
course  stable  depth-to-draf t  ratio.1  With  this  design  approach,  the  question 
remains  as  to  how  much  performance  is  lost  by  not  adjusting  the  gains  to  be 
optimal  at  each  ship  operating  condition.  Notice  also  that  the  model  of  the 
ship  (Fe)  included  in  the  Kalman  filter  implementation  in  accordance  with 
eq.  (47)  should  be  updated  to  reflect  the  ship  characteristics  in  each 
operating  condition  if  these  are  known.  An  adaptive  path  controller  is  one 
which  can  detect  the  existing  ship)  characteristics  and  automatically  adjust 
the  controller  and  filter  gains  and  the  Kalman  filter  model  of  the  ship 
accordingly . 

To  illustrate  the  effect  of  an  "incorrect"  water  depth,  we  simulated  the 
Tokyo  Mar u  under  the  control  of  the  controller  defined  in  Table  5  while  being 
subjected  to  the  design  lateral  current  disturbance  shown  in  Fig.  4  in  deep 
water  (H/T=*>) .  The  resulting  lateral  offset  response  is  shown  in  Fig.  18  as 
the  solid  line.  This  response  is  comparable  with  that  shown  in  Fig.  1  i  for  the 
same  disturbance  in  the  "design"  water  depth  H/T=1.89.  Also  shown  in  Fig.  18 
as  the  dashed  line  is  the  response  which  would  result  if  the  controller  had 
been  redesigned  (or  adapted)  to  be  optimal  at  H/T  -  in  lieu  of  H,/T=1.89.  The 
loss  in  performance  due  to  the  use  of  a  nonadaptive,  constant-gain  and  constant 
filter  ship  model  design  can  be  seen  by  comparing  the  two  responses  in  Fig.  18. 
The  maximum  lateral  offset  during  the  startup  of  the  system  increases  from 
50.2  m.  to  63.5  m. ;  the  maximum  offset  during  the  ramp  change  in  disturbance 
level  increases  from  13.8  m.  to  18.5  m.  While  the  transient  response  is 
degraded,  the  controller  does  retain  the  desired  inteqtul  control  propet ty  and 
returns  the  ship  to  the  desired  path  in  the  presence  ot  a  constant,  disturbance. 


The  rudder  anqle  response  correspondinq  to  Fiq.  18  is  shown  in  Fig.  19;  the 
yawing  moment  estimate  N'  is  shown  in  Kiq.  20.  Comparinq  Fiq.  20  with  Fiq. 

4,  the  filter  produces  incorrect  estimates  of  the  constant  disturbances.  The 
effect  of  using  nonoptimal  filter  qains  and  incorrect  characteristics  in  the 
Kalman  filter  ship  model  Fc  is  to  produce  an  equilibrium,  mean  yawing  moment 
estimate  of  about  0.00199  durinq  the  first  15  ship  lengths  in  lieu  of  the 
0.00103  shown  in  Fig.  4.  This  bias  error  in  the  yawing  moment  estimate  and  a 
similar  error  in  the  lateral  force  estimate  degrade  the  transient  performance 
but  the  controller  has  the  very  desirable  property  of  producing  the  same 
equilibrium,  mean  state  as  would  the  optimal  controller.  For  example,  both 
the  H/T=1.89  design  and  the  optimal  H/T=°°  design  yield. 
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-  H/T=1.R9  - H/T-'” 

Figure  18.  Lateral  Offset  Response  to  Design  Lateral  Current  - 
Designs  for  H/T=1.R9  and  H/T~‘"  Operating  at  H/T="'. 


Figure  20.  Yawing  Moment  Estimate  with  Design  Lateral  Current. 
Design  for  H/T=1.89  Operat  ing  at  H/T='". 
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when  the  ship  is  subjected  to  the  yawing  moment  and  lateral  force  shown  in 
Fig.  4  for  t'<15  while  in  deep  water  (H/T=<”) .  The  equilibrium  rudder  angle  is 
nonzero  here  because  the  "design"  disturbance  in  Fig.  4  is  not  exactly  equiva¬ 
lent  to  a  lateral  current  when  applied  to  the  ship  in  any  water  depth  other 
than  the  H/T=1.89  for  which  it  was  derived. 

For  a  path  controller  to  remain  optimal  at  each  operating  condition,  the 
controller  gains  C,  the  Kalman-Bucy  filter  gains  K,  and  the  ship  model  Fe  used 
in  the  filter  n q.  (47)  must  be  adjusted.  The  change  of  the  ship  characteris¬ 
tics  Fs  or  Fe  with  water  depth  is  shown  in  Table  2.  To  illustrate  the  effect 
of  water  depth  on  the  gains.  Table  f>  shows  t.ho  optimal  gains  for  the  Tokyo 
Marti  at  H/T=°°.  These  gains  produce  the  lateral  offset  response  shown  by  the 
dashed  line  in  Fig.  18.  The  gains  in  Table  6  for  H/T=™  can  be  compared  with 
those  in  Table  5  for  H/T=1.89.  The  controller  feedback  gains  on  '  and  Y' 
show  the  greatest  variations. 


T 

controller  gains  C 

Kalman 

-Bucy  filter 

gains  K 

5.8752 

4.6624 

0.8842 

0.0001 

2.4635 

19.4819 

15.1 516 

-0.2571 

1.9204 

3.5844 

2.8169 

-10.9043 

2.4252 

0.0034 

-0.4099 

4.6587 

-0.7955 

0.0000 

0.0000 

0.0000 

1126.44 

0.1394 

0.2548 

0.0509 

82.614 

-0.0878 

-0.1596 

0.5079 

Table  6.  Optimal  Gains  for  Tokyo  Maru  at  H/T= '  and  Fn=0.116  - 
Random  Walk  Disturbance  Model 

In  Table  2,  it  can  be  seen  that  ten  parameters  in  the  linear  state 
variable  equations  of  motion  of  the  ship  eq.  (11)  vary  with  water  depth. 

These  same  parameters  will  also  vary  with  draft,  trim,  etc.  In  considering 
the  use  of  an  adaptive  control  system  it  is  useful  to  establish  how  important 
the  variation  in  each  of  those  parameters  is  to  the  response  of  the  system. 

If  the  variation  of  a  parameter  has  little  influence  on  or  correlation  with 
the  response,  it  will  be  difficult  to  establish  that  parameter  from  the 
response.  Fortunately,  if  a  parameter  has  little  influence  on  the  response 
it  may  not  be  necessary  to  estimate  the  parameter  accurately  or  it  may  not  be 
necessary  to  adapt  for  the  changes  in  that  parameter  at  all.  In  our  previous 
work,*  we  studied  the  sensitivity  of  the  ship  path  control  response  to  the  ten 


parameters  which  vary  with  depth.  This  study  was  for  somewhat  di t  for  ant  path 
controller  designs  for  the  Tokyo  Maru  utilizing  first-order  shaping  filter 
disturbance  models.  These  results  are  still  useful,  however,  in  this  context 
and  are  summarized  in  Table  7.  In  this  computer  experiment,  a  number  of 
controllers  were  designed  and  then  their  root  mean  square  (RMS)  responses  to 
the  design  disturbance  and  noise  levels  were  established  using  the  OPTSYS 
computer  program.  These  results  were  then  combined  into  a  single  RMS  cost 
J  which  is  proportional  to  the  cost  function  used  in  the  optimal  control 
system  design;  i.e., 

J  =  p'2  +  5.881B  n,?  +  iS ' 2  +  ,  (70) 

where  signifies  the  RMS  value  of  each  quantity.  This  RMS  cost  provides 

an  appropriate  single  measurement  of  the  effectiveness  of  each  controller. 

case 

1.  all  parameters  adapted;  optima]  at  H/T=» 

2.  all  parameters  adapted  except  y^  (  ,  f  ?t. ,  >  ,  or  f 
|  3.  all  parameters  adapted  except  f  , ,  or 

4.  all  parameters  adapted  except  f 

5 .  all  parameters  adapted  except  y 

6.  all  parameters  adapted  except:  f 

7.  all  parameters  adapted  except  f , 

•1.  no  paramete  rs  adapted;  optimal  tor  ll/'t  d.HU 

9.  only  parameters  f  ,  f.jr  ,  and  , adapt  ed 

all  nonadapted  parameters  for  Il/T  l.  ei 

Table  7.  RMS  Cost  J  for  Tokyo  Maru  with  Various  Connollers  at  H/T-'j  - 
Kirst-Order  Shaping  filter  Disturbance  Model' 

Shown  in  linos  1  and  R  of  Table  7  are  tin  RMS  cost  J  for  the  controllers 
designed  to  be  optimal  at  M/T="  and  H/'T  1  .  M9 ,  respectively,  when  both  are 
operated  at  H/T=iV  .  The  RMS  cost  can  be  seen  to  degrade  t  tom  the  optimal 
9.00681  to  0.0978  t  (  +  1S1.)  if  a  noriadapt  ive  design  were  used  in  deep  water 
with  the  design  developed  for  the  ship's  loast-course-st able  depth  of  H/T-1.89. 
To  show  the  importance  of  each  individual  parameter  we  also  designed  a  series 
of  controllers  using  all  character istics  for  H/T=®  except  a  single  parameter 
which  was  set  at  its  value  for  H/T=1.H9.  This  product's  the  response  of  a 
controller  which  could  adapt  for  all  parameters  except  the  one  hold  at  its 
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H/T=1.89  value  and  shows  the  sensitivity  of  the  response  to  that  parameter. 


In  lines  2  through  7  of  Table  7,  it  can  be  seen  that  only  f^,  ^25'  anc*  Y21 
have  a  significant  effect  on  the  RMS  response.  Coefficient  f  of  8‘  in  the 
r1  equation  is  by  far  the  most  significant  parameter.  In  line  9  of  Table  7 
we  show  the  RMS  cost  produced  by  a  desiqn  which  would  adapt  only  for  changes 
in  f  ^25'  aiK*  ^21  keeping  the  remaining  seven  parameters  fixed  at 

their  value  at  H/T=1.89.  This  RMS  cost  is  within  0.31.  of  that  achieved  with 
the  optimal  design  shown  in  line  1.  These  results  arc  very  significant  and 
show  that  it  is  not  essential  that  the  adaptive  path  control  systems  to  be 
studied  in  the  following  chapters  adapt  for  all  ten  of  the  parameters  which 
vary  with  water  depth,  draft,  trim,  etc. 

.  To  further  investigate  the  point  made  in  the  previous  paragraph,  we 
designed  a  "partially-adapted"  controller  using  all  the  assumptions  of  Section 
3.3  except  that  the  parameters  and  ^25  were  taken  for  H/T=°"  from  Table  2 
while  the  remaining  eight  parameters  were  given  their  values  at  H/T=1.89.  We 
chose  not  to  include  Y21=^26  ^23  and  ^25  because  this  parameter  appears 

as  a  product  with  N1  in  eq.  (59)  and  is  thus  very  difficult  to  establish 
on-line.  The  resulting  controller  and  Kalman-Bucy  filter  gains  are  shown  in 
Table  8.  These  gains  can  be  compared  with  those  designed  to  be  optimal  at 
H/T=l . 89  and  shown  in  Table  5  and  with  those  designed  to  be  optimal  at  H/T=“> 
and  shown  in  Table  6.  They  can  be  seen  to  approach  the  gains  shown  in  Table  6. 
The  gains  in  Table  8  represent  those  of  a  controller  which  would  adapt  only 
for  changes  in  and  f  at  a  time  when  the  ship  is  operating  in  deep  water. 


f and  f25  for  H/T=» ;  other  f  for  H/T=1.B9 

Table  8.  Gains  for  Partial ly- Adapted  Design  for  Tokyo  Maru  at  Fn=0.11fi  - 

Random  Walk  Disturbance  Model 
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To  illustrate  the  effectiveness  of  the  partially-adapted  design  given  in 
Table  8,  we  simulated  the  Tokyo  Marti  under  the  control  of  this  controller 
while  being  subjected  to  the  design  lateral  current  disturbance  shown  in  Fig. 

4  in  deep  water  (H/T=“) .  The  resulting  lateral  offset  response  is  shown  in 
Fig.  21  as  the  solid  line;  this  is  directly  comparable  to  the  response  shown 
for  the  nonadapted  design  in  Fig.  18.  The  response  with  the  optimal  design 
for  H/T=<»,  the  fully-adapted  design, is  shown  in  Fig.  21  as  the  dashed  line 
for  comparison.  The  maximum  lateral  offset  during  the  startup  of  the  system 
is  53.2  m.  compared  to  50.2  m.  with  full  adaptation  and  63.5  m.  with  no 
adaptation.  The  maximum  offset  during  the  ramp  change  in  disturbance  level 
is  13.9  m.  compared  to  13.8  m.  with  full  adaptation  and  18.5  m.  with  no 
adaptation.  Adapting  only  for  changes  in  f  and  f  thus  provides  transient 
response  very  close  to  that  provided  by  the  fully-adapted,  optimal  design. 

The  rudder  angle  response  corresponding  to  Fig.  21  is  shown  in  Fig.  22;  the 
yawing  moment  estimate  N '  is  shown  in  Fig.  23.  Notice  that  the  equilibrium, 
mean  yawing  moment  estimate  during  the  first  15  ship  lengths  is  about  0.00126 
compared  with  0.00199  with  the  nonadapted  design  in  Fig.  20  and  the  exact 
value  of  0.00103  which  would  be  produced  by  the  ful ly-adapted ,  optimal  design. 


-  adapted  only  for  f  ,  ,  and  f  ,,  - optimal  at  ll/T=- • 

Figure  21.  Lateral  nffset  Kespon;;.  to  Design  Lateral  Current  - 
Partially-Adapted  Design  anil  Design  for  H/T- Operating  at:  ll/T~  ■  . 
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Figure  22.  Rudder  Angle  Response  to  Design  Lateral  Current 
Partially-Adapted  Design  Operating  at  H/T=”r. 
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Figure  23.  Yawing  Moment  Estimate  with  Design  lateral  Current 
Partially-Adapted  Design  Operating  at  h/T--™>. 


Summarizing  this  chapter,  wo  have  shown  that  effective  ship  path 
controllers  can  be  developed  using  a  random  walk  model  for  the  yawing  moment 
and  lateral  force  disturbances.  This  controller  provides  effective  control 
with  short-term,  essentially  zero-mean  disturbances  as  shown  in  Fig.  11  and 
with  more  long-term,  nonzero-mean  disturbances  as  shown  in  Fig.  13.  We  then 
showed  that  this  controller  experiences  a  degradation  in  transient  response 
when  used  at  an  "incorrect  depth"  or  in  any  other  condition  causing  the  ship 
characteristics  to  change  from  the  values  used  in  the  controller  design. 

At  an  "incorrect  depth"  this  controller  is,  however,  capable  of  producing  the 
same  equilibrium,  mean  response  as  would  the  optimal  design.  Finally,  we 
showed  that  adapting  for  changes  in  as.  few  as  two  parameters  in  the  equations 
of  motion  can  eliminate  almost  all  of  the  degredation  in  transient  response 
experienced  with  changing  ship  characteristics.  In  the  following  two  chapters, 
we  investigate  two  different  approaches  for  the  design  of  such  adaptive  ship 
path  controllers  and  evaluate  their  performance  capabi 1 i t ies  and  limitations. 


4.  Weighted  Least-Squares  Parameter  Estimator 

In  this  section,  we  present  the  derivation  and  evaluation  of  a  Weighted 
Least-Squares  (WLS)  parameter  estimator  which  could  be  used  in  the  gain  up¬ 
date  loop  of  an  adaptive  ship  path  control  system.  This  section  is  based 
upon  Cuong's  Ph.D.  Dissertation.  The  WI.S  approach  requires  no  probabilistic 
assumptions  about  the  uncertainties  of  the  measurements  and  the  process  dis¬ 
turbances.  It  treats  the  observed  measurement  data  in  a  deterministic  rather 
than  probabilistic  sense.  The  unknown  system  parameters  are  chosen  so  that 
the  response  of  a  model  of  the  ship  dynamics  provides  a  "best  fit”  to  a  finite 
record  of  sensor  measurements.  This  treatment  of  measurement  data  is  called 

“limited  memory"  or  "moving  window".  The  moving  window  concept  was  originally 

48 

proposed  by  Jazwinski  to  handle  the  divergence  problem  of  Kalman  filters. 
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Dunn  and  Montgomery  utilized  Jazwinski ' s  concept  in  a  batch  processing 
manner  to  do  the  on-line  parameter  identification  for  the  NASA  F8-DFBW  air¬ 
craft.  This  approach  is  conceptually  simple  which  made  it  a  logical  first 
choice  in  our  parameter  identification  investigation. 


Derivation  and  Development. 


The  source  of  information  for  the  parameter  estimation  algorithm  is 
the  measurement  vector  .  Sensor  outputs  are  sampled  and  the  results  are 
buffered  in  the  memory  of  a  control  computer.  After  LW  samples  have  been 

collected,  the  parameter  estimation  algorithm  starts  processing  this  infor¬ 

mation.  Meanwhile,  additional  measurement  data  can  be  ignored  until  it  is 
time  to  buffer  another  window  of  data.  The  timing  sequence  is  illustrated 
in  Fig.  24  where  each  graduated  mark  on  the  time  axis  represents  a  sampling 
instant  and: 

•  LU  is  the  length  of  the  parameter  estimation  cycle, 

•  LW  is  the  length  of  the  data  window, 

•  LB  is  the  length  of  the  batch  update  calculation  period. 

A  cycle  of  a  typical  parameter  estimation  process  start:;  wi t  h  a  stabiliza¬ 
tion  period  during  which  the  ship  is  run  under  the  control  of  recently 
acquired  parameter  estimates  and  gain  matrices  C  and  K.  This  period  allows 
transients  caused  by  the  abrupt  changes  in  the  parameter  estimates  and  gain 
matrices  to  die  out.  This  period  is  followed  by  the  data  window  loading 


period  in  which  LW  samples  of  the  measurement  vector  z_'  are  collected  and 
stored  in  the  estimator  memory.  The  cycle  is  completed  by  the  batch  update 
calculation  period  in  which  the  parameter  estimation  and  gain  updates  are 
performed.  This  period  produces  a  new  set  of  parameter  estimates  and  new 
values  for  the  control  gain  C  and  the  Kalman  filter  gain  K  if  the  parameter 
estimates  are  judged  to  be  physically  reasonable. 


new  parameter  estimates  £  and  gains 

C  and  K  are  incorporated  into  _ 

controller 


.  I  time 

batch  1 


stabilization  | 


window 

loading 


update 
calculat ion 


Figure  24.  Timing  Sequence  for  Weighted  Least-Squares  Parameter  Estimator. 

After  another  stabilization  period,  a  new  data  batch  is  buffered  and 
then  processed  by  the  algor ithm.  In  this  manner  L.W  memory  locations  per 
measurement  channel  are  needed  for  the  moving  window.  When  the  states 
change  very  slowly  or  the  time  constants  of  the  system  are  long,  the  window 
has  to  bo  very  long.  This  can  represent  a  large  memory  requirement.  The 
parameter  estimation  task  can,  however,  be  processed  on  a  time-available 
basis. 

In  the  control  loop  development,  tiie  equations  of  motion  of  the  ship 
were  expressed  by  the  5-state  system  eq.  (11).  The  rudder  angle  A  was  in¬ 
cluded  as  a  state  in  order  to  introduce  a  realistic  t imo  response  for  the 
steering  system.  We  noted  in  Section  1.2  ttiat  the  rudder  angle  time  response 

is  known  exactly  if  the  rudder  time  constant,  the  initial  condition  '  ,  and 

o 

the  rudder  command  history  are  known.  Therefore,  A(t)  will  be  considered  a 
known  function  in  this  development.  The  system  equations  then  become, 
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The  disturbances  N'  and  Y'  are  considered  to  be  stochastic  processes  in 

eq.  (71)  but  this  assumption  has  no  effect  on  the  WLS  development. 
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Cuong  has  shown  that  for  a  WLS  parameter  estimation  algorithm  to 
be  effective  in  estimating  ship  parameters  all  states  must  be  available 
for  measurement.  The  measurement  vector  to  be  used  in  the  gain  update 
loop  is  therefore  taken  as, 

z'  =  x  +  v' ,  (73) 


where  v'  is  a  Gaussian  white  measurement  noise  vector  which  is  ignored 
within  the  WLS  algorithm.  In  general,  the  parameter  estimation  loop 
requires  a  more  accurate  level  of  measurements  than  does  the  control  loop. 
In  the  simulations  presented  later  in  this  section,  the  reference  level  of 
the  measurement  noise  vector  v'  is  as  shown  in  Table  9.  This  represents  a 
very  small  level  of  noise  but  is  not  necessarily  the  highest  level  which 
could  be  tolerated. 
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Table  9.  Reference  Measurement  Noise  for  Parameter  Lstimator  v' 


It  was  demonstrated  in  Section  3.4,  that  accurate  estimates  of  only 
two  of  the  10  coefficients  in  eq.  (71),  and  f  ,  are  essential  to  an 

effective  adaptive  control  loop.  Our  first  approach  was  therefore  to  select 
a  parameter  vector  of  minimum  dimension,  i.e., 


(P  /  P  »  P  /  P.]T=  ( f  2  3 '  f  25 '  N’ '  y')T' 


4  q 


in  order  to  minimize  the  computational  load.  Cuong  has  shown  that  this 
parameter  choice  is  a  condition  of  underparameterization  in  which  the  WLS 
algorithm  converges  to  an  estimate  of  £  which  is  completely  impossible 
physically.  Even  though  accurate  estimates  of  the  other  coefficients  in 
the  open-loop  dynamics  matrix  f^  are  not  really  needed  by  the  control  loop, 
these  coefficients  must  be  included  in  the  parameter  set  to  be  estimated 
by  the  parameter  estimator.  If  they  are  not,  the  effects  of  errors  in  these 
coefficients  are  interpreted  as  changes  in  f^  and  f  .  The  parameter  vec¬ 
tor  to  be  estimated  by  the  parameter  estimator  must  therefore  be  taken  as, 
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Recall  that  the  disturbances  N"  and  Y*  art1  estimated  by  the  control  loop 
Kalman  filter  and  these  estimates  are  actively  used  for  control  purposes. 

The  unknown  disturbances  must  also  be  included  in  the  parameter  set.  If 
they  are  not,  the  effect  of  the  disturbances  is  interpreted  as  changes  in 
the  system  parameters.  The  disturbance  estimates  produced  by  the  parameter 
estimator  are  only  a  by-product  of  the  estimation  process  and  are  not  used 
further . 


For  use  in  the  WLS  parameter  estimator,  the  ship  is  modeled  by  the 

T 

state  vector  x  =  [i K  ,  r'  ,  h'  ,  n'  1  and  the  following  equations: 
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where  the  values  of  the  coefficients  in  the  matrix  are  assigned  constant 
assumed  values  indicated  by  (...).  These  "known"  coefficients  are  assigned 
their  value  at  the  ship's  least-course-stable  depth-to-drai t  ratio  H/T=1.89. 


With  p  =  N'  and  pu 

values  of  M '  and  Y' . 

the  y.  coefficients 
i] 

of  N'  and  Y'  include 
in  the  values  of  the 
acceptable . 


=  Y' ,  the  parameter  estimator  estimates  the  average 

These  estimates,  however,  are  not  very  accurate  lecause 
are  kept  constant  at  assumed  values.  The  estimates 
the  effects  of  these  errors.  Since  we  are  not  interested 
disturbances,  per  se ,  this  arrangement  is  clearly 


At  the  beginning  of  the  parameter  estimation,  the  computer  memory  holds 

the  system  states  at  the  start  of  the  data  window  and  the  rudder  angle  S' 

and  measurement  zj  histories  during  the  data  window.  Using  any  set  of  the 

parameters  £,  the  buffered  state  at  the  start  of  the  window,  and  the  buffered 

rudder  angle  history,  the  ship  model  eg.  {71)  is  simulated  through  the  same 

time  period  as  the  window  to  produce  x  ^  where  k  is  the  sample  time  index. 

This  yields  a  set  of  model  measurement?;  Hx  ,  which  fill  the  data  window 

-m,  k 

and  a  Weighted  Least-Squares  cost  function  or  fit  error  .7  can  be  defined  as? 
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The  weighting  matrices  can  be  used  to  adjust  for  dimensional  differences 

among  the  measurements  or  can  be  used  to  place  more  emphasis  on  matching 
selected  measurements  or  on  matching  the  most  recent  measurements.  The 
weighting  matrix  is  assumed  to  be  diagonal  so  the  cost  function  J  can  be 
expanded  to  give, 
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The  WLS  parameter  estimation  problem  i?;  then  to  vary  the  parameter  estimate 
£  to  produce  the  parameter  estimate  vector  £*  which  yields  a  minimum  of  this 
cost  function.  This  causes  the  response  of  the  ship  model  as  represented 
by  the  model  measurements  to  be  a  least-squares  fit  to  the  measured  response. 

The  minimization  of  the  Weighted  Least-Squares  cost  function,  eq .  (77), 
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can  be  accomplished  iteratively  using  any  of  a  number  of  numerical  algor it  inns . 

0 

Dunn  and  Montgomery  used  a  modified  tJewton-kaphson  algorithm.  it  ir  diffi¬ 
cult  and  expensive  in  this  problem,  however,  to  numerically  estimate  the 
gradient  of  J  needed  by  this  algorithm.  This  is  especially  true  if  an  on¬ 
line  application  is  desired.  In  this  situation,  a  dot Lvative-free  algorithm 
is  very  attractive.  It  is  also  desirable  to  make  efficient  use  of  previously 

computed  function  values.  After  preliminary  work  with  a  number  of  alternative 

methods,  the  DUD  (Doesn't  Use  Derivatives)  algorithm  was  selected  for  use 
-  M 

here.  This  algorithm  was  developed  by  Ralston  and  was  presented  by  Ralston 
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and  Jennrich  in  1978.  Cuong  extended  the  DUD  algorithm  to  the  Weighted 
Least-Squares  problem  and  added  additional  search  stopping  conditions  needed 
for  its  application  to  the  ship  parameter  estimation  problem.  The  DUD 
algorithm  does  not  require  the  gradient  of  J  as  the  name  implies.  Once  the 

algorithm  is  started,  only  one  evaluation  of  J  is  needed  in  each  optimiza¬ 

tion  cycle  performed  under  normal  conditions.  DUD  has  proven  to  be  an 
effective  and  efficient  search  algorithm  for  this  problem. 

As  noted  above  the  diagonal  weighting  matrices  in  eq.  (77)  can 
be  used  to  assign  different  levels  of  importance  to  the  various  measure¬ 
ments.  They  can  be  varied  with  the  time  index  k  to  account  for  a  lower 
confidence  in  old  data.  They  can  also  be  used  to  achieve  a  non-dimensional- 
ization  of  the  various  errors  in  fitting  the  model  response  to  the  data. 

In  our  work,  the  W,  were  chosen  to  be  constant  with  time.  Recall  that  in 
k 

the  solution  of  the  Linear-Quadratic-Gaussian  state  estimation  problem, 
the  innovation  is  "weighted"  by  R  \  where  R  is  the  spectral  density  of 
measurement  noise.  It  is  reasonable  to  use  a  parallel  approach  here  and 
select  W^=W=R  \  This  accomplishes  a  desirable  non-dimensionalization 
and  reflects  the  relative  noise  l.vels  in  each  measurement.  Using  this 
approach,  the  DUD  algorithm  was  found  to  require  about  171,  fewer  evalua¬ 
tions  of  J  to  converge  compared  with  runs  using  an  equal  weighting;  i.e., 

I  where  I  is  the  identity  matrix. 

The  simulation  of  the  ship  model,  eg.  (75),  begins  from  tin'  buffered 
state  of  the  system  at  the  start  of  the  data  window  loading  period.  We  have 
found  that  when  realistic  noise  exists  in  the  system,  it  is  proforrablo  to 
use  the  estimate  of  the  state  at  this  point,  produced  by  the  control  loop 
Kalman  filter  rather  than  the  noisy  measured  states.  The  state  estimate 


is  an  expected  value  of  the  state  at  the  start  of  the  window  loading  period 
and  is  not  strongly  affected  by  the  noise  at  that  time  which  probabilisti¬ 
cally  might  be  large.  It  is  undesirable  in  general,  however,  to  use  fil¬ 
tered  data  in  parameter  estimation  so  all  of  the  data  in  the  moving  window, 
except  the  initial  values,  are  taken  directly  from  sensor  measurements. 

This  ensures  the  maximum  information  content  in  the  data  base. 

We  noted  above  that  the  parameter  estimates  £*  or  £  produced  by  the 
DUD  algorithm  in  the  WLS  parameter  estimator  are  only  used  to  calculate 
updated  controller  gains  C  and  Kalman  filter  gains  K  if  the  estimates 
are  judged  to  be  physically  reasonable.  In  our  work,  we  performed  this 
validity  check  by  noting  the  expected  range  of  each  parameter  given  in 
Table  2  and  the  sensitivity  study  results  discussed  in  Section  3.4  which 
showed  that  the  control  loop  needs  accurate  estimates  of  only  the  two 
parameters, f  and  f  The  validity  check,  therefore,  utilizes  the  fol¬ 

lowing  logic  : 

•  if  either  f^^  or  f  is  beyond  the  permissible  range,  as  given 
in  Table  10,  no  update  is  made;  i.e.,  the  controller  remains 
unchanged; 


•  if  f23  and  f  are  within  the  permissible  ranges  but  some  of  the 
remaining  estimates  are  beyond  the  permissible  ranges,  only  f 
and  f^  are  updated  in  eq.  (59) ; 


if  all  estimates  are  within  the  permissible  ranges, 
are  updated  in  eq.  (59) . 


all  the  f. 
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If  an  update  is  to  be  made,  the  new  feedback  gains  C  and  filter  gains  K 


are  calculated  as  described  in  Section  3.3.  The  results  are  then  incor¬ 


porated  into  the  controller. 


Parameter 

Lower  limit  p 

t 

Upper  limit  £^ 

f22 

-2.50 

-1.10 

f23 

2.10 

12.00 

f25 

-1.30 

-0.30 

f32 

0.0 

0.42 

f 

33 

o 

r-"' 

1 

-0.  30 

f35 

-0.22 

0.0 

based  upon  experimental  range  +  303.. 

Table  10.  Permissible  Ranges  of  Parameter  Estimate;-,  for  the  Tokyo  M.un 
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The  development  of  the  WLS  parameter  estimator  is  complete  at  this 
point  except  for  the  addition  of  an  open-loop  excitation  or  input  dither 
signal  which  will  be  shown  to  be  necessary  in  the  next  section.  Before 
proceeding,  however,  it  is  appropriate  to  review  some  of  the  alternative 

parameter  set  choices  which  were  found  to  be  unsuccessful  prior  to  our 

4  'j 

choosing  eq.  (74).  Cuong  discusses  these  in  detail  so  only  brief 
comments  will  be  included  here.  The  following  parameter  choices  were 
considered : 


p7  =  Y21N\  p8  =  y22y\  p9  =  y31n-,  p10 


y  Y  1  ; 

32 


’21'  p8  “  N’ '  P9 


^  32 '  P10  Y  '  r22 


*22'  Y31 


p7  >2l'  p8  Y32'  N  ^  '  Y  Y  '  '22  ^22 '  *31  Y3l' 

In  the  first  approach,  the  r  and  8  equations  included  terms  p  +  p  and 

7  8 

py  +  p  ,  respectively.  The  WLS  algorithm  quickly  produced  physically  im¬ 
possible  results  with  the  terms  in  these  sums  having  very  large,  equal 
magnitudes  but  opposite  sign.  In  the  second  approach,  the  r  and  (1 

equations  involved  products  of  parameters;  i.e.  p  p  and  p  p  ,  respectively. 

7  o  i  l'J 

The  coefficients  y  and  were  given  constant,  assumed  values.  The  DUD 

algorithm  failed  to  converge  in  this  problem.  The  third  approach  would  use 
the  yawing  moment  and  lateral  force  disturbance  estimates  from  the  control 
loop  Kalman  filter  in  the  parameter  estimator.  This  approach  was  not  actually 
tried  but  was  discarded  because  it  would  produce  a  highly  undesirable-  coupling 
between  the  state  estimation  filter  and  the  parameter  estimator. 

4.2  Time- Invariant  Parameters  Performance 


In  this  section,  we  evaluate  the  performance  of  a  WLS  ship  parameter 
estimator  when  the  parameters  are  time-invariant.  With  the  WLS  parameter 
estimator  each  parameter  estimation  cycle  is  a  completely  independent 
process.  The  only  information  used  from  the  previous  cycle  is  the  old 
parameter  estimate  which  is  used  as  the  initial  estimate  in  the  DUD  algo¬ 
rithm  search.  This  initial  estimate  affects  the  convergence  rate  of  the 
search  but  has  negligible  effect  on  the  resulting  parameter  estimate'. 
Therefore,  time- invariant  parameters  in  this  context  requires  that  the 


parameters  be  constant  only  during  each  data  window  loading  period.  The 
parameters  could  vary  from  one  window  loading  period  to  the  next.  The 
stepped  bottom  profile  shown  in  Fig.  25  was  therefore  used  in  our  simula¬ 
tions  to  evaluate  the  effectiveness  of  the  WLS  parameter  estimator  with 
locally  time-invariant  parameters.  The  timing  sequence  was  selected  so 
that  the  window  loading  period  would  not  overlap  one  of  the  discrete  changes 
in  water  depth.  Ship  parameters  can  of  course  change  due  to  many  causes; 
i.e.,  draft,  trim,  water  depth,  etc.  Rather  than  change  the  parameters  in 
a  completely  arbitrary  and  random  manner,  however,  we  chose  water  depth 
as  the  independent  variable  to  use  to  study  the  effect  of  physically  possible 
parameter  changes . 


Figure  25.  Stepped  Bottom  Profile 

The  signal-to-noise  ratio  is  critical  to  the  effectiveness  of  a  WhS 
ship  parameter  estimator.  The  signal  here  is  the  measurement  vector  z. ' 
which  is  contaminated  by  the  noise  v'  which  will  be  at  least  as  large  as 
shown  in  Table  9.  For  the  measurements  to  contain  information  about  the 
dynamic  characteristics  of  the  ship,  it  is  necessary  that  the  ship  trajec¬ 
tory  include  motion  which  is  large  compared  with  the  noise  level.  This 
requirement  is  in  direct  conflict  with  precise  ship  path  control  which 
seeks  to  eliminate  any  motion  from  the  desired  track.  The  information 
content  of  the  data  window  can  be  increased  by  increasing  the  amount  of 
ship  motion  and  by  lengthening  the  data  window.  In  off-line  systems  iden¬ 
tification  of  ship  parameters  the  "data  window"  is  usually  many  minutes 


r 
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in  which  the  ship  travels  many  ship  lengths.  It  is  usually  considered  neces¬ 
sary  for  the  record  length  to  be  at  least  a  few  times  the  longest  time  constant 
of  the  system  for  successful  identification  using  noisy  measurement.  The 
tests  are  conducted  with  constant  ship  characteristics  and  in  unrestricted 
waters  so  that  significant  ship  motion  can  be  used.  In  path  control  in 
restricted  waters,  ship  motion  must  be  restricted,  calculations  must  be  per¬ 
formed  on-line,  and  the  ship  parameters  may  change  significantly  over  even 
a  ship  length.  Parameter  estimation  in  this  context  is  very  difficult  and 
requires  very  low  levels  of  acceptable  noise  and  engineering  tradeoff  in  design. 

The  parameters  f  which  the  WLS  parameter  estimator  must  estimate  are 
the  open-loop  parameters  of  the  ship.  During  path  control,  however,  the 
dynamics  of  the  ship  are  governed  by  the  closed-loop  properties  of  the  ship 
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and  the  control  loop.  Further,  Astrom  and  Wittenmark  have  noted  the  fol¬ 
lowing  conditions  as  necessary  for  the  convergence  of  a  least-squares  para¬ 
meter  estimator: 

•  the  system  must  be  persistently  excited  by  the  input  signal: 

•  the  input  signal  must  be  independent  of  any  disturbances. 

The  estimation  of  open-loop  ship  parameters  in  the  presence  of  disturbance's 
and  the  path  control  are  basically  incompatible  and  cannot  be  undertaken 
simultaneously.  Thus,  an  adaptive  ship  path  controller  must  devote  itself 
alternately  to  parameter  estimation  using  open- loop  rudder  commands  and  then 
to  path  control.  If  the  rudder  commands  are  generated  by  the  feedback  con¬ 
troller  during  the  data  window  loading  period,  the  input  will  be  corre¬ 
lated  with  the  disturbances  and  it  will  not  be  possible  t o  determine  all 
parameters.  Further,  the  controllers  used  here  will  generally  not  allow  the 
level  of  ship  motion  needed  for  the  parameter  estimation.  Adaptive  path  con¬ 
trol  must  therefore  accept  increased  ship  motion  during  the  window  loading 
period.  The  improved  performance  from  tiaving  better  parameter  estimates  must 
offset  the  effects  of  the  periods  of  open-loop  control. 

The  open- loop  input  signals  used  in  adaptive  controllers  art'  often 
called  dither  signals.  They  usually  have  a  zero  mean  to  reduce  the  long 
term  effect  on  the  system.  In  highly  dynamic  systems  such  as  aircraft  , 
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the  dither  signal  could  be  superposed  upon  the  pilot's  input  commands 
and  continuous  small  dither  signals  have  been  proposed.  In  the  ship  path 
control  problem,  the  signal-to-noise  ratio  is  such  that  the  dither  signal 
magnitude  must  be  so  large  that  continuous  use  would  be  unacceptable.  Al¬ 
ternate  periods  of  open-loop  control  using  the  dither  signal  with  the 
control  loop  turned  off  followed  by  closed-loop  control  without  the  dither 
signal  are  used  here.  The  purpose  of  the  dither  signal  is  to  excite  all 
modes  of  the  system  and  to  maintain  them  at  large  enough  amplitudes  for 
the  parameter  estimation  to  be  successful.  The  signal  should  therefore 
include  all  the  open-loop  natural  frequencies  of  the  system. 

In  this  investigation  the  dither  signals  were  chosen  to  be  square 
wave  rudder  commands  of  various  amplitude  and  period.  An  alternative  might 
be  the  Pseudo-Random  Binary  Sequences  (PRBS)  used  in  systems  identification 
of  ship  parameters.  In  general,  the  larger  amplitude  dither  signals 
excite  the  ship  more  and  improve  the  effectiveness  of  the  parameter  esti¬ 
mation.  With  a  physical  limit  on  rudder  angle  of  about  35°,  we  have  used 
a  dither  rudder  command  6^'  magnitude  of  0.5  rad.  in  all  the  simulations 
presented  below.  With  short  dither  periods,  the  rudder  does  not  have  time 
to  reach  this  commanded  value  before  the  sign  of  the  command  is  reversed. 

A  magnitude  of  0.5  rad.  already  stretches  the  validity  of  the  linear 
system  equations  so  higher  values  were  not  considered.  The  initial  dither 
signal  (at  the  beginning  of  the  window  loading  period)  was  given  a  sign 
opposite  to  the  rudder  angle  existing  just  prior  to  that  time  in  order  to 
maximize  its  effect  on  the  ship.  Since  the  square  wave  dither  signal  is 
used  for  only  the  finite  window  loading  period  (not  continuously) ,  it 
can  excite  frequencies  below  the  fundamental  frequency  of  the  signal. 

4.2.1  Effect  of  Data  Window  Length  and  Dither  Period 

Various  simulations  of  the  Tokyo  Maru  were  performed  to  evaluate 
the  performance  of  the  WLS  parameter  estimator.  In  our  earliest  work,  we 
simulated  the  Tokyo  Maru  under  the  control  of  an  adaptive  control  system 
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using  a  WLS  parameter  estimator  when  subjected  to  only  0.001  times  the 
reference  noise  level  given  in  Table  0  while  sailing  over  the  stepped 
bottom  profile  shown  in  Fig.  25.  A  short  data  window  of  only  0.2  ship 
length  was  used  with  this  low  noise  level;  the  parameter  estimation  was 
performed  once  per  ship  length.  Starting  with  an  initial  offset  from 
the  desired  path  of  n'  =  rj‘  =  0. 1638  (1  beam)  and  using  no  dither  signal, 
the  WLS  parameter  estimator  failed  to  work.  Since  the  simulation  started 
with  a  reasonably  large  initial  condition,  there  was  large  ship  motion 
in  the  first  two  ship  lengths  and  the  measurements  thus  contained  a  rea¬ 
sonable  amount  of  information  about  the  system.  The  estimator  converged 
in  the  first  cycle  and  gave  poor  estimates  which  were  in  the  permissible 
ranges  given  in  Table  10.  By  the  second  cycle,  the  ship  motion  had  died 
out  to  a  large  extent  and  there  was  much  less  information  contained  in  the 
response.  As  a  result,  even  though  the  a  1  g  ">r  i  t  lun  search  converged,  the 
estimates  were  outside  the  permissible  rang*  •  By  th.  third  cycle,  the 
path  controller  had  suppressed  most  .  <  *  ;,<•  -  hi;  motion  an  i  the  shij  was 
moving  almost  in  a  straight  line.  7L>  .mail  amount  t  information  which 
could  be  extracted  from  this  mi  id  m  ••••  -.  t  •  .  i  1 ;  !  .  masked  by  the 

measurement  noise  (even  though  ?:!••  .  '  *'•  ;  •  •  •  nee  level).  The 

DUD  algorithm  failed  to  rnnv  t  ;•• .  1 1  •  -  •  ■  m  m  .*  i  •>,  a  longer  data 

window  would  be  of  no  value  in  'hi  •  <  1  1  ■  ‘  g  a  ’.  it.her 

signal  with  a  period  'if  hi;  !•  r.  ••  ••  ’  ••  •  •  •  ,  ,  .•  •  },«  lata  window 

length  to  0. 5  ship  length  (I  1  am  '.-i  .  *  imur.ir 

quickly  converged  to  pst  imat  e  whi  •  w-  * .  ■  oat.  *  •  •  i  •  .  ggt.il  ir.ird 

figures  in  all  12  update  ;  <-r  io  I  .  Th<  ■■  i  *  ul  ♦  ••  *». ;  it  :vd  the  ;ieo  ity 

of  a  dither  signal  and  valid)*  ml  t  i  :■  .  '  .  *  i  r  |twu. 

With  the  use  of  the  reference  1,-vl  ;  im  a  ur'-m.’.i  n  a  si  given  in 

Table  the  WLS  parameter  estimator  will  n /  t  r<  due.-  ae.-eptable  results 
if  the  data  window  lo.ding  p.  t  1  od  is  i  i .  ■  t «  ...e  1  ana  t  ;,e  dither  signal  is 
used.  The  effectiveness  is  also  improved  I  y  istim:  ing  tin  ]engt  ti  of  t  tie 
dither  period.  Table  11  summarizes  t  !..  results  of  tin  first  update  cycle 
of  simulations  of  the  Tokuo  M.it  sailing  >v«  •  i  the  st  ej  peg  U,t  tom  shown  in 

Fig.  25  under  ttie  control  of  an  adaptive  path  controller  using  the  WLS 


parameter  estimator  when  subjected  to  the  reference  level  of  noise.  Results 
are  shown  for  various  data  window  lengths  and  dither  signal  periods.  The 
initial  value  of  each  parameter  is  only  used  in  the  first  cycle;  in  subse¬ 
quent  cycles  the  most  recent  valid  estimate  is  used  as  the  starting  point 
for  the  DUD  algorithm  search.  With  a  data  window  length  and  dither  period 
of  0.5  ship  length,  the  results  failed  the  validity  check  and  the  algorithm 
failed  to  converge  at  all  after  the  depth  change  at  t'=3.  Results  are  also 
shown  for  a  data  window  length  of  5.0  ship  lengths  (1000  samples)  and  dither 
periods  of  0.5,  1.0,  and  5.0  ship  lengths.  In  these  cases,  only  one  update 
period  of  6.0  ship  lengths  was  simulated.  In  general,  the  parameter  esti¬ 
mates  are  reasonably  good  and  improve  with  the  dither  signal  period. 


parameter 

initial 

value 

H/T=l . 89 

correct 

value 

H/T=» 

WLS  parameter  es 
window  length  (sh 
dither  period (shi 

timates  using 
ip  lengths)/ 
o  lengths) 

0.5/0. 5 

5. 0/0. 5 

5. 0/1.0 

r 

5 . 0/2 . 5 

f22 

-1.7657 

-1.9515 

1.2835* 

-1.8179 

-1.8323 

-1.8582 

f23 

5.7359 

3.1591 

-9. 8155* 

2.8715 

2.8902 

2.0648 

'2, 

-0.88074 

-1.0410 

-0.93100 

-1.0146 

-1.0229 

-1.0216 

f  32 

0.17199 

0. 31507 

0.95165* 

0.36378 

0. 32584 

0.31693 

f  33 

-0.52766 

-0.63651 

-3.1628* 

-0.72829 

-0.64829 

-0.63519 

f  35 

-0.15607 

-0.16163 

-0.13573 

-0.15646 

-0. 16053 

-0.16101 

offset  n' 

at  end  of  data 

window 

-0.0005 

-0.2854 

-0. 5706 

-1.435 

number  of 

evaluations  of 

cost  J 

68 

60 

107 

66 

*  estimate  outside  permissible  range  in  Table  10 
Table  11.  Effect  of  Data  Window  Length  and  Dither  Period  on  WLS  Parameter 


Estimates  Without  Disturbances. 

The  results  shown  in  Table  11  for  t  tie  5  ship  length  data  window  arid 
0.5  ship  length  dither  signal  probably  represent  about  the  best  estimate:, 
that  could  be  achieved  with  a  WLS  parameter  estimation  in  a  ship  path  con¬ 
troller.  The  offset  due  to  the  dither  signal  at  the  end  of  the  data  window 
is  about  83m  or  1.74  times  the  ship  beam.  Shorter  window  lengths  and  shorter 
dither  periods  would  reduce  this  value  but  produce  less  accurate  parameter 
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estimates.  Simulations  results  with  a  2  . 5  shi;  length  data  window 
are  shown  in  the  next  section.  Since  we  showed  in  S.-dion  '.4  that  th.e  con¬ 
trol  loop  only  really  needs  accurate  estimates  of  f  and  f  ,  -the  accuracy 

2  ?  2  5 

of  the  estimates  of  only  these  two  parameters  should  be  the  basis-  of  judqing 
the  effectiveness  of  a  WI,S  parameter  estimator  designed  for  use  in  an  adaptive 
path  controller.  A  design  tradeoff  would  be  necessary  to  achieve  the  l»st 
overall  performance  of  the  path  controller.  The  longer  dither  periods  could 
not  be  used  in  a  path  controller  but  are  of  interest  if  a  WLS  parameter 
estimator  were  to  be  used  for  system  identification  in  unrestricted  waters. 

It  is  interesting  to  note  in  Table  11  that  the  short  dither  period 
allows  good  estimates  of  the  parameters  in  the  r-equation  (f  ^  f i ' anf^ 
but  relatively  poorer  estimates  of  the  parameters  in  the  h-equat ion 
(f  ,  f^,and  The  longer  dither  periods,  however,  significantly  im¬ 

prove  the  estimator's  ability  to  estimate  the  parameters  in  the  .-equation. 
This  can  be  explained  by  consideration  of  t  lie  system  time  constants  and 
eigenvectors.  The  open-loop  eigenvalues  and  eigenvectors  of  the  5-state 
model  of  the  Tokyo  Maru  include  two  zero  eigenvalues  and  a  fluid  associated 
with  the  rudder  time  constant.  The’  other  two  eigenvalue.-;  and  the  associated 
eigenvectors  for  deep  water  are  as  follows: 

=  -0.0992,  i,  =  10.08  ship  lengths?, 

4  4 

f  =  [-0.0933,  0.0093,  0.0054,  O.'h'f.H,  0.0]'1, 

1  -  -2.489,  i  -  0.402  ship  length, 

5  5 

■'  =  [-0.3669,  0.91  32,  -0.1553,  -0.0850,  0]'. 

-r> 

Eigenvector  i  has  a  time  constant  of  '’.4  ship  length  and  is  dominated  by 
the  second  component  (r  ').  A  dither  r  i  <  •<  1  of  0.5  ship  length  is?  effective 
in  exciting  this  mode  and  thus  allows  effect  i  ve  imat  ion  of  the  para¬ 
meters  in  the  r-cquation.  Eigenvector  '  ^  i  ••  dominat'd  by  the  fourth  com¬ 
ponent  (’•').  Although  ;  '  is  only  a  weak  component  of  e  i  <n  nvect  <>t  ' 


it  is  directly  fed  into  the  n-equation  and  n  is  the  main  component  of 
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With  a  time  constant  of  10  ship  lengths,  this  mode  is  more  effectively 
excited  by  the  longer  dither  periods  thus  making  the  parameters  in  the 
i’-equation  more  identifiable. 

4.2.2  Performance  with  Bias  Disturbances 

In  a  ship  path  controller  application,  the  WLS  parameter  estimator 
would  also  have  to  estimate  the  yawing  moment  and  lateral  force  distur¬ 
bances  acting  on  the  ship  as  included  in  eq.  (74).  To  evaluate  the  per¬ 
formance  of  the  WLS  parameter  estimator  with  time-invariant  disturbances, 
we  simulated  the  Tokyo  Maru  under  the  control  of  an  adaptive  ship  path 
controller  using  this  algorithm  while  being  subjected  to  a  constant 
1-knot  lateral  current  as  defined  in  Fig.  4  for  t'<15.  This  disturbance 
was  held  constant  for  18  ship  lengths;  the  reference  measurement  noise 
level  from  Table  9  was  used.  The  controller  and  parameter  estimator  were 
intialized  with  ship  parameters  for  H/T=1.89  but  the  ship  was  simulated 
to  be  actually  operating  in  deep  water  (H/T="») .  The  initial  conditions 
for  these  simulations  were  the  steady-state  conditions  for  the  control 
loop;  i.e.  , 

<>’  =  =  e,‘  =  S'  =  -0. 1428, 

r'  =  r’  =  n'  =  n’  =  0.0, 

6’  =  6'  =  0.1310, 

N'  =  0.001992  , 

Y '  =  0.003437  . 

The  disturbance  shown  in  Fig.  4  is  not  exactly  equivalent  to  a  lateral 
current  when  the  ship  is  at  any  depth-to-draft  ratio  other  than  H/T=1.89 
for  which  it  was  derived  so  the  steady-state  rudder  angle  is  not  zero. 

The  constant  disturbance  simulations  used  a  parameter  estimation  cycle 
of  6  ship  lengths  (LU=1200  samples)  and  three  combinations  of  data  window 
length  and  dither  period.  The  results  of  the  first  parameter  estimation 
cycle  of  these  simulations  are  shown  in  Table  12.  The  results  with  window 


length  of  2.5  ship  lengths  (500  samples)  and  a  0.5  ship  length  dither  period 
are  excellent.  Results  with  a  5  ship  length  data  window  and  two  dither  periods 
are  shown  as  a  reference  point  for  the  effectiveness  of  the  WLS  parameter  es¬ 
timator.  These  designs  would  be  impractical  in  a  path  control  application 
due  to  the  excessive  offset  from  the  desired  path.  These  paths  would  also 
clearly  violate  the  validity  of  the  linear  ship  model.  The  effects  of  the 
assumed  constant  values  for  the  y . .  coefficients  in  eq .  (75)  on  the  estimates 
of  the  disturbances  N'  and  Y'  can  be  clearly  seen  in  Table  12.  In  addition 
to  allowing  a  more  reasonable  offset  from  the  path,  the  reduction  of  the 
window  length  from  5  to  2.5  ship  lengths  reduces  the  computation  cost  of 
each  evaluation  of  the  cost  J  by  about  one-half  so  the  shorter  window  length 
results  in  about  half  the  data  storage  and  computational  load  with  only  a  small 
loss  in  parameter  estimate  accuracy. 

]  initial  ~  correct  Twls  paiameter  estimates  using  window! 
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Table  12.  Effect  of  Oata  Window  length  and  Dither  Period  on  WLS  Parameter 
Estimates  wit):  Constant  Disturbance: 

;v>  further  illustrate  the  performance  of  t  he  WI.S  parameter  estimator 
wit),  1  2.5  ship  length  data  window  and  a  0.5  ship  length  dither  period  when 
subjected  to  a  constant  disturbance,  the  ;ifii  j  trajectory  for  t  lie  simulation 
shown  in  Table  12  is  shown  in  Pig.  70.  Pig.  26  •;):■, ws  t  he  lateral  offset 
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response  of  the  ship  for  the  first  three  update  periods  or  18  ship  lengths 
The  dither  signal  has  a  zero  mean  rather  than  a  mean  value  corresponding 
to  the  equilibrium  rudder  angle  for  the  constant  disturbance  when  the 
controller  is  using  the  initial  parameters.  As  a  result,  the  first  dither 
period  gives  the  ship  a  strong  motion  to  the  right  ( n ' > 0 ) .  This  produces 
improved  estimates  in  Table  12  compared  with  the  corresponding  results 
in  Table  11.  As  soon  as  the  controller  is  updated  with  nearly  correct 
parameters  at  t'=6.o,  the  equilibrium  rudder  angle  is  essentially  zero 
and  subsequent  dither  periods  have  less  of  an  effect  on  the  ship's 
offset  from  the  desired  track.  The  maximum  offset  is  almost  one-half 
ship  length  during  the  second  update  cycle.  The  simulation  represents 
a  type  of  "start-up  transient"  for  the  adaptive  control  system  so  the 
lateral  offset  in  the  third  update  period  more  closely  represents  the 
offset  which  the  dither  signal  would  cause  in  subsequent  update  cycles. 


Figure  26.  Lateral  Offset  Response  to  Constant  Disturbance  —  WLS  Para¬ 
meter  Estimator. 

The  parameter  estimate  results  for  the  WLS  parameter  estimator  simula¬ 
tion  shown  in  Fig.  26  are  shown  in  Table  13.  The  estimates  were  acceptable 
in  each  update.  The  Kalman  filter  ship  equations  and  the  gains  C  and  K 
were  therefore  updated  at  the  end  of  each  eycl".  Since  the  initial  parameter 
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estimate  for  the  second  and  third  estimation  cycles  was  the  estimate 
produced  by  the  previous  estimation  cycle,  the  number  of  evaluations 
of  the  cost  function  decreased  significantly  after  the  first  update. 

The  accuracy  of  the  estimates  of  the  6-equation  parameters  appears  to 
decrease  as  the  ship  motion  is  reduced.  In  general,  the  parameter  esti¬ 
mates  are  very  good.  Recall  that  the  disturbance  estimates  are  never 
used.  A  shorter  data  window  length  might  continue  to  produce  parameter 
estimates  of  acceptable  accuracy  with  a  smaller  offset  from  the  desired 
path.  Further  design  tradeoff  studies  could  establish  the  best  compromise. 
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19 

final  cost  JxlO^ 
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update  time 
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Table  13.  Parameter  Estimates  from  Constant  Disturbance  emulation 


4. 3  Time-Varying  Parameters  Performance 

In  this  section,  we  evaluate  the  performance  .  l  a  w  :■  nip  j>ar.»me*  ».-r 
estimator  when  the  parameters  are  time-varying.  In  1  he  i  icvu'Uh  :  e*-t  j,>n  , 
we  considered  the  performance  when  the  parameters  wen  '-oi.  .tint  during  •  vl. 
data  window  loading  period  which  might  have  to  lie  is  h.ng  as  .  •  .hit  lengths. 


The  parameters  could  however,  change  from  one  parameter  estimation  cycle  to 
the  next.  Here  we  consider  situations  in  which  the  parameters  change  contin¬ 
uously  during  the  data  window  loading  periods.  We  continue  to  use  the  water 
depth  as  the  independent  variable  which  defines  the  changing  parameters.  We 
consider  two  cases.  First,  we  consider  the  situation  in  which  the  water  depth 
is  changing  continuously  so  that  the  ship  coefficients  are  changing  in  a  con¬ 
tinuous  manner.  We  then  consider  the  situation  in  which  the  disturbances  in 
the  parameter  vector,  eq.  (74),  change  continuously  during  the  data  window 
loading  period  as  they  would  when  the  passes  another  ship  or  fixed  obstruction. 

4.3.1  Time-Varying  Ship  Coefficients 

To  evaluate  the  performance  of  the  WLS  parameter  estimator  with  continu¬ 
ously  varying  ship  coefficients,  we  simulated  the  Tokyo  Maru  under  the  control 
of  an  adaptive  path  controller  using  the  WLS  parameter  estimator  while  sailing 
over  a  downward  sloping  bottom.  This  simulation  was  performed  without  any 
measurement  noise  using  a  very  short  data  window  loading  period  of  only  0.2 
ship  length.  The  update  cycle  was  0.3  ship  length.  The  bottom  profile  was 
constructed  to  be  level  at  H/T=1.89  during  the  first  data  window.  The  esti¬ 
mator  was  initialized  with  parameter  values  at  H/T=1.89  and  this  startup  period 
was  introduced  to  allow  the  system  to  stabilize  prior  to  experiencing  a  change 
in  parameters.  From  time  t*  =  0.35,  which  was  the  beginning  of  the  second  up¬ 
date  period,  the  bottom  began  a  downward  slope  in  which  the  depth-to-draf t  ratio 
increased  linearly  at  a  rate  of  0.105/update  or  0.350/ship  length.  The  results 
of  the  first  10  updates  in  this  simulation  are  summarized  in  Table  14.  In 
many  updates,  the  estimates  were  far  from  the  true  values.  In  many  instances, 
the  algorithm  produced  parameter  values  which  were  beyond  the  permissible  ranges 
defined  in  Table  10.  The  estimates  began  to  improve  somewhat  after  the  initial 
transient  caused  by  the  abrupt  change  in  bottom  slope  at  t'=0.35.  Net  ice  that 
the  estimates  do  not  approximate  the  mean  value  of  each  parameter  during  the 
data  window  loading  period  as  would  be  desired.  Considering  that  the  simula¬ 
tion  summarized  in  Table  14  is  without  measurement  noise  and  used  only  a  0.2 
ship  length  data  window,  we  conclude  that  the  WLS  parameter  estimator  can  not 
effectively  estimate  time-varying  ship  coefficients.  With  realistic  measure¬ 
ment  noise,  the  data  loading  window  period  would  have  to  be  at  least  two  ship 
lengths  as  shown  in  Section  4.2.  The  bottom  slope  used  m  the  simulation  sum¬ 
marized  in  Table  14;  i.e.  about  one  draft  in  three  ship  lengths,  is  not  un¬ 
realistic!'/  high. 
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The  failure  of  the  WLS  parameter  estimator  to  effectively  estimate  time- 
varying  ship  coefficients  can  be  explained  by  considering  the  formulation  of 
the  algorithm.  The  true  ship  response  can  be  considered  as  the  response  of  a 
system  of  time-varying  linear  differential  equations.  The  coefficients  in 
these  equations  vary  continuously  during  the  time  period  of  the  data  window. 

The  ship  model,  on  the  other  hand,  is  the  time-variant  system,  eq.  (75),  and 
we  seek  the  parameter  vector  g*  which  minimizes  eq.  (77).  We  expect  this  para¬ 
meter  estimate  to  be  within  the  permissible  ranges  defined  in  Table  10  and  fur¬ 
ther  would  hope  that  the  estimates  are  the  mean  values  of  each  parameter  during 
the  data  window  loading  period.  There  is  no  guarantee  at  all  that  the  constant 
parameter  £*  which  minimizes  eq.  (77)  will  meet  these  objectives.  The  results 
shown  in  Table  14  show  that  even  without  measurement  noise  and  using  a  unrealis¬ 
tically  short  data  loading  window  period  the  estimates  do  not  meet  these  objec¬ 
tives. 

4.3.2.  Time-Varying  Disturbances 

Even  when  the  ship  coefficients  are  essentially  constant  during  the  data 
window  loading  period,  it  is  possible  that  the  ship  will  be  subjected  to  exter¬ 
nal  disturbances  which  vary  over  this  period.  To  evaluate  the  performance  of 
the  WLS  parameter  estimator  in  this  type  of  situation,  we  simulated  the  Tokyo 
Maru  under  the  control  of  an  adaptive  path  controller  using  the  WLS  parameter 
estimator  while  sailing  past  another  ship.  The  bottom  was  assumed  to  be  constant 
deep  water  (H/T=°°) .  The  ship  was  given  zero  initial  conditions  x=x=0;  the  con¬ 
troller  and  parameter  estimator  were  initialized  for  parameters  at  H/T=1.89. 

The  reference  measurement  noise  from  Table  9  was  utilized.  The  ship  was  simu¬ 
lated  for  15  ship  lengths.  It  was  simulated  to  pass  another  ship  as  defined 

I 

by  Fig.  3  beam-to-beam  at  t=7.  The  WLS  parameter  estimation  cycle  was  5  ship 
lengths.  The  data  window  length  was  2.5  ship  lengths;  the  dither  period  was 
0.5  ship  length.  With  an  update  period  of  5  ship  lengths,  three  update  cycles 
were  completed  during  the  simulation.  The  effects  of  the  passing  ship  were 
felt  from  t'=5  through  t'=8.5  so  the  disturbance  was  present  during  tie  first 
part  of  the  window  loading  period  for  the  second  update. 


A  summary  of  the  results  of  the  three  update  cycles  of  the  simulation  is 
presented  in  Table  15.  During  the  first  update,  the  parameter  estimates,  pai- 
ticularly  f  and  are  needed  by  the  control  loop,  are  close  to  the 


true  values.  The  disturbance  estimates  N'  and  V  are  within  the  noise  levels 
corresponding  to  these  quantities.  At  the  end  of  this  update  cycle,  the  para¬ 
meters  in  the  Kalman  filter  ship  model  won  updated.  'Jew  gain  matrices  C  and 
K  were  calculated  and  updated  in  the  system.  !>ai  inn  the  second  update  period 
the  passing  ship  disturbance  starts  to  act  at  the  ship  is  running 

under  the  command  of  the  updated  controller.  At  •'  -  ,  tne  disturbance.'  reach 

their  maximum  values  of  about: 

N'xlO5  =  -20.4, 

Y'xlO5  =  48.0. 

The  second  data  window  was  loaded  from  t'~-7.4:5  to  t*  -9.  >•'.  : Airing  this 

period,  the  disturbances  declined  r-'om  their  peak  values  to  ::ero  at  t'=8.50. 

The  DUD  algorithm  search  was  very  slow  to  converge  in  the  second  cycle  and  the 

results  were  bad.  Three  parameters  were  outside  tin-  permissible  ranges  defined 

in  Table  10.  Since  the  estimate  f  _  was  beyond  th<  Permissible  range,  there 

2  j 

was  no  update.  During  tint  third  updat.  period,  the  disturbance  was  gone. 

The  parameter  estimates  were  good  and  the  controller  was  updated  a  second  time. 
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Vie  conclude  from  this  simulation  that  the  WLS  parameter  estimator  cannot 
estimate  the  parameters  when  there  is  a  significant  time-varying  distur¬ 
bance  present  during  the  window  loading  period.  During  the  second  update 
sho'-  v  n  Table  15  the  algorithm  tried  to  approximate  the  time-varying 
disturbances  by  constant  quantities  N*  and  Y'.  In  the  process,  it  adjusted 
the  other  parameter  estimates  to  minimize  the  cost  J.  Since  the  relative 
effect  of  tr  ’ sturbances  on  J  is  much  larger  than  the  changes  due  to  the 
system  parameters,  the  final  estimates  of  the  system  parameters  were  sub¬ 
stantially  off.  The  system  protected  itself,  however,  through  the  check  on 
the  validity  of  the  parameters.  When  the  parameter  estimates  in  the  second 
update  failed  the  validity  check,  no  update  was  made  and  the  controller  con¬ 
tinued  to  operate  with  the  gains  obtained  in  the  first  update.  The  slow 
rate  of  convergence  in  the  second  update  is  also  an  indication  that  time- 
varying  disturbances  were  present  during  the  data  window  loading  period. 

Thus,  the  estimation  process  could  be  stopped  without  an  update  after  a 
prescribed  number  of  evaluations  of  the  cost  function  as  further  protection 
against  invalid  results. 

4 . 4  Computational  Requirements 

The  WLS  parameter  estimator  would  place  a  fairly  large  dynamic  data 
storage  and  computational  load  on  an  onboard  computer.  The  parameter  esti¬ 
mation  and  gain  update  calculations  could,  however,  be  performed  on  a  batch 
basis  as  time  is  available.  In  developing  the  simulation  program  used  in 
our  work,  we  made  no  serious  effort  to  conserve  storage  locations  or  CPU 
time.  The  program  was  therefore  very  expensive  to  run.  The  Euler  integra¬ 
tion  step-size  used  in  the  various  simulations  presented  above  was  taken  as 
0.005  ship  lengths  or  about  .24s.  For  convenience,  this  was  also  used  as 
the  sample  time.  In  a  practical  application  of  a  WLS  parameter  estimator 
on  a  ship  such  as  the  Tokyo  Maru,  the  sample  time  could  easily  be  extended 
to  perhaps  Is.  With  a  data  window  of  2.5  ship  lengths  at  the  12  knots  used 
in  our  simulations,  the  data  window  would  require  120  locations  for  each  of 
the  four  measurements  and  the  rudder  angle  plus  an  additional  4  locations 
for  the  initial  state.  The  DUD  algorithm  as  used  here  would  require  another 
4057  stoiage  locations.  This  would  give  a  total  number  of  dynamic  storage 
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locations  of  4661  for  the  implementation  of  the  WLS  parameter  estimator  ] 

A  WLS  ship  parameter  estimator  would  therefore  be  feasible  but  it  would  ] 

place  fairly  heavy  storage  and  CPU  time  demands  on  an  onboard  computer.  If  ,  3 

dynamic  storage  were  a  limiting  factor,  alternative  search  algorithms  could  \ 

be  considered.  For  example,  the  Helder  and  Mead  STMI'LKX  algoritlun  was  used 
in  some  of  our  earlier  work.  This  algorithm  generally  required  more  iterations 
than  DUD  to  converge  but  would  require  only  81  additional  dynamic  storage 
locations.  This  would  give  a  total  number  of  dynamic  storage  locations  of 
only  685.  The  reduced  storage  would  be  offset  by  a  greater  CPU  time  require¬ 
ment.  :  i 
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5.  Minimum  Variance  Parameter  Estimator 


In  this  section,  we  present  the  derivation  and  evaluation  of  a  Mini¬ 
mum  Variance  Parameter  Estimator  (MVE)  which  could  be  used  in  the  gain  up¬ 
date  loop  of  an  adaptive  ship  path  control  system.  This  section  is  based 
upon  Cuong's  Ph.D.  Dissertation.  The  MVE  is  a  probabilistic  digital  fil¬ 
tering  technique  which  can  be  used  to  determine  a  set  of  parameters  of  the 
ship  equations  of  motion.  The  filter  gains  are  chosen  to  minimize  the  trace 
of  the  estimate  error  covariance  matrix.  It  is  a  recursive  filtering  tech¬ 
nique  where  the  new  parameter  estimate  is  formed  from  the  previous  estimate 
plus  a  proportional  constant  times  the  innovation.  The  filter  gains  and  pro¬ 
pagation  of  the  error  covariance  can  be  calculated  recursively  based  on  pre¬ 
vious  values  and  recent  measurements. 

The  derivation  of  this  parameter  estimator  parallels  the  derivation  of 
the  state  estimation  Kalman  filter.  It  is  more  complicated,  however,  because 
the  measurements  are  contaminated  with  both  additive  and  multiplicative  noise. 

Vi 

This  method  was  first  proposed  by  Kotob  and  Kaufman  for  a  linear  system 
with  measurement  noise,  time-invariant  parameters,  and  no  process  disturbances. 
They  used  this  filter  to  estimate  the  system  parameters  of  the  F-8  DFBW  air¬ 
craft-  They  showed  it  to  be  more  efficient  and  effective  than  an  Extended 
Kalman  Filter  and  a  Weighted  Least-Squares  Parameter  estimator  in  three  simu¬ 
lated  performance  tests.  They  also  proposed  an  extension  for  time-varying 
parameters.  Here  we  extend  the  MVE  approach  to  accommodate  zero-mean  process 
disturbances  and  time-varying  parameters. 

5. 1  Derivation  and  Development 

The  Minimum  Variance  Estimator  is  a  linear  parameter  identification 
scheme  which  has  three  inherent  restrictions.  First,  the  equations  of  motion 
must  be  linear  in  the  estimated  parameters  jd.  Second,  in  order  to  separate 
the  state  estimation  problem  from  the  parameter  identification  problem,  it 
is  required  that  all  states  be  available  for  measurement.  Finally,  the  al¬ 
gorithm  as  developed  here  can  accommodate  only  zero-mean  process  disturbances. 

The  ship  path  control  problem  must  be  formulated  to  satisfy  the  restric¬ 
tions  of  the  MVE.  First,  the  linearity  requires  that  there  be  no  coupling 
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among  the  estimated  parameters.  This  condition  is.  automat  ionl  ly  .-..it  i  at  io«l 
if  we  choose  the  parameter  vector  to  be: 


li  = 


T 

If22'f23'fii5'f  32'f33'f3s' 


(78) 


The  separation  condition  demands  that  we  utilize  the  4-state  system  because 
we  can  measure  only  £',r'  and  n' .  Tlie  external  yawing  moment  N '  and 

side  force  Y1  are  not  available  through  direct  measurement.  Any  attempt 
to  use  their  estimates  from  the  Kalman  Filter  in  parameter  estimation  will 
represent  a  coupling  between  state  and  parameter  estimators  which  we  want 
to  avoid.  Fortunately,  the  precise  knowledge  of  the  disturbances  is  not 
necessary  in  this  case  since  we  only  need  the  covariances  ot  the  disturbances 
in  the  calculation  of  estimator  gain. 


As  mentioned  earlier,  the  MVE  is  a  recursive  filtering  technique  which 
makes  it  well  suited  to  real-time  operation.  A  possible  timing  sequence  for 
algoritlim  is  shown  in  Fig.  27.  The  parameters  are  not  identified  at 
every  measurement  sample  time.  In  this  case,  the  parameters  are  estimated 
at  intervals  of  f  =10.  (The  derivation  of  this  algorithm  requires  that 
{’  ■  1 . )  A  preliminary  study  was  carried  out  to  compare  performance  of  an  esti¬ 
mator  with  with  that  of  an  estimator  with  f=10.  Although  the  higher 
identification  frequency  estimator  increased  the  computational  load  five 
times,  it  did  not  present  any  clear-out  superiority  compared  with  the  lower 
frequency  estimator.  For  {'=10  and  the  sample  time  assumed  in  this  work, 

estimator  would  have  about  2.4  seconds  for  each  update  each  cycle .  This 
represents  a  realistic  computational,  time  allotment.  Being  a  recursive 
technique,  the  MVE  has  a  very  small  dynamic  data  storage  t oqu i rcmwit. 

At  each  recursive  step,  the  algorithm  needs  two  consecutive  measurements 

vector;--  ,  and  z'  ,  the  recent  rudder  angle  ,,  the  last  parameter  esti 
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Figure  27. 


Timing  Sequence  for  Minimum  Variance  p,iramet<  i 


le  t  imat  >i 


formation  requires  only  34  storage  locations.  The  entire  algorithm  could 
be  implemented  with  less  than  100  storage  locations. 


r 


The  Minimum  Variance  parameter  estimation  algorithm  is  derived  for  a 
discrete  time  system  of  equations.  The  continuous  ship  system  equations 
from  Section  4.1;  i.e., 

x  =  Fx  4  G6*  +  !'w  ,  (72) 

must  therefore  be  transformed  into  a  system  of  difference  equations.  These 
differential  equations  can  be  represented  by, 


Fx  4  G6 '  4  Tw  , 

—  k  k  —  k 


or, 


xk  +  1  =  [T  4  AtF]xk  4  AtG6 1  k  4  Atrw  . 


(79) 


If  we  let  F'  =  I  4  AtF,  G'  =  AtG,  and  C  =  Atr,  the  dynamics  of  the  ship 
can  then  be  represented  by  the  difference  equations, 


or, 


x  ,  i 

-  k4l 


=  F'x  ,  4  G'6V  4  T'w 
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The  disturbances  are  now  represented  by  two  discrete  stochastic  sequences 

f  N ' k }  and  {V  }  which  we  assume  to  be  independent.  The  y.  .  coefficients 

are  given  constant,  assumed  values  y.  .  as  in  Section  4.  All  four  states 
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are  assumed  to  be  available  for  measurement;  i.e.. 


z 

-  k 


^k+  **k' 


(82) 


where  the  noise  (v1  }  is  a  Gaussian  white  sequence  with  covariance  matrix  K. 

Fill  owing  an  approach  similar  to  that  used  in  Section  3.3,  the  parameters 
can  re  modeled  as  Brownian  motion  processes.  In  this  discrete  formulation, 
the  parameter  vector,  eq.  (78),  is  therefore  modeled  as, 
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£k  =  ^k-£  +  m) 

whore  '  «  5  is  a  Gaussian  white  sequence  with  diffusion  coefficient  matrix  3. 

The  parameter  i  defines  the  frequency  of  idontif ication .  It  is  required  that 
i!  >  I  throughout  this  work.  The  diffusion  coefficient:;  can  he  selected  with 
east  using  the  properties  of  the  Brownian  motion  process  as  developed  in 
Section  3. 3. 

In  order  to  estimate  the  unknown  parameters  in  F1  and  G1 ,  it  is  necessary 
to  rearrange  eq.  (81)  so  that  expressed  explicitly  in  terms  of 

and  ,  tne  unknown  and  known  parameter  vectors,  respectively.  In  this  par¬ 
ticular  problem,  the  known  parameter  vector  is  simply  =[l]  .  The  revised 
.5/3 ten  equations  become, 
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whet*  x,  -i>:  ,x  „ ,  X  ,  x,  .]  .  Note  that  c,  arid  0,  are  matrix  function: 
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of  x,  .  Substituting  x,  -v*  ,  into  C,  and  H,  ,  these  matrices  can  he  do- 
—  x  —  k  k  —  k  k  k 

composed  into  deterministic  and  stochastic  component.:;  which  depend  on  ?. '  an<! 
V  :j  .ai  a'  e  ly  ;  i  .  e  .  , 


•(xk>  -  C(z'k)  -  C(v'k)  *  Ck  -  Ck  . 


where  we  def  \  , 
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vation  which  is  corrupted  by  the  multiplicative  noise  tern  r  c  and  t  be 

V 

additive  noise  terms  v'  and  Us,.  Shifting  the  index  )  v  fine  step  buck  - 

—  k  +  1  k- k 

ward,  eq .  ('ll)  becomes, 

~  A.  V  V 

,  R  r  r,  ,  -  r  u  „  -  b  s  tv’  t-  !'*w  , 

-k  k- 1  *-  k  - 1  k- 1  *-k-{  k-1  k-1  k  k-1 

where  we  have  used  p,  =  p,  „  because  eq.  (B3)  approximates  the  stochas- 

*-  k-  L  k-t 

tic  process  (£»(t)}  by  the  stochastic  step  function  p  I 
We  can  now  define  the  innovation  vector  at  time  k  as, 

ik  "t'k  ■  tn" 

where  the  second  term  on  the  right-hand  side  is  the  conditional  expect  ion 
of  z 1  with  respect  to  the  whole  measurement  history;  i  , 
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Noting  that, 


V.  k - 1  ^k-t!  +  Dk-l  -k-1 


-1  -k 


-*k  'k-11  '"k-1  £k-f  *  °k- 


1  "k-1  ' 


,(■  disturbances  ’  w  are  restricted  to  zero-mean  processes . 
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ibst  ;  toting  eq.  (02)  into  eq.  (04i  yields, 
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Using  an  approach  parallel  to  the  Kalman  filter  derivation,  we  can 
assume  that  the  parameter  estimate  can  be  given  by  the  following  scheme, 

Ek  =  Ek_^  +  \  Ik  - 

where  the  gain  matrix  is  to  be  established.  Subtracting  eq.  (83)  from 
eq.  (96),  we  get, 

=  £k-£  +  \  lk-n  k-t  ■  (-17) 

If  we  multiply  eq.  (97)  by  its  transpose,  we  get, 

~~T'V  ~T  TT  T  ~  TT  *T 

Ek  Ek  -  Ek_f  Ek_{*  \L*Lk  \  *  Ek_e  Vt  *  Ek -th\  *  k-f 

-  2k.elX  -  Shi-t  -  h-e^l-e  -  Ek-fCf  •  <*>> 

Defining  conditional  expections  as  follows: 

W  =  E[Yklzk-i]  - 
pk  =  E2t£k£Tk  1  < 


J  =  Ez  [Vk’  ’ 

the  conditional  expectation  of  both  sides  of  eq.  (98)  with  respect  to  th< 
whole  measurement  history  ^ yields , 


pk  =  Pk-£  +  KkEz[Ik-X-k]Kk  +  *  +  Ez  [£k-f-k  1Kk  +  Vz  ll-vh- 
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^k-£  ik  ]\  -  Vz  [ikiz]  "  Ez[^k-£ik-e]  ■  Kzl^k-al-(]  ■  (99> 


Cuong  has  evaluated  the  various  conditional  expectations  on  the  right- 
hand  side  of  eq.  (99)  for  the  case  when  (>1.  if  these  results  are  sub¬ 
stituted  into  eq.  (99) ,  we  arrive  at  the  following  expression  for  the 
propogation  of  the  estimate  error  covariance, 

Pk  =  Pk-Z+  *  -  ■Sc  Vl  \-<  '  \-l  ~Ch  <  +  VS-1  Vf  K-1  +  Req'k-X 


where  wo  define  an  equivalent  noise  cava;  ianee  matrix  P  as, 

oq 


R  u  ,  =  ,  +  B  4-  R  +  f’Q  V 

eq, k-1  k-1  -w 


and  where  R  is  the  measurement  noise  covariance  matrix  and, 


A  v  T  vf 

^  ~  Ez  [C^  ^  ^  ^  ^  '  we^9*ltec*  covariance  matrix, 


v  T  vt 

B  ~  "z  ''Jk-  L  -k-l-k-1  °k-l 


0  =  E  [w  w  ]  ,  disturbance  covariance  matrix, 

w  z  -k-l-k-1 

Elements  of  these  matrices  will  be  developed  further  below.  Minimiza¬ 
tion  of  the  trace  of  the  estimate  error  covariance  P  will  produce  a  minimum 
variance  estimate  of  the  parameters  £  ,  .  Stationary  conditions  for  the  mini¬ 
mization  of  the  trace  of  P^  are  obtained  by  setting  all  derivatives  of  eq. 
(100)  with  respect  to  the  elements  of  equal  to  zero;  i.e., 

/>.  /Vp  /\rp 

OK,  [C.  ,  P.  »  C,  .  +  R  ,  .  ]-  2P.  „  C  =  0  . 

k  k-1  k-(  k-1  eq,k-l  k-t  k-1 

Rearranging  this,  we  get  the  following  expression  for  the  filter  gain: 
a  at  ,  -1 

-  P,  ,  ?,  .  [C.  .  P.  t  C.  +  P.  ,1  .  (102) 

.<  k-<  k-1  k-1  k-(  k-1  eu,k-l 

'si:.;  rq.  (102)  to  simplify  oq.  (100),  we  arrive  at  the  following  form  of 
the  error  covariance  matrix, 

A 

>\  =  r,  .  *  :  -  k.  c.  .  p,  „  .  n°n 

k  k-t  k  k-1  k-( 

If  t  i.e  unknown  parameters  are  t  ime- invariant ,  the  parameter  model  ii  f  fusion 
.-Off  f  ic.i  ent  I  is  zero  in  this  expression. 

T)‘.e  B  matrix  is  defined  below  eq.  (101)  as  a  conditional  export  at  ion 

V  _  V 

involving  the-  product  D  which  in  tins  case  is  just  ^  given  by 

•  •  :.  (BP) .  Kq .  (8B)  can  be  used  in  the  definition  of  B  to  yiel  1, 


•  c.tr  (■;•  •  ; 


The  matrix  P  i  ; ,  therefore* ,  a  constant  matrix  which  depend:,  ent  i  rely  •  *n  t  h< 


measurement  noise  variances  a.  and  sampling  time  At.  The  disturbance  power 
T  1 

matrix  PQ  P  can  be  developed  next  using  the  definition  of  n  below  eg. 
w  w 

(101).  We  have  assumed  that  (N'^l  arid  {y1  )  are  independent  stochastic 
processes  with  zero  means.  We  can  assume  that  these  disturbances  can  he 
modeled  as  white  noise  sequences.  The  O  matrix  then  becomes, 


0  =  E  (w  w  ,  ,  ] 

^w  z  —  k-1  —  k-1 


-*L 

At 


Y 

At 


where  and  q^  are  the  spectral  densities  of  the  corresponding  continuous 
white  noise  disturbances  N '  and  Y' .  Using  this  result  the  disturbance 
power  matrix  becomes, 

0  0 


q  r,T=  At 
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°  (qNY21Y31  +  V22732) 
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(qN~2i~31  +  qY~22~32) 

(qNY31  +  V32  > 
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The  disturbance  power  matrix  is,  therefore,  also  a  constant  matrix. 


The  weighted  covariance  matrix  U  is  defined  below  eq.  (101)  as  a 

v 

conditional  expectation  involving  the  vector  £  which  is  given  by, 


Vi  ^k-e 


0 

At  (Vk-l,2  Pk-^,1  T  Vk-1,3  Pk-£, 2 

At  (Vk-l,2  Pk-f,4‘+  Vk-1,3  Pk-C,5) 

0 


The  presence  of  unknown  parameters  Pk  ^  ^  within  this  matrix  presents  com¬ 
putational  difficulties.  These  parameters  are  known  to  be  within  a 
known  finite  range  so  for  the  purposes  of  evaluating  ,  it  is  reasonable 

to  assume  that  the  parameters  are  Gaussian  random  variables  with  known  means 
and  variances;  i.e., 


f’  f  P  i  1  = 


-  2  2 
E[  (pi  “  Pi>  )  =  ai 


i  =  1 ,2, . . .6. 


(105) 
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The  mean  can  be  taken  as  the  mean  of  the  probable  ranqe  of  the  parameter 
and  x.  might  be  taken  as  half  of  the  probable  range  of  the  parameter. 

Assuming  that  the  various  components  of  the  measurement  noise  are  orthogonal, 
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,  therefore,  also  a  constai.L 

Uk_1=  U 

which  c.in  be 

calculated 

in  advance. 
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.-  entire  equivalent  me. 

.xt-ment  noise 

covariai.ce  can  therefore,  lx:  calculated  in  advance.  If  the  unknown  parameters 

are  time -invariant,  the  variances  i.  arc  zero  in  eg.  (IOC).  The  U  and  [< 

L  oq 

matrices  might  then  be  updated  periodically  using  the  latest  estimates  p.  in 

'  1 

lieu  of  tic-  assumed  mean  value--.  t>.  . 

'  i 

Kotob  and  Kaufman  nott  1  that  a  substantial  bias  was  observed  to  appear 
in  the  parameter  estimates  using  their  time-invariant,  parameters  MV!::  algo¬ 
rithm.  Tlvse  biases  were  caused  by  the  dependence  of  the-  additive  and  mul¬ 
tiplicative  noise  terms  in  eg.  (01).  They  developed  a  "bias  reduction" 
scheme  to  reduce  or  eliminate  this  bias.  We  will  follow  their  basic  approach 
here  to  develop  an  bias  error  correction  procedure  for  the  general  time- 
varying  parameters  case.  Post.mult  iplying  both  sides  of  eg.  (10.1)  by  the  ex¬ 
pression  inside  the  bracket,  it  becomes, 

z'  ^  T  /vr 

Y  i  ,  r  +  K  R  -  P  „  c 

k  k-1  k-1  k-1  k  eg,k-l  k-t  k-l. 

Rearranging  this  expression,  we  arrive  at : 
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to  give, 


/\T  -1 

result  into  eq.  (103)  and  then  postmultiply  by  C  R 

k  1  eq , k  1 
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(108) 


Using  eq.  (94)  and  eq.  (108)  in  eq.  (96)  yields. 
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In  order  to  reduce  the  bias  in  estimating  £,  an  extra  term  can  be  added  to 
the  aba' 
choose , 


the  above  expression.  This  term  should  contain  and  p^-$and  ^  so  we  can 
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where  the  gain  matrix  G  is  yet  to  be  determined.  Equation  (109)  then 

K“1 

becomes , 

aa  at-I^ot-i^a  a 

4  =ik-f  +  [Pk"*]  [Ck-lReq,k-l-k"  Ck-lReq,k-lCk-l£k-f  +  Gk-lEk-C 

(110) 

We  can  now  take  the  expectation  of  both  sides  of  eq.  (110)  and  impose  the 
additional  condition  that: 

A  A 

lim  E  [jt^J  <y  lim  E  [o^ ^  1 
k  — ™  k — 


This  assumes  that  after  initial  transients  die  out  and  the  estimator  is 
tracking  a  change  in  parameters,  the  parameter  change  between  the  parameter 
estimates  is  small.  If  the  parameters  are  time-invariant  this  condition 
is  exact.  With  time-varying  parameters,  the  effectiveness  of  this  assump¬ 
tion  must  be  evaluated  by  simulation.  Noting  that  [P  -  <*>]  is  always  nonzero, 

k 

the  expectation  of  quantities  inside  the  bracket  in  eq.  (110)  must  be  zero; 
i  .  e .  , 


E  [C 


T 

k-1 


R 


eq, k-1 


Ck-1  Req,k-1  Ck-1  E-k-l  +  Gk-1  ^k-£ 


0.  (Ill) 


Equation  (111)  can  now  be  used  to  establish  the  gain 


sion , 


A  v 

c  =  c  -  c 

k-1  k-1  k-1  ' 


The  o xpre s - 


can  be  used  in  eq.  (92)  and  this  can  be  substituted  into  eq.  (Ill)  .  p, 

A  V  A 

«»Xi  rosiion  C  -  C  +  C  can  lx;  used  to  eliminate  r  and  t  ? 
k-1  k-1  k-1  k-1 
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transpose  from  the  terms  involving  p^_^and  p^_^  •  Utilizing  the  fact  that 

ck-i 3,1(1  Vi-  \-i'and  Vi*  vk and  Vi-  Vi and  Vi* 30(1  Sc-i 3,1(1  K-i 

are  orthogonal  #  the  expectation  in  eq.  (Ill)  then  reduces  to, 

EI£k-<>  •  Elek-lRix-lS-l1’El^1  '  lE,Ck-lReJ,k-lCk-l'  *  E1*k-lE^,k-A-l 

-  Gk-1  K 


For  an  estimate  to  be  unbiased,  it  is  necessary  that, 

lim  E[£  J  =  lim  E[n  J 
k-«n  k-*°°  K 


which  will  be  true  if, 

Vt  -1  V 

3.  ,  =  E[C  ,R  .  ,C,  .] 
k-i  k-1  eq, k-1  k-1 


(112) 


This  term  is  a  generalized  measure  of  noise  and  disturbance  power. 


The  matrix 
expressions  for 
yields , 


'k-1 


= (At) ' 


G  can  be  developed  explicitly  for  this  problem  using  the 
vk_  -1 

C,  ,  and  R  ,  ,  =  [t . . ] .  The  expectation  term  by  term 
k-1  eq,k-l  13 


t22°2 

0 

0 

C32^2 

0 

0 


0 

t22°3 

0 

0 

t  ,a2 
32  3 

0 


0 

0 

0 

0 

0 


t  o2 
23  2 

0 

0 

‘33°  2 


0 

0 


t23°3 

0 

0 


,  2 
t  0 
33  3 


0 

0 

0 

0 

0 

0 


(113) 


The  matrix  G,  ,  is  thus  a  constant  matrix  which  can  be  calculated  in  advance 

k_L  -1  2 
knowing  the  elements  of  R  and  the  measurement  noise  variances  a.  . 

eq  1 

The  final  parameter  estimate  equation  is  obtained  by  adding  the  bias  error 
correction  term  to  eq.  (96)  giving, 


£k=  *k-e  +  h 


-k 


+  [P,  -  4>]  G. 
k 


k-1  £k-£  * 


(114) 


This  completes  the  derivation  of  the  Minimum  Variance  parameter  estimator. 
The  complete  algorithm  is  summarized  below. 
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5.2  Algorithm  Summary 

Since  the  derivation  of  the  Minimum  Variance  parameter  estimator  has 
been  long  and  complicated,  it  is  useful  at  this  point  to  present  am  algo¬ 
rithmic  summary  of  a  typical  recursive  cycle.  The  equivalent  noise  spectral 

density  P  can  be  calculated  and  stored  in  advance  using  eqs.  (101),  (104), 
eq 

(105),  and  (106).  If  the  unknown  parameters  are  time-invariant ,  the  variances 
are  zero  in  eq.  (106).  Likewise,  the  generalized  noise  power  matrix 
can  be  calculated  and  stored  in  advance  using  eq.  (113).  Let  k  be  the  value 
of  the  time  index  at  the  next  parameter  estimate.  The  parameter  estimation 
cycle  is  provided  with  the  previous  parameter  estimate  ^  and  the  current 
estimate  error  covariance  matrix  ^  which  were  stored  from  the  previous 
cycle.  The  recursive  estimation  cycle  proceeds  as  follows: 

1.  At  t'=k-l,  a  set  of  measurements  z_'  ^  and  the  rudder  angle 

1  ar°  sampled; 

A  A 

2.  These  data  are  used  to  calculate  and  ^  using  eqs.  (85)  and 

(87) ,  respectively; 


3.  At  t=k,  a  second  set  of  measurements  is  sampled 

4.  The  pseudomeasurement  vector  is  calculated  from  its  definition; 


i.e.  , 


z .  =  z' 

—  k  — k 


Dk-1  ^k-1  : 


5.  The  filter  gain  matrix  is  calculated  using  eq.  (102);  i.e.. 


A  A«p 

K.  =  P.  .  C.  .  (C.  .  P.  ,  C.  .  +  R  ) 
k  k-t  k-1  k-1  k-t  k-1  eq 


(102) 


6.  The  updated  error  covariance  matrix  is  calculated  using  eq.  (103) ; 

i.e. , 

Pk  =  Pk-C  +  *  “  *k  Vi  Pk-£  ' 

where  diffusion  coefficient  matrix  4  is  zero  in  the  time -invariant  parameters 


7.  The  innovation  vector  is  calculated  using  eq.  (94);  i.e., 


^  A 

2  .  “  C, 


5-k“  ck-i  Zy-l 


8.  The  new  parameter  estimate  p^  is  calculated  using  eq.  (114);  i.e., 


2k  “  ^k-£  +  \h  +  (V  Gk-1  Zk-l 


(114) 


and  again  $  is  zero  in  the  time-invariant  parameters  case. 


This  completes  one  parameter  estimation  cycle  of  the  Minimum  Variance 
Parameter  estimator.  Compared  with  the  the  Weighted  Least-Squares  method, 
it  is  very  elegant  and  powerful.  It  is  potentially  capable  of  estimating 
time-varying  parameters  in  real  time.  Recalling  that  the  computational 
requirement  for  the  WLS  method  was  quite  large,  the  load  imposed  by  the 
MVE  would  be  small  by  comparison.  As  developed  here,  the  MVE  cannot, 
however,  accommodate  bias  disturbances  which  can  be  handled  by  the  WLS 
approach.  We  evaluate  the  performance  by  the  MVE  in  the  next  section. 


5. 3.  Time- Invariant  parameters  Performance 

In  this  section,  we  evaluate  the  performance  of  the  MVE  ship  parameter 
estimator  when  the  parameters  are  time-invariant.  Simulations  were  conducted 
with  and  without  zero  mean  disturbances.  From  Section  4.2  we  would  expect 
that  for  the  MVE  to  be  effective  with  disturbances,  the  system  must  be  persis¬ 
tently  excited  by  an  input  signal  which  is  independent  of  the  disturbances. 
Simulations  have  confirmed  this  requirement.  We  therefore  continue  here  to 
use  alternate  periods  of  (1)  open-loop  control  using  a  dither  signal  with  the 
controller  turned  off  and  (2)  closed-loop  control  using  the  controller.  We 
continue  to  use  a  square  wave  input  rudder  command  6^  of  0.5  radian  ampli¬ 
tude  and  0.5  ship  length  period  for  the  dither  signal  in  this  section. 


The  need  to  excite  the  system  can  be  illustrated  further  by  considering 

equations  (85),  (87),  (90),  (94), (102),  and  (114).  The  MVE  parameter  estimator 

eq.  (114)  shows  that  the  estimator  improves  its  estimates  primarily  through 

the  term  K^I^  •  The  estimator  gain  eq.  (102)  shows  that  the  gain  is  directly 
at 

proportional  to  C  .  The  innovation  vector  I  is  shown  in  eqs.  (94)  and  (90) 

/N  ^  ^  ^ 

to  depend  upon  C  ,  D ,  and  the  measurements  z'  .  From  their  definitions 
k-1  k-1  .  K 

A  A 

in  eqs.  (85)  and  (87) ,  ^  and  ^ ,  respectively,  depend  directly  upon  the 

measurements  z^  As  a  result,  if  the  ship  is  closely  controlled  along  the 
desired  path  so  that  the  measurements  tend  to  zero  both  and  1^  tend  to  zero 
and  the  estimator  cannot  improve  its  estimates.  For  the  parameter  estimator 
to  be  effective,  the  measurements  must  contain  non-zero  states  which  are  large 
compared  to  the  measurement  noise  levels.  Periods  of  open-loop  control  using 
a  dither  signal  are  therefore  needed. 


To  complete  the  design  of  the  MVE  parameter  estimator  for  time-invariant 

parameters  for  the  Tokyo  Maru  a  number  of  design  choices  must  be  made  to  define 

the  equivalent  noise  covariance  matrix  R  defined  by  eq.  (101),  define  the 

eq 

generalized  noise  power  matrix  G  defined  by  eq.  (113) ,  and  to  initialize  the 

error  covariance  matrix  which  is  propogated  in  accordance  with  eq.  (103). 

Recall  that  with  time-invariant  parameters,  the  diffusion  coefficient  matrix  4> 

associated  with  the  parameter  model  eq.  (83)  is  zero.  We  used  the  reference 

measurement  noise  level  from  Table  9  to  establish  the  measurement  noise 

covariances  needed  to  calculate  the  matrices  K,B,U,  and  G.  The  discrete 

measurement  noise  standard  deviations  a*,  in  Table  9  were  squared  to  form  the 
2  3 

noise  variances  o needed  in  the  diagonal  measurement  noise  covariance 
matrix  R,  the  B  »dtrix  defined  in  eq. (104),  the  weighted  covariance  matrix 
y  ^=U  defined  in  eq.  (106) ,  and  the  generalized  noise  power  matrix  G^  j=G 
defined  in  eq.  (113).  With  time-invariant  parameters,  the  variances  ol  in 
the  weighted  covariance  matrix  U  are  zero.  To  complete  the  definition  of  U, 
we  chose  the  "mean”  values  of  the  parameters  needed  in  eq.  (106)  to  be  near 
the  midpoint  of  the  permissible  ranges  of  the  parameters  defined  in  Table  10; 
i.e. , 


22  s 

Pl= 

-1.750, 

23! 

P2= 

6.250, 

32S 

P4= 

0.175, 

33 1 

p5= 

-0.785. 

To  complete  the  definition  of  the  equivalent  noise  covariance  matrix 

T 

R  ,  the  disturbance  power  matrix  F'o  P  defined  by  eq.  (105)  must  be  es- 
eq  w 

tablished.  The  coefficients  in  eq.  (105)  were  assigned  assumed,  constant 
values  for  H/T=1.89  from  Table  2.  The  remaining  required  quantities  are 
q'  and  q* ,  the  assumed  spectral  densities  for  the  continuous  white  noise 
disturbances  N'and  Y' ,  respectively.  The  choice  of  these  quantities  has 
a  major  impact  on  the  effectiveness  of  the  MVE  parameter  estimator.  If 
these  quantities  are  selected  to  be  too  large,  the  disturbance  power  matrix 
tends  to  be  "large”  and  the  equivalent  noise  covariance  matrix  R  also 
tends  to  be  "large”.  It  can  be  seen  in  eq.  (108)  that  the  parameter  estimator 
gain  matrix  is  proportional  to  R~*  .  If  R^  is  "large",  the  estimator  gain 
matrix  is  "small"  and  the  parameter  estimator  will  change  its  estimates  very 
slowly.  This  is  true  even  when  the  innovation  vector  I  contains  significant 


information  . 


The  speed  of  response  of  the  parameter  estimator  is,  therefore, 

directly  controlled  by  the  choice  of  and  q^.  We  initially  assumed  these 

quantities  to  have  the  values  given  in  eqs.  (67)  and  (68),  respectively. 

These  values  are  slightly  larger  than  the  spectral  densities  derived  from 

the  passing  ship  disturbance  shown  in  Fig.  3  and  given  in  eqs.  (53)  and  (54) . 

49 

Cuong  has  shown  that  using  these  values  the  MVE  ship  parameter  estimator 
response  is  far  too  slow  to  be  of  any  value  in  an  adaptive  ship  path  control 
application.  We  therefore  chose  these  quantities  to  be  based  upon  approxi¬ 
mately  0.1  times  the  passing  ship  disturbance  shown  in  Fig.  3,  i.e., 

q-  =  2.00  x  10~10,  (115) 

hn 

q'  =  1.25  x  10~9  .  (116) 

hY 

These  assumed  disturbance  spectral  densities  produce  a  MVE  parameter  estimator 
which  responds  at  an  acceptable  rate.  These  assumptions  complete  the  defini¬ 
tion  of  the  disturbance  power  matrix  and  therefore  also  the  equivalent  noise 

covariance  matrix  R  and  the  generalized  noise  power  matrix  G. 
eq 

The  final  quantity  to  be  selected  is  the  initial  estimate  error  covariance 
Pq.  This  assumption  is  not  critical  to  the  effectiveness  of  the  MVE  estimator. 
It  affects  the  initial  rate  of  convergence  somewhat  but  its  effect  dies  out 
fairly  fast.  To  obtain  an  initial  Pq  it  is  reasonable  to  assume  that  all  ini¬ 
tial  parameter  errors  are  independent  Gaussian  random  variables  with  zero 
means.  Hence,  all  off-diagonal  elements  of  Pq  are  zero.  The  diagonal  terms, 
which  are  the  error  variances,  remain  to  be  chosen.  They  can  be  selected 
based  upon  the  expected  range  of  each  parameter.  In  the  simulations,  we 
used  the  following  initial  estimate  error  covariance  which  assumes  that  the 
error  covariances  are  the  square  of  the  difference  between  the  respective  para- 


meter  values  at 

H/T=°°  and 

H/T=l. 

89  shown 

in  Table 

2. 

.03452 

0 

0 

0 

0 

0 

0 

6.640 

0 

0 

0 

0 

P  = 

0 

0 

.02568 

0 

0 

0 

.  (117) 

o 

0 

0 

0 

.02047 

0 

0 

0 

0 

0 

0 

.01185 

0 

0 

0 

0 

0 

0 

.00003091 

— 

To  evaluate  the  performance  of  the  MVE  parameter  estimator  with  time- 
invariant  parameters,  we  simulated  the  Tokyo  Mar u  under  the  control  of  an 
adaptive  path  controller  using  the  MVE  parameter  estimator  defined  above. 

The  simulations  presented  here  represent  a  type  of  startup  condition  for  the 
controller.  The  Kalman  filter  ship  model,  path  controller  gains  C,  Kalman 
filter  gains  K,  and  the  parameter  estimator  were  all  nitialized  for  deep 
water  ship  characteristics.  The  ship  and  Kalman  filter  were  given  zero 
initial  conditions  x=*x=0.  At  the  start  of  the  simulations,  the  depth-to-draft 
ratio  was  abruptly  changed  to  H/T=1.89  or  viewed  differently  the  adaptive 
path  controller  was  turned  on  with  the  depth-to-draft  ratio  H/T=1.89.  The 
depth-to-draft  ratio  was  then  held  constant  to  produce  time-invariant  ship 
characteristics.  The  reference  measurement  noise  from  Table  9  was  utilized. 

The  dither  signal  was  a  square  wave  rudder  command  6'  of  0.5  radian  amplitude 

c 

and  0.5  ship  length  period.  The  dither  signal  was  utilized  only  for  the  first 
2.5  snip  lengths  of  each  6.25  ship  lengths.  This  corresponds  to  an  LW=2.5 
ship  length  data  window  and  LU*6.25  ship  length  update  period  for  the  WLS  para¬ 
meter  estimator.  The  MVE  parameter  estimator  frequency  was  chosen  as  10  sample 
times  or  about  2.4  seconds. 

The  Tokyo  Maru  was  simulated  under  the  control  of  the  MVE  parameter 

estimator  adaptive  path  controller  while  being  subject  to  (1)  no  external 

disturbances,  (2)  the  passing  ship  disturbance  shown  in  Fig.  3  with  the  ship 

passing  beam-to-beam  at  t'=7,  and  (3)  fairly  small,  continuous  white  noise 

disturbances  with  o  *  2x10  5  and  o  =5x10  5.  The  simulation  conditions  and 
N  Y 

the  approximate  ship  path  for  the  simulation  with  the  passing  ship  dis¬ 
turbance  are  summarized  in  Fig.  28.  The  ship  reached  a  maximum  offset  of 
n’=. 01738  or  about  one  beam  at  t'=3.8.  Since  the  simulation  represents  a 
startup  transient  for  the  controller,  the  effect  of  the  subsequent  dither 
periods  is  generally  less.  The  parameter  estimate  results  for  the  first 
10  ship  lengths  of  the  three  simulations  are  summarized  in  Table  16.  since 
there  are  200  estimates  for  each  parameter  in  the  10  ship  lenqths.  Table  16 
presents  only  the  initial  guess  and  the  minimum,  maximum,  mean,  and  final 

values  of  the  estimates  for  each  parameter.  The  RMS  error  in  each  estimate, 
200 

"  RMS(p.)  =  y  <p  -  p,)2  1/2 ,  (118) 

k=l 
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Figure  28.  Simulation  Conditions  and  Approximate  Ship  Path  —  Constant 
Bottom  Example. 


is  also  shown  as  a  basis  for  evaluating  the  effect  of  the  disturbances  on 
the  estimator  performance.  Recall  that  the  control  loop  requires  effective 
estimates  of  only  the  two  most  sensitive  parameters  f23  and  f^,..  T*ie  esti~ 
mates  of  these  two  parameters  and  f  are  good.  The  other  parameter  estimates 
are  slow  and  generally  ineffective.  The  disturbances  do  not  significantly 
degrade  the  performance  and  in  some  cases  the  estimates  improve  due  to  the 
additional  ship  motion. 


i 


i 


Disturbance  Parameter 


Parameter  Estimates  from  MVE  Parameter  Estimator  —  Constant  Bottom  Examples 


The  estimates  of  parameters  f  , f  ,  and  f^  from  the  simulation  with  the 
passing  ship  disturbance  are  shown  in  Figures  29,  30,  and  31,  respectively. 

The  true  values  are  shown  by  constant,  dashed  lines.  The  simulation  results 
with  the  reference  noise  level  from  Table  9  are  shown  as  the  solid  lines.  The 
simulation  results  with  0.001  times  of  the  reference  noise  level  are  sketched 
approximately  as  long,  broken  dashed  lines  in  Figures  29,  30,  and  31  for  refer 
ence.  The  estimate  of  parameter  f  shown  in  Fig.  29  converges  to  near  its 
true  value  in  about  3  ship  lengths.  The  estimate  of  parameter  f  shown  in 
Fig.  30  converges  to  near  its  true  value  in  only  one  ship  length.  After  about 
t'=4  it  varies  about  the  low  noise  estimate  which  appears  to  have  a  small  bias 
error.  These  estimates  appear  to  be  more  than  acceptable  for  use  by  an  adap¬ 
tive  path  controller.  The  estimate  of  parameter  f  shown  in  Fig.  31  is 

slower  to  converge  and  shows  two  jumps  at  the  end  of  each  dither  period.  The 
effect  of  the  dither  signal  can  be  seen  fairly  clearly  in  the  low  noise  esti¬ 
mates  in  Figures  29,  30,  and  31.  The  additional  ship  motion  produced  by  the 
passing  ship  disturbance  which  acts  from  about  6<t'<8  appears  to  help  the 
parameter  estimates  move  closer  to  their  true  value. 


TIME 


true 

value 


Figure  29.  Estimate  of  f23  by  MVE  Parameter  Estimator  —  Constant  Bottom  with 
Passing  Ship  Example. 


Estimate  of  f25  by  MVE  Parameter  Estimator  —  Constant  Bottom  with 
Passing  Ship  Example 


The  simulations  presented  in  this  section  show  that  the  MVE  ship 

parameter  estimator  can  be  effective  in  estimating  the  two  parameters 

needed  by  an  adaptive  path  controller,  f„,and  f„,  when  these  parameters 

2  3  25 

are  constant.  This  performance  does  not  deteriorate  with  significant 
zero-mean  disturbances.  The  performance  with  time-varying  parameters  is 
evaluated  in  the  next  section. 


5.4  Time-Varying  Parameters  Performance 


In  this  section,  we  evaluate  the  performance  of  the  MVE  ship  para¬ 
meter  estimator  when  the  parameters  are  time-varying.  As  in  Section  4.3.1, 
we  use  the  water  depth  as  the  independent  variable  which  defines  the  changing 
parameters.  We  simulated  the  Tokyo  Mar u  under  the  control  of  an  adaptive  path 
controller  using  the  MVE  parameter  estimator  while  sailing  over  two  bottom 
configurations.  The  controller  design  was  identical  to  that  developed  for 
the  Tokyo  Maru  in  the  previous  section  except  that  it  was  extended  to  ac¬ 
commodate  the  time-varying  parameters.  Two  design  changes  are  necessary. 

2 

First,  non-zero  parameter  variances  must  be  utilized  when  the  weighted 

covariance  matrix  U,  ,  is  calculated  in  accordance  with  eq.  (106) .  We 
k-1 

utilized  the  square  of  half  the  range  of  each  parameter  in  Table  2  as  the 


assumed  variance.  Finally,  a  non-zero,  diffusion  coefficient  matrix  <t> 
must  be  included  in  the  parameter  estimator  equation,  eq.  (114),  and  in  the 
estimate  error  covariance  propogation  equation,  eq.  (103).  This  diagonal 
matrix  is  composed  of  the  diffusion  coefficients  for  the  Brownian  motion 
processes  which  were  used  to  model  the  parameters  in  eq.  (83) . 


The  diffusion  coefficient  matrix  41  can  be  established  using  the 

properties  of  the  Brownian  motion  process  developed  in  Section  3.3.  There 

we  showed  that  the  power  spectral  density  of  the  continuous  white  noise  in 

the  Brownian  motion  models  for  the  disturbances  could  be  approximated  by, 

2 

q  =  AtS  ,  (66) 

where  S  is  the  expected  rate  or  change  or  slope  of  the  continuous  process 
being  modeled  and  At  is  the  integration  stepsize  used  in  the  Euler  integra¬ 
tion  discretization  of  the  continuous  process.  Now  in  the  derivation  of  the 
MVE  parameter  estimator  we  modeled  the  parameters  by  discrete  models,  eq. (83), 
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where  the  discrete  time  step  is  the  estimation  period  or  £At.  The  element 

*  if  the  diffusion  coefficient  matrix  is  the  variance  of  the  white  sequence 
ii 

ir.  eq.  (83)  so  using  eq.  (64)  we  have, 

*ii  =  5  =&  ' 

or  using  eq.  (66) ,  (119) 

_  £2 

'ii  i  * 

We  approximated  the  $  using  eq.  (119)  by  assuming  that  the  parameters  could 
diffuse  the  difference  between  their  values  at  H/T  =  1.89  and  H/T  =  00  in 
Table  2  in  the  nondimensional  time  of  about  1.5  ship  length.  The  results 
are  summarized  in  Table  17  for  £=  10. 


difference 

slope  S 

♦ii 

Pl=  f22 

0.196 

0.1314 

1.726  x  10-3 

P2=  f23 

2.716 

1.8218 

3.319  x  10_1 

P3=  f 25 

0.169 

0.1133 

1.284  x  10~3 

P4=  f 32 

0.151 

0.1011 

1.023  x  10-3 

P5=  f 33 

0.114 

0.0770 

5.924  x  10”4 

P6=  f 35 

0.00586 

0.00393 

-6 

1.545  x  10 

_ 

Table  17.  Design  Diffusion  Coefficient  Matrix 


We  simulated  the  Tokyo  Maru  under  control  of  an  adaptive  path  controller 
using  the  MVE  parameter  estimator  while  sailing  over  a  rising  bottom  and  a 
falling  bottom.  These  simulations  were  performed  (1)  without  external  distur¬ 
bances  and  (2)  with  the  passing  ship  disturbance  shown  in  Fig.  3  with  the 
ship  passing  beam-to-beam  at  t'=7.  The  simulations  were  begun  with  zero 
initial  conditions;  i.e.,  x=x=0.  The  dither  signal  had  a  0.5  rad.  magnitude 
and  0.5  ship  length  period.  The  periods  of  open-loop  control  using  the 
dither  signal  and  closed-loop  control  using  the  path  controller  were  the 
same  as  used  in  the  previous  section.  The  reference  measurement  noise  level 
from  Table  9  was  utilized.  For  these  simulations,  the  initial  parameter 
estimates  were  assumed  known  and,  therefore,  the  initial  estimate  error 
covariance  PQ  was  taken  as  zero  in  lieu  of  using  eq.  (117) . 


For  tho  rising  bottom  example,  the  bottom  was  assumed  to  start  at  a 
depth-to-draf t  ratio  H/T=10.0  (corresponding  to  deep  water  characteristics). 
The  depth  then  decreased  one  draft  per  ship  length  until  t'=8.7  at  which 
r.oi'.t.  it  reached  H/T=1.30  and  remained  constant.  This  corresponds  to  an 
upward  bottom  slope  of  5.5  percent.  The  simulation  conditions  and  the 
approximate  ship  path  for  the  simulation  with  the  passing  ship  disturbance 
are  summarized  in  Fig.  32.  Even  though  the  bottom  was  smoothly  rising, 
the  parameter  profiles  were  more  irregular.  The  resulting  parameter  his¬ 
tories  show  two  regimes;  i.e. , 

•  a  slowly  varying  region  from  H/T=10.0  to  2.5; 

•  a  quickly  varying  region  from  H/T=2.5  to  1.30. 

This  example  was,  therefore,  very  demanding  of  the  MVE  parameter  estimator. 

If  the  MVE  is  designed  for  the  more  slowly  varying  parameters,  the  estimator 
is  too  slow  to  track  the  variation  in  the  second  phase.  On  the  other  hand, 
if  it  is  designed  for  large  parameter  variations,  its  estimates  are  more 
sensitive  to  noise  and  random  disturbances.  The  MVE  design  studied  here, 
as  reflected  in  the  parameter  variation  slopes  used  to  obtain  $  in  Table  17 
is  more  appropriate  for  the  parameter  variation  in  the  first  phase  of  this 
example. 

In  the  rising  bottom  example,  the  parameter  estimator  was  initialized 
with  the  correct  deep  water  parameters.  The  simulation  was  performed  for 
15  ship  lengths.  The  results  of  the  first  10  ship  lengths  of  the  two  rising 
bottom  simulations  are  summarized  in  the  upper  half  of  Table  18.  The  table 
shows  the  initial  value,  final  value,  and  true  final  value  for  each  para¬ 
meter.  The  RMS  error  in  each  estimate  calculated  using  eq.  (118)  is  also 
shown  as  a  basis  for  evaluating  the  effect  of  the  passing  ship  disturbance 
on  the  estimator  performance.  Again  recall  that  the  control  loop  requires 
effective  estimates  of  only  the  two  most  sensitive  parameters,  f  and  ^25' 

In  general,  the  presence  of  the  passing  ship  disturbance  helps  to  improve 
the  accuracy  of  most  parameter  estimates.  This  behavior  can  be  explained 
by  noting  that  the  ships  pass  beam-to-beam  at  t’=7  so  the  disturbance, 
occurring  mostly  during  the  duration  of  the  second  dither  signal,  increased 
the  ship  motion.  It,  therefore,  increased  the  information  content  in  the 

innovation  vector.  It  can  be  seen  from  these  simulations  that  f  ,  and  f„_ 

33  35 

are  insensitive  parameters;  their  estimates  tend  to  stay  near  the  initial 

values.  Parameter  f  ,  is  also  fairly  insensitive.  Parameter  f  is  slightly 

3*. 
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more  sensitive.  Parameters  f23  and  f  are  the  most  identifiable  parameters 
as  expected  form  the  sensitivity  study. 


i  par  i  path  i 
J  est.  |  control  ' 

i  >  I 

onnm  mnnr 

dither  signal  :  amp|.=  0.5  rad.  ,  period  =  0.5 


nondimensional  tinne  scale  :  ship  length 


yyT  =  i.3  0 


Figure  32.  Simulation  Conditions  and  Approximate  Ship  Path  —  Rising  Bottom 
Example. 

A  A 

The  paths  of  the  parameter  estimates  f  and  f^  produced  by  the  MVE 
parameter  estimator  are  shown  as  solid  lines  in  Fig.  33  and  Fig.  34,  respec¬ 
tively,  for  all  15  ship  lengths  of  the  rising  bottom  simulation  with  the 
passing  ship  disturbance.  The  correct  parameter  values  are  shown  as  dashed 
lines  on  these  figures.  The  assumed  parameter  variation  slopes  used  to  estab¬ 
lished  the  diffusion  coefficient  matrix  $  in  Table  17  are  also  shown  as  solid 
straight  lines  for  reference.  The  very  rapid  parameter  change  at  about 

A  A 

t ' *7. 5  is  clearly  seen.  The  estimates  and  f25  follow  the  true  parameter 
values  fairly  well  during  the  slowly  varying  portion  of  the  simulation.  The 
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effect  of  the  passing  ship  disturbance  can  be  seen  beginning  at  about  t'=5.5. 
The  passing  ship  disturbance  ends  at  t'=8.5.  The  MVE  is  unable  to  track  the 
rapid  rise  in  parameters  beginning  at  t'=7.5.  It  is  able  to  follow  f  well 
up  to  the  very  rapid  change  which  begins  at  about  t'=8.5.  At  this  point,  the 
parameter  change  becomes  more  rapid  and  the  dither  signal  is  stopped. 

The  estimate  of  f  is  worse  in  this  period.  The  parameter  estimator  was 
still  reacting  to  the  passing  ship  when  the  rapid  change  in  the  parameters 
began  to  occur  and  required  a  few  ship  lengths  to  catch  up.  Recall  that  the 
ship  is  very  stable  at  the  final  water  depth  and,  therefore,  the  dither  sig¬ 
nal  excites  less  ship  motion  and  the  parameters  are  specially  hard  to  estimate 
in  the  final  condition. 


Figure  33.  Estimate  of  f  by  MVE  Parameter  Estimator  —  Rising  Bottom  with 
Passing  Ship  Example 


For  the  falling  bottom  example,  the  water  depth  was  assumed  to  in¬ 
crease  such  that  the  depth-to-draft  ratio  increased  linearly  from  H/T=1.30 
to  h/T-2.50  in  10  ship  lengths.  This  is  a  much  less  drastic  change  in 
ship  characteristics  than  present  in  the  rising  bottom  example.  The  simu¬ 
lation  conditions  and  the  approximate  ship  path  for  the  simulation  with  the 
passing  ship  disturbance  are  summarized  in  Fig.  35.  The  parameter  estimator 
was  initialized  with  the  correct  ship  parameters  and  the  simulation  was  run 
for  the  10  ship  lengths.  The  results  of  the  two  falling  bottom  simulations 
are  summarized  in  the  lower  half  of  Table  18.  In  general,  the  RMS  errors 
are  greater  with  the  passing  ship  disturbance.  The  final  value  of  the 
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Figure  35.  Simulation  Conditions  and  Approximate  Ship  Path  —  Falling 
Bottom  Example. 
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estimate  of  f  is,  however,  closer  the  true  value  with  the  disturbance. 

As  observed  before,  f„,  and  f  are  the  needed  and  the  most  identifiable 

23  25 

parameters. 

A  A 

The  paths  of  the  parameter  estimates  f  and  produced  by  the  MVE 
parameter  estimator  are  shown  as  solid  lines  in  Fig.  36  and  Fig.  37,  respec¬ 
tively,  for  the  falling  bottom  simulation  with  the  passing  ship  disturbance. 

The  correct  parameter  values  are  shown  as  dashed  lines  on  these  figures. 

The  assumed  parameter  variation  slopes  used  to  establish  the  diffusion  coef¬ 
ficient  matrix  4  in  Table  17  are  also  shown  as  solid  straight  lines  for 
reference.  The  estimates  follow  the  true  values  fairly  well.  At  the  begin¬ 
ning  of  the  simulation,  the  ship  is  so  stable  that  the  dither  signal  has 
little  effect.  As  a  result  there  is  little  ship  motion  and  very  little  in¬ 
formation  available  to  the  parameter  estimator.  The  estimates  tend  to  stay 
near  their  initial  values.  After  about  2  to  3  ship  lengths,  the  ship  motion 

A  A 

becomes  large  enough  that  the  estimates  f  and  f^  start  tracking  the  true 
parameter  profiles.  The  effect  of  the  passing  ship  disturbance,  which  is 
felt  from  5.5<t'<8.5,  is  clearly  evident. 

In  general,  these  simulations  show  that  a  MVE  ship  parameter  estimator 
can  be  effective  in  estimating  time-varying  ship  characteristics.  The  dif¬ 
fusion  coefficient  matrix  4  selection  must,  however,  be  consistent  with  the  rate 
of  parameter  change  encountered  by  the  ship.  The  estimator  performance  de¬ 
grades  during  periods  of  external  disturbances,  but  in  the  long  run  the  ship 
motion  caused  by  the  disturbances  can  result  in  improved  parameter  estimates. 

5.5  Parameter  Identif iabil ity  and  Measurement  Importance 

A  useful  by-product  of  the  study  of  the  MVE  ship  parameter  estimator 
is  a  rough  indication  of  the  identif iabil ity  of  each  parameter  and  the  rela¬ 
tive  importance  of  the  various  measurements  to  the  estimation  of  the  unknown 
parameters.  The  numerical  values  of  the  optimal  parameter  gain  matrix 
give  this  information.  The  principal  part  of  the  parameter  estimating  scheme 
is  given  by  eq.  (96);  i.e., 

Ek-  =  £k-£+  V  h" 


where , 


~k’  ”  'VWVk'  ' 


(96) 


At  any  time  step  k' ,  each  element  k  of  the  gain  matrix  K.«  represents  an 

mn  K. 

optimal  weighting  factor  on  the  innovation  i  with  respect  to  the  parameter 

A 

estimate  p  .  Furthermore,  the  innovation  i,  relates  to  the  measurement  of 

V '  * i 2  to  the  measurement  of  r' ,  i^  to  the  measurement  of  B',  and  i^  to  the 

measurement  of  n'.  Hence,  the  magnitudes  of  k  ,,k  .k  ,  and  k  ,  give  an 

ml  m2  m3  m4 

indication  of  the  relative  importance  of  these  four  measurements  to  the  esti 

mate  p  . 
m 

Since  the  various  measurements  have  different  magnitude  scales,  the 
most  revealing  form  of  is  to  "normalize”  each  column  of  the  gain  matrix 

by  the  ■  xisting  magnitude  of  the  associated  measurement.  Likewise,  since 
each  unknown  parameter  has  its  own  magnitude  scale  we  can  further  "normalize 
each  row  of  the  gain  matrix  by  the  existing  magnitude  of  the  associated  para 
meter.  A  typical  result  is  shown  in  Table  19.  In  general,  the  magnitudes 
of  the  rows  of  Table  19  compared  vertically  indicates  the  relative  identi- 
fiability  of  each  of  the  associated  parameters.  Copf f icients  f^>  an<3 

f,^  are  therefore  the  most  easily  estimated  parameters.  The  magnitudes  in 
the  columns  of  Table  19  compared  horizontally  indicate  the  relative  impor¬ 
tance  of  each  of  the  associated  measurements  to  the  estimation  of  each  of 
the  associated  parameters.  This  indicates  roughly  whore  money  should  be 
pi  ‘  » o  improve  the  accuracy  of  the  measurements .  Since  the  control  loop 

•  ■  is  effective  estimates  of  only  f  and  f,r,  an  adaptive  ship  path  con¬ 
trol  let  requires  high  accuracy  and  low  noise  in  the  measurements  of  r’  and 

1 .  it  is  encourag ing  to  note  that  the  measurement  of  the  lateral  offset 

•  1  .s  relatively  unimportant  in  this  context. 


»  ■u..?i.ited  associated  measurement 
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Tat le  19.  Typical  "Normalized"  Optimal  Parameter  list imator  Gain  Matrix 
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5.6  Computational  Requirements 

As  noted  in  the  introduction  to  this  section,  the  MVE  ship  parameter 
estimator  would  require  less  than  100  dynamic  storage  locations  in  the 
implementation.  The  estimator  would  also  require  very  little  CPU  time. 

The  subroutine  used  in  our  work  was  not  specifically  designed  to  minimize 
operations  but  was  timed  to  require  an  average  of  only  12  milliseconds  per 
parameter  estimation  cycle  on  the  University  of  Michigan  Amdahl  470/V7 
computer.  This  very  small  load  on  an  onboard  computer  would  make  the  MVE 
parameter  estimator  an  attractive  candidate  for  the  parameter  estimator 
component  of  an  adaptive  ship  path  control  system.  In  such  arrangement,  the 
MVE  algorithm  could  be  readily  implemented  in  a  time-sharing  mode  on  a  main 
and  general-purpose  computer  onboard  a  ship.  It  could  also  be  programmed 
in  a  dedicated  minicomputer.  This  arrangement  might  be  advantageous  when 
an  independent  and/or  portable  real-time  parameter  estimator  is  desired. 
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6.  Conclusions 


The  research  summarized  in  this  report  has  concerned  the  development  and 
evaluation  of  techniques  which  could  be  utilized  in  adaptive  path  control 
systems  for  surface  ships  in  restricted  waters.  In  this  section,  we  will 
review  some  of  the  principal  conclusions  and  observations  based  on  this 
work . 

•  The  example  vessel  used  in  this  work  was  the  Tokyo  Maru,  a  290m  tanker, 
for  which  sufficient  hydrodynamic  data  was  available  from  the  work  of  Fujino. 
The  characteristics  of  this  ship  undergo  significant  changes  depending  upon 
its  operating  condition  and  environment.  This  type  Of  ship  is  either  course 
stable  or  unstable  depending  upon  water  depth.  Further,  the  longest  time 
constant  of  this  ship  is  very  long  compared  with  the  time  scale  of  importance 
in  path  control  in  restricted  waters.  The  example  vessel,  therefore,  reflects 
many  of  the  most  challenging  requirements  for  an  adaptive  ship  path  control 
system. 

•  The  basic  concept  of  a  two-loop  adaptive  path  control  system  as  shown  in 
Fig.  L  appears  feasible.  The  inner  or  control  loop  performs  the  state 
estimation  and  control  generation  continuously  as  a  high  priority  task.  The 
outer  or  gain  update  loop  can  estimate  the  system  parameters  and  calculate 
updated  controller  and  state  estimator  gains  at  a  slower  rate  or  batchwise. 

In  general,  the  two  loops  must  function  independently. 

•  A  very  effective  control  loop  can  be  designed  using  random-walk,  or  more 
correctly  Brownian  motion,  models  for  the  unknown  yawing  moment  and  lateral 
force  disturbances  which  act  upon  the  ship.  These  controllers  are  effective 
with  short-term,  zero  mean  disturbances  and  with  more  slowly  varying,  non¬ 
zero  mean  disturbances.  These  controllers  have  the  integral  control  property 
necessary  to  accommodate  bias  disturbances  without  a  constant  offset  from  the 
desired  track.  Severe  disturbance  transients  such  as  the  entrance  to  the 
SNAME  H-10  Panel  ABC  harbor  (Fig.  5)  would,  however,  require  some  type  of 
supplemental,  anticipatory  controller.  The  random  walk  disturbance  model 
control  loops  can  be  designed  as  statistical  steady-state  Kalman-Bucy  filters 
and  optimal  state  feedback  controllers  using  Potter's  algorithm  provided  a 
small  numerical  parameter  is  included  in  the  design  equations  to  eliminate 
otherwise  zero  eigenvalues.  Their  performance  is  not  severely  degraded 


by  changing  ship  coefficients.  When  designed  using  incorrect  ship  parameters, 
the  transient  performance  of  the  control  loop  is  degraded  but  its  mean  steady- 
state  condition  is  unchanged. 


•  The  Brownian  motion  process  is  a  useful  concept  which  can  be  used  as  a 

model  for  disturbances  or  processes  which  are  constant  (bias)  or  contain 
unknown,  slowly  varying  components.  The  properties  of  the  Brownian  motion 
process  can  be  used  to  provide  a  p;  actical  design  approach  for  establishing 
the  diffusion  coefficient  of  these  models.  Instead  of  estimating  both  a 
correlation  time  and  a  variance  to  establish  the  power  spectral  density  for 
a  first-order  shaping  filter  model,  the  single  modeling  parameter,  the 
diffusion  coefficient,  can  be  estimated  by  eq.  (65),  (66)  or  (119)  using 

the  expected  rate  of  change  or  diffusion  of  the  process  being  modeled. 

•  An  adaptive  ship  path  controller  needs  effective  on-line  estimation 

of  only  two  of  the  ten  coefficients  of  the  state  equations  which  change  with 
ship  conditions  and  environment.  Adapting  for  changes  in  only  f23  and  f25 
allows  the  recovery  of  the  most  important  part  of  the  transient  response 
performance  which  is  lost  by  not  using  exact  ship  characteristics.  Coeffi¬ 
cient  f23  is  the  coefficient  of  i)’  and  f2j  is  the  coefficient  of  6’,  respec- 
tively,  in  the  r-equation.  This  fact  should  be  recognized  in  the  evaluation 
of  the  effectiveness  of  parameter  estimators  intended  for  use  in  an  adaptive 
ship  path  controller. 

•  The  degree  of  difficulty  of  the  adaptive  control  problem  depends  on  the 
signal -to-noise  ratio.  For  a  surface  ship,  this  ratio  is  very  small  when 
viewed  on  the  time  scale  of  importance  in  restricted  waters.  As  a  result, 
the  adaptive  path  control  of  surface  ships  becomes  a  very  difficult  problem. 

•  An  open-loop  dither  signal  is  required  for  successful  ship  parameter 
estimation.  During  the  estimation  period,  the  control  command  has  to  be 
independent  of  disturbances  and  the  ship  motion  has  to  be  maintained  at  a 
sufficient  amplitude.  Adaptive  path  control  can  be  achieved  by  performing 
the  parameter  estimation  phase  and  the  path  control  phase  alternately.  This 
creates  a  major  tradeoff  in  design  since  the  improved  path  following  perfor¬ 
mance  achieved  by  having  improved  estimates  of  the  ship  coefficients  must 
offset  the  loss  in  path  following  performance  needed  to  allow  the  estimation 
of  the  characteristics.  A  "robust"  ship  path  controller  design  which  is 


-103- 


t 


insensitive  to  errors  in  the  ship  equations  of  motion  would  likely  prove 
superior . 

•  The  degree  of  estimation  accuracy  depends  strongly  on  the  relative 
levels  of  information  and  noise  (process  and  measurement) .  The  amount  of 
information  can  be  increased  significantly  by  lengthening  the  dither  signal 
period  and  the  length  of  time  it  operates.  This,  however,  is  restricted 
due  to  the  path  control  consideration.  It  was  found  that  the  square  wave 
rudder  command  dither  signal  of  0.5  rad.  amplitude  and  0.5  ship  length 
period  operating  for  over  2  ship  lengths  gave  reasonable  results  with  the 
linear  ship  model  studied  here.  Design  tradeoff  studies  would  be  needed 

to  establish  the  best  tradeoff.  In  general,  we  used  low  levels  of  measurement 
r.cise  out  demanded  more  parameter  estimate  accuracy  than  really  needed  by 
the  control  loop. 

•  In  very  shallow  water,  the  Tokyo  Maru  becomes  so  stable  that  the  dither 
signal  used  here  has  a  relatively  small  effect  on  the  ship.  As  a  result, 
there  is  less  ship  motion  and  a  lower  signal-to-noise  ratio.  Parameter 
estimation  is,  therefore,  much  more  difficult  in  this  situation. 

•  The  parameters  f 22 >  ^2  3'  and  ^25  r-e(luation  are  closely  associated 

with  the  shorter  time  constant  of  the  ship  which  is  about  0.4  ship  length. 

The  parameters  f32r  £33*  and  £35  in  the  (1-equation  are  more  closely  associated 
with  the  longest  time  constant  of  the  ship  which  is  about  10  ship  lengths. 

As  a  result,  a  short  (0.5  ship  length)  dither  period  can  be  effective  in 
excitinq  the  ship  in  the  shorter  time  constant  mode  allowing  reasonable 
estimation  of  parameters  f22»  f 2 3 '  and  ^25-  The  short  dither  period  is  much 
less  effective  in  exciting  the  longer  time  constant  mode  making  estimation 
of  parameters  f32>  ^33>  and  f  35  less  effective.  Longer  dither  periods  which 
would  be  needed  to  more  effectively  excite  the  longer  time  constant  mode  would 
be  incompatible  with  a  path  control  objective.  Fortunately,  the  control  loop 
only  really  needs  accurate  estimates  of  f23  and  fj^  to  achieve  an  effective 
level  of  adaptation. 

•  The  choice  of  parameters  and  the  availability  of  measurements  are  of 
critical  importance  to  the  Weighted  Least-Squares  (WLS)  parameter  estimator. 
Given  a  data  window  length  acceptable  in  a  path  control  application,  the  WLS 
cost  function  is  extremely  sensitive  to  the  measurement  noise  and  the?  system 
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model.  We  found  that  all  states  must  be  measured  and  that  all  six  of  the 
coefficients  fpj  in  the  open-loop  dynamics  matrix  and  control  distribution 
matrix  in  eq.  (71)  must  be  included  in  the  parameter  vector.  Since  the  cost 
function  is  much  more  sensitive  to  external  disturbances  than  ship  coefficient 
changes,  the  unknown  external  d  1st  ur  buncos  must  also  be  included  in  the  para¬ 
meter  vector.  The  WLS  algorithm  could  not  effectively  handle  linear  combina¬ 
tions  of  parameters  or  products  of  parameters  so  it  was  necessary  to  assign 
the  coefficients  y^j  in  the  disturbance  distribution  matrix  assumed,  constant 
values . 

•  With  a  2.5  ship  length  data  window,  the  WLS  parameter  estimator  can 
effectively  estimate  the  coefficients  fpj  when  the  ship  characteristics  and 
external  disturbances  are  essentially  constant  during  the  data  window  loading 
period.  A  major  design  tradeoff  exists  in  the  choice  of  the  data  window 
length.  With  a  longer  data  window  and  therefore  a  longer  period  of  open-loop 
dither  signal  operation,  the  information  content  of  the  data  base  is  increased 
allowing  better  parameter  estimate  accuracy  but  the  ship  is  given  a  greater 
offset  from  the  desired  path.  The  improved  control  loop  performance  achieved 
with  more  accurate  parameter  estimates  must  offset  the  loss  in  path  following 
performance  caused  by  the  dither  signal.  With  our  experience  to  date,  we 
must  conclude  that  the  performance  of  the  control  loop  designs  studied  in 
Chapter  3  are  good  enough  with  incorrect  ship  parameters  that  we  would  not 
expect  an  adaptive  capability  using  a  WLS  parameter  estimator  to  prove 
worthwhile . 

•  The  WLS  parameter  estimator  could  not  effectively  estimate  the  coeffi¬ 
cients  f^j  when  either  the  ship  characteristics  or  the  external  disturbances 
changed  significantly  during  the  data  window  loading  period.  In  path  control 
situations  in  restricted  waters,  both  would  be  expected  to  change  during  at 
2.5  ship  length  data  window.  The  WLS  parameter  estimator  studied  here  could, 
however,  be  of  value  in  other  applications  where  constant  ship  characteristics 
and  disturbances  could  be  reasonably  expected. 

•  The  Doesn't  Use  Derivatives  (DUD)  algorithm  used  here  to  minimize  the 
WLS  cost  function  proved  to  be  an  effective  and  efficient  search  algorithm. 

In  retrospect,  however,  it  requires  an  excessive  amount  of  dynamic  data 
storage.  Thus,  it  is  more  suited  for  an  off-line  application  than  for  an 
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on-line  application  as  part  of  an  adaptive  system.  Tt  did  prove  effective 
for  the  purposes  of  this  study.  Alternative  search  algorithms  with  more 
reasonable  storage  requirements  are  available  for  on-line  applications.  The 
stopping  conditions  of  any  algorithm  would  have  to  be  tuned  to  be  compatible 
with  the  noise  content  of  the  data  window. 

•  The  Minimum  Variance  (MVE)  parameter  estimator  can  effectively  estimate 

the  time-varying  coefficients  f23  and  f25  needed  by  an  adaptive  path  controlle 
when  the  ship  is  not  subjected  to  bias  disturbances.  The  algorithm  is  re¬ 

cursive  and  highly  efficient  both  in  dynamic  data  storage  and  computation 
time.  The  MVE  algorithm  requires  ship  motion  to  obtain  information  about 

the  ship's  characteristics  so  the  tradeoff  associated  with  the  use  of  an 
open-loop  dither  signal  exists  with  this  method  as  it  would  with  any  alter¬ 
native.  The  method  does  show  promise  and  is  worthy  of  further  consideration 
for  possible  use  in  an  adaptive  ship  path  controller. 

•  The  principal  weakness  of  the  MVE  parameter  estimator  as  studied  here  is 
it;:  inability  to  accommodate  bias  disturbances.  We  initially  tried  to  incor¬ 
porate  consideration  of  these  disturbances  by  modeling  the  disturbances  as 
crovr.ian  motion  processes  in  lieu  of  white  noise  in  the  derivation  of  the 
disturbance  power  matrix  r'Qwr'  .  This  approach  produces  a  disturbance 
power  matrix  which  grows  continuously  with  time  and  this  causes  the  parameter 
estimator  gain  matrix  to  asymptotically  approach  zero.  An  alternative,  which 
could  not  be  included  in  the  scope  of  this  study,  would  be  to  follow  an 
approach  parallel  to  that  used  with  the  WL.S  parameter  estimator.  The  unknown 
yawing  moment  and  lateral  force  disturbances  could  be  included  as  additional 
parameters;  i.e.  the  parameter  vector  could  be  chosen  as  eq.  (74)  in  lieu  of 
oq.  (78)  and  the  entire  MVE  algorithm  could  be  rederived.  We  now  feel  this 
should  be  successful  and  intend  to  study  this  possibility  in  future  work. 

•  The  three  major  design  choices  in  the  development  of  the  MVE  algorithm 
are  the  estimation  frequency  C,  the  assumed  disturbance  power  spectral 
densities  q^  and  q^.,  and  the  parameter  model  diffusion  coefficient  matrix  i. 
The  derivation  requires  that  the  parameters  not  be  estimated  more  frequently 
than  every  two  sample  times;  i.e.  £> 2.  Our  work  showed  little  value  in 
estimating  the  ship  parameters  more  often  than  £=10  or  about  2 . 4  seconds .  An 
even  longer  estimation  period  would  probably  be  acceptable.  The  choice  of 


the  disturbance  power  spectral  densities  directly  affects  the  speed  of 
response  of  the  MVE  alqorithm.  If  they  are  assumed  too  high,  the  MVE 
parameter  estimator  is  much  too  slow  to  respond  to  be  of  value  in  a  path 
control  application.  The  parameter  model  diffusion  coefficient  matrix  4> 
cm  be  selected  using  the  properties  of  the  Brownian  motion  process  as 
reflected  in  eq.  (119).  This  choice  also  affects  the  response  rate  of  the 
MVE  parameter  estimator.  If  the  <J>^£  are  too  small,  the  estimator  cannot 
follow  more  rapidly  changing  parameters.  If  the  <pa  are  too  large,  the 
estimator  will  be  overly  sensitive  to  noise  and  disturbances. 

•  The  consideration  of  the  optimal  MVE  estimator  gain  matrices  shows  that 
in  order  to  effectively  estimate  the  coefficients  f23  and  f25  needed  by  an 
adaptive  path  controller  the  measurements  r  and  8  are  the  most  important. 

The  lateral  offset  from  the  path,  which  would  be  the  most  difficult  to 
treasure  accurately,  is  of  much  less  importance.  The  development  of  an 
adaptive  ship  controller  should,  therefore,  emphasize  accuracy  and  noise 
reduction  in  the  measurements  of  yaw  rate  and  lateral  velocity. 
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